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Abstract 

The non-relativistic energy levels of ortho-positronium are calculated in the 
quadrupole and octupole approximations for the interaction potential. For this 
O ' purpose, the RST eigenvalue problem of angular momentum is illustratively solved 

for the quantum numbers jo = 0, 1, 2, 3, 4 and ^^^jz = ±1- This eigenvalue problem 
admits ambiguous solutions for < \^^^jz\ < jo whereas the solutions for ^^'^jz = 
and ^^^jz = ijo are unique. In order to attain some (at least approximative) so- 
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I lutions of the energy eigenvalue problem one tries a factorized ansatz for the wave 

■ function and thus splits off the angular problem (with its ambiguous solutions) from 

I the residual radial problem. The latter does, as usual, finally fix the energy eigen- 

values. But it is just by this procedure that the ambiguity of the angular problem 

^ I is transferred to most of the energy levels which thereby become doubled. The 

H ■ 

_Cy _' corresponding doubling energy amounts to (roughly) one percent of the total bind- 

ing energy and is, however, of purely electric origin, since magnetism is completely 
neglected. Indeed, the charge distributions of both positronium constituents (i.e. 
electron and positron) do inherit their ambiguity from the ambiguous solution of 
the angular eigenvalue problem (-^ charge "dimorphism"); and naturally the di- 
morphic configurations must then possess slightly different interaction energies of 
the electrostatic type. ^ 



PACS Numbers: 03.65.Pm - Relativistic Wave Equations; 03.65.Ge - 
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Theories and Applications; 03.75.B - Matter Waves 
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I Introduction and Survey of Results 

The present paper aims at the further elaboration of Relativistic Schrddinger The- 
ory (RST), i.e. a fluid-dynamic version of relativistic few-particle quantum theory [1,2]. 
This fluid-dynamic approach is applied here in order to study the non-relativistic energy 
spectrum of ortho-positronium; and it is found that there does occur a somewhat strange 
phenomenon, i. e. the doubling of the number of the energy levels. But since such a fluid- 
dynamic view on the physics of elementary particles is rather unusual, it may perhaps 
appear desirable to first premise some conceptual and epistemic remarks about its rela- 
tion to the conventional quantum theory which is strongly dominated by the probabilistic 
point-particle picture of quantum matter. 

As desirable and instructive as such a comparison of our fluid-dynamic concepts to 
the conventional point-particle logic may be, there unfortunately arises a certain difficulty 
here because the point-particle proponents seem still to disagree about their fundamental 
notions. 

The birth of quantum mechanics, some nine decades ago, has first been accompanied 
by considerable confusion and many controversies about the epistemic status of the new 
theory, see ref. [3] for a historic account of the notorious Bohr-Einstein debate. But now 
that almost a century has passed away, and a considerable technical and cultural progress 
could be made thanks to the overwhelming applications of the new theory (think, e. g., of 
nuclear and laser technology, semiconductors, materials science, supernovae astrophysics, 
big bang cosmology etc.) one should nowadays think that the teething troubles of the 
theory's childhood have been successfully overcome already long time ago. Indeed, it is 
hardly conceivable how such an immense civilizing progress could have been made without 
possessing a correspondingly profound theoretical basis. In this sense, one expects that 
present-day quantum theory has been anchored in the mean time on solid fundaments, 
comparable to the situation with classical electromagnetism or thermodynamics, and thus 
will be able to provide us with hopeful future developments. 
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However, this impression is grossly fallacious. Indeed, a superficial inspection of some 
modern presentations of the fundamentals of conventional quantum theory is sufficient 
in order to ascertain that the controversies about the true meaning of this theory and 
its conceptional framework are becoming now rather more violent than dying out. This 
regrettable fact is brought to light mostly on the occasions of celebrating the very birth of 
quantum mechanics, e.g. Heisenberg's matrix mechanics of 1925. Usually, the inventor (or 
discoverer) of a new fruitful theory (such as, e.g., Galilei [4], Newton [5] or Maxwell [6]) 
enjoys later on much reputation by the scientific descendants; and the corresponding first 
texts about the original new idea are frequently celebrated by historic reprints. However, 
concerning Heisenberg's matrix mechanics of 1925, the situation is rather ambiguous: 
whereas some people do hold it adequate to celebrate Heisenberg's invention (discovery) 
by dedicating a book to his original proposal of matrix mechanics [7], others argue that 
this proposal "consisted essentially in the introduction of novel algorithms (i. e. his matrix 
mechanics) and even renowned and formalistic physicists do nowadays confirm that this 
paper is actually unintelligible" (ref.[8], p. 39; authors' translation). 

A further subject of never-ending controversies concerns even the central concept of 
quantum theory: i.e. the wave function (or more generally: quantum state). Is the wave 
function to be associated to some real property of the observed system and thus does 
exist also when the system is not observed? Or does the wave function merely reflect the 
state of information of the observer? In the latter case, the wave function is thought to 
"collapse" whenever the observer's state of information undergoes a sudden change (e. g. 
by reading off the pointer of some measurement device). A competent opinion looks as 
follows: "In the first place, it is often argued that the wave function itself should not be 
regarded as giving an objective description of the world (or part of it) but as providing 
information merely of 'one's state of knowledge' about the world. This view I really cannot 
accept." (R. Penrose in ref. [9], p. 121). Another opinion is the following: "The claim, 
that the wave function be not real but describe 'information', reminds me of arguments 
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which are widespread in homeopathy" (ref. [8], p. 79; authors' translation). 

Some co-authors do also frankly concede that they disagree even about the fundamen- 
tal quantum concepts: "Although we have tried to write a 'coherent' book, our reader will 
soon notice that our conceptions vary on some basic notions. Characteristic are our differ- 
ent opinions on the relevance of mathematical concepts for the interpretation of quantum 
mechanics and hence different inclinations to make active use of these concepts in physical 
arguments." [10] 

It should be evident that such a loose bedrock is not suited to erect on it a continual 
scientific and technological progress; and it seems that already the fathers of the new 
theory (especially Bohr and Heisenberg) were aware of such a grave drawback of their new 
theory and therefore they tried to cobble together some minimalistic interpretation of their 
quantum formalism (i.e. the kopenhagen interpretation, see, e.g., ref. [9]) which should 
suffice to handle with the quantum formalism "for all practical purposes" . However, the 
Bohr-Heisenberg interpretative proposal of their formalism evokes nowadays some harsh 
critique: "// a physicist declares that he consistently applies the kopenhagen interpretation, 
then this does mean nothing else than that he consistently always refers to the same kind 
of conceptual inconsistency - but this does not yield a consistent theory^'' (ref.[8], p. 51; 
authors' translation). 

Wave-Particle Duality 

One of the key concepts of the kopenhagen interpretation does refer to the comple- 
mentarity principle, with the wave-particle duality representing the most conspicuous 
exemplification thereof. Indeed, the progress of the past decades was made by extensive 
use of both the particle and the wave concept in like manner. Both concepts have been 
thought to be unrenouncable for the description of elementary matter and its interactions. 
This fact becomes immediately evident by merely glimpsing at the titles of some of the 
competent text books, see ref.s [ll]-[25]. On the other hand, there seem to arise now seri- 
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ous arguments which point to the possibihty that 'particles' (in the classical and quantum 
sense) do not really exist but merely are some kind of fiction (if not illusion) : " There are 
no quantum jumps, nor are there particles" [8,26] . That might be true. But even if the 
long-lived concept of a 'particle' (classical and quantum) ultimately should turn out as 
a mere fiction, then it is surely a very useful fiction since it advanced the progress of 
the past decades considerably. If now a new time should dawn with the particle concept 
being superseded by some novel, more powerful concept (perhaps decoherence [8,10]?, or 
emergence [27]?) then that novel concept will also possess the epistemic status of being 
not more than a fiction\ Indeed, it is hold quite generally that the limited epistemic 
abilities of a human being admits to describe 'reality' at most in terms of fictions [28] 
which, however, must then of course be required to be free of logical and observational 
contradictions. 

* 

Thus it seems to us that the notorious quantum controversies are nothing else than 
some kind of useless struggles for domination among different fictions. In place of opposing 
these different fictions against each other in an irreconcilable way, one should rather con- 
sider them as helpful complementations of each other, so that they together can mediate 
to us a more integral picture of what is really going on in the quantum world. 

Bearing this situation in mind, one may doubt whether the ardently desired and unique 
'theory of everything^ [29] will ever be discovered (or invented, resp.) in order to supersede 
Bohr's complementarity principle. So long as we must be content with theories having a 
loose end, we should rather recall Bohr's viewpoint on this question to our mind. Indeed, 
it seems that Bohr was the first who noticed clearly that the behaviour of elementary par- 
ticles could not be satisfactorily described in terms of one single logical framework. More 
concretely, he thought that the (self-suggesting) probabilistic point-particle description 
of an elementary particle should be complemented by a fluid-dynamic wave description, 
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namely in order to better understand and manage the notorious wave-particle duality. 
This original idea of Bohr is now known as the 'complementarity principle^ [9] (see also 
the corresponding remarks in the precedent paper [1]). According to this principle (if it 
is understood in the right sense) one needs two logical systems (or "pictures" ) , excluding 
one another to a certain degree, in order to become able to predict and 'understand' the 
outcomes of the considered experiments. 

For instance, in the well-known two-slit interference experiment for electrons one be- 
comes forced to find some plausible explanation why the electrons do form the well-known 
interference pattern on the detection screen after having passed the two-slit arrangement. 
The solution of this problem is based upon the postulate that the electrons do pass the 
two-slit region in form of a wave (with appropriate wave length) ; but their arrival on the 
detection screen occurs in form of point-like particles [30] . 

In order to manage intelectually such a highly ambiguous phenomenon it seems wise to 
us to resort to Bohr's complementarity principle and develop two complementary logical 
systems (or 'pictures'), i.e. the probabilistic point-particle picture and the fluid-dynamic 
wave picture. This means that we have to specify two mathematical formalisms together 
with their associated physical interpretations so that we can 'understand' and predict the 
outcomes of physical experiments and/or observations. Whether or not a future theory of 
everything will be able to incorporate simultaneously both pictures (or create a completely 
new picture) must be left unclarified for the time being. 

Probabilistic Point- Particle Picture 

The necessity for developing additionally also a fluid-dynamic picture of the quantum 
objects becomes perhaps most clear when one reconsiders the somewhat artificial manner 
in which the point-particle theory tries to explain the undeniable wave-like aspects which 
become evident in the notorious double-slit experiment. 

A closer inspection of the interference pattern and the associated double-slit geometry 
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in terms of the particle picture may suggest the following conclusion: The confining of the 
particle position to a spatial volume measured by a certain linear dimension (slit distance 
Ax, say) must induce in some (mysterious) way a momentum uncertainty {Ap^, say) so 
that both point-particle quantities will obey the Heisenberg uncertainty relation 

Ax-Ap,>^, (I.l) 

being needed for explaining the appearance of the interference pattern (see any textbook 
about elementary quantum mechanics which treats the double-slit experiment in the point- 
particle picture, e.g. [31]). Furthermore, since the interference pattern on the detection 
screen is built up by a huge number of particle-like impacts, one thought that one must 
describe such a phenomenon in terms of a statistical theory. For such a purpose, the 
appropriate mathematical formalism appeared to be the Hilbert space of quantum states 
1^ > with the classical observables (such as position x and momentum p) being transcribed 
to the corresponding operators (here x and p) acting over that Hilbert space. The particles 
passing through the double slit and striking at the detection screen do then build up 
a statistical ensemble which is to be characterized by the statistical operator (density 
matrix) p, with tr p = 1. Thus, let the quantum state of a point particle, located point- 
like at ^ (0 ^ ^ Ax), be denoted by the Hilbert space vector |^> 

^|e>=eie>, (1-2) 

then the statistical operator p for such an ensemble of particles with unknown (but well 
determined) locations in the double-sht region (0 < ^ < Ax) is adopted as 

Aa; 

p^j dC^Olexel- (1-3) 



Here, the probabihty w{^) for the space point C (0 ^ C ^ ^^) being occupied by some 
particle is to be normalized to unity as usual 

Ax 

jdiw{i)^l. (L4) 
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Finally, it remains to demonstrate how the desired Heisenberg uncertainty relation 
(II. ip can arise from such a mathematical construction. For this purpose, one merely has 
to interprete the position and momentum uncertainties Ax and Ap as the corresponding 
mean square deviations ( "variances" ) 

{Axf = tr [p (x- < X >Y] (I.5a) 
{Apf = tr [p (p,- < >)'] , (I.5b) 

and additionally one postulates the Heisenberg commutation relations for the operators 
X and p as follows: 

[x,pj =ih-l . (1.6) 

But the crucial point is now that this probability construction fll.2p - fll.6p cannot bring 
forth that desired uncertainty relation fll.ip which one would like to see being validated 
in order to become able to explain the appearance of the interference pattern. Here, 
the problematic constituent of the whole probability construction is the assumption (11. 3p 
which means that anyone of the point particles does occupy a unique (albeit unknown) 
point of space. Namely, such a collection of pure point particles does not develop the 
right lateral pressure for generating a transverse momentum px across their direction of 
flight (y-axis, say) which is necessary in order to produce just the observed interference 
pattern on the detection screen. Indeed, if this assumption fll.Sp is rejected and replaced 
by the assumption that any "point particle" is in a (exotic) state > which is some 
^^superposition" of the states |^ > (11. 2^ . being themselves due to the classical localization 
of a real point particle, then it is no problem to deduce the wanted uncertainty relation 
fll.ip . Namely, the statistical operator p becomes now 

p=|^><^| (1.7) 
with the superposition |\l/> of the classical configurations |^ > being defined through 

Ax 
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Obviously, this is just the place where the point-particle picture has to introduce fluid- 
dynamic concepts, since the particle's spatial presence is (in some mysterious way) now 
smeared out over some finite region. And consequently, the variances (11.5ap - (]1.5b|) emerge 
now in the position representation as 

Ax 

(A.)^^/...rW(.-<.>fW.) (I.9a) 



Ax 

{Ap^f = Jdx rix) (p.- <p,>fiP{x) , (I.9b) 



where the mean values <x>, <Px> are defined as usual through 

Ax 

<x> = J dx ip* {x) X ip{x) (1. 10a) 



Ax 

<Px> =F J dx ip*^^) Pxtp{x) (1. 10b) 



h d 

Px = -TT- 

l OX 

Now, whenever the considered system (here: a point particle confined to the one- 
dimensional interval [0, Aa;]) is in such a superposition |\1' > of classical configurations 
|.^> and therefore can be described by a "wave function" iplx) (=<x|\I^>), then a quite 
general theorem ensures that an inequality of the following type must hold (Schwarz' 
inequality, e. g. [32]) 
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(Ax)'-(Ap)^>(^|<C>|) , (1.11) 



where < C > is the expectation value of the commutator C 

C =[x,p,]=zhl . (1.12) 



Obviously, the result (II.lip - flI.12p for the variances Ax, Apx is now just the wanted 



uncertainty relation fll.ip . 
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It appears natural that such a successful invention of a mathematical formalism (to- 
gether with the physical interpretation of the inherent mathematical quantities) will be 
celebrated as a great progress in ^understanding^ the working in the micro- world. So much 
the more the formalism can be extended in order to deal with all possible physical situa- 
tions. It is merely necessary to set up a Hilbert space whose basis vectors may be selected 
in such a way that they describe the classically realizable configurations (|^>, say); but 
other basis systems are also possible. The Hilbert space for the many-particle systems 
is the tensor product of the one-particle Hilbert spaces. Next, one writes down the gen- 
eral state vector |\l/>, if desired as a superposition of the classical states |^ >, cf. (II. 8|) : 
and then one lets this state vector |\E' > evolve in time according to the time-dependent 
Schrodinger equation 

ih^^\^>t=H\^>t . (1.13) 
For the stationary states one puts 

|^>t= e"^^*/'' |^>o (1.14) 

which then yields the time-independent Schrodinger equation 

H\'^>o= E\'^>o (1.15) 

which, e. g., admits to determine the energy spectrum {E} of bound systems. This quanti- 
zation formalism may not only be applied to mechanical systems but also to field systems 
and then entails the well-known quantum field theory. But despite its overwhelming suc- 
cess the latter quantum-formalism continues to own the epistemic status of a probabilistic 
point-particle fiction; and therefore the question must be asked whether there do perhaps 
exist other 'fictions' which, after thorough elaboration, do also yield a comprehensive and 
consistent picture of the micro-world? 
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Complementarity and Environment 



As already mentioned at various occasions, the precedent point-particle picture seems 
to us to be plagued by a certain deficiency: i.e. the very concept of a "point" particle. 
In the original classical sense, a point particle does occupy some point x of 3-space and 
possesses a definite momentum p. Concerning the precise position, the corresponding state 
vector 1^ > in the one-particle Hilbert-space obeys the eigenvalue equation fll.2|) . However, 
by means of the set (11. 3 P of such state vectors one cannot explain the interference pattern 
on the detection screen! Therefore one brings into play the conception of a coherent 
superposition |\E' > f ll.8|) of all classically possible positions ^ (0 < ^ < Ax) so that this 
range Ax of all those possible locations together with the range of all possible momenta Ap 
can obey the required uncertainty relation (ll.ip . It seems that such a strange superposition 
of classical locations has been considered because one felt being forced to take account of 
the particle's finite extension" within the logical framework of point particles! Indeed, 
this auxiliary construction may be understood in the sense that the point particle becomes 
somehow "smeared" over the whole interval Aa; and thus appears as a kind of droplet of 
some esoteric fluid; whereas on other occasions it may be treated as a true point particle! 
It surely does not come as a surprise that such a self-contradictory and mysterious idea 
of a point particle has provoked the question: "/s this now profound thoughtfulness or the 
ultimate degree of craziness?" (ref.[8], p. 52; authors' translation). 

Most textbook authors prefer to evade the notion of a point particle with unknow- 
able position and momentum and they rather try to explain the double-slit experiment 
completely in terms of the wave picture [30]. But then the problem becomes urgent anew 
when one resolves microscopically the interference pattern, being thought to arise in the 
wave picture by the superposition of different wave trains. However, one finds the pattern 
being composed of point-like impacts on the detection screen. So one feels oneself being 
thrown back again to the notorious phenomenon of wave-particle duality. 

In such a confused situation it seems reasonable to suppose that the dominance of 
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either the particle aspect or the wave aspect is brought forth cooperatively by the object 
itself plus its surroundings. Such a holistic theory (if feasible) could then interpolate 
between the pure particle-like and the pure wave-like aspects of the considered physical 
situation. But when one cuts the system off its environment (=^ closed system) the 
balance of the particle-like and wave-like aspect goes lost and one has to decide whether 
one wants to neglect the particle aspect or the wave aspect. In many situations, one of 
both aspects will be the dominant one and therefore suggests itself to become preferred 
over the other one for the description of the isolated system. In other instances, both 
aspects may be mutually counterbalanced; and if one nevertheless cuts the system off its 
surroundings and favours only one of the complementary aspects, then there will arise 
those notorious "quantum paradoxes" , e. g. the mentioned double-slit phenomenon. In 
recent time, there were put forth some proposals of treating the system in combination 
with its surroundings (see, e.g., ref.s [8,10,27]) but presently a final judgement seems 
premature. 

Historically, the method of neglecting the system's environment, together with the 
preference of the probabilistic point-particle picture, has been the prevailing world view 
and the corresponding success is impressive, see ref.s [ll]-[25]. But, if concentrating on 
the system alone by neglecting its surroundings, one should not despise the possibility of 
describing the isolated system also in fluid-dynamic terms according to the complementary 
wave picture. One cannot exclude that the more comprehensive view on "system plus 
environment" can also be built up by starting from the side of the fluid-dynamic approach. 

A proposal of the latter kind is the "Relativistic Schrddinger Theory" [1,2], i.e. a 
fluid-dynamic version of relativistic few-particle quantum mechanics. The subsequent 
investigation is an application of this theory to ortho-positronium to be conceived as a 
closed system [para-positronium has been the subject of the preceding paper [2]). More 
concretely, the present paper is a study of the question to what extent the non-relativistic 
RST spectrum of ortho-positronium does coincide with that of the conventional quantum 
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theory (which is a probabihstic point-particle theory). The main result of this study is 
that most of the conventional energy levels become duplicated as an indirect consequence 
of the specific way in which the angular momenta of the constituents (i. e. electron and 
positron) are to be composed in RST to the total angular momentum. Unfortunately, 
the corresponding energy eigenvalue problem is in RST much more complicated than in 
the conventional quantum theory so that we have to resort to appropriate perturbation 
techniques which we develop up to the "octupole" approximation [33] . Thus our obtained 
result of level duplication emerges here in the quadrupole and octupole approximations 
and therefore might eventually be an artefact of those perturbation orders, i.e. possibly not 
being present in the exact RST solution of the eigenvalue problem for ortho-positronium. 
However, the main effect of level doubling occurs in the quadrupole approximation (Fig. 
VILA on p IlOSp and the relative doubling energy amounts to more than 1%; on the other 
hand, the octupole correction of the levels in the quadrupole approximation amounts to 
less than 1 percent (table on p. 11451) . This hints at level doubling being a true effect in RST 
which receives merely some minor corrections from the higher-mult ipole approximations. 

* 

It is true, the origin of this level duplication can be traced back to the fact that 
ortho-positronium has unity of spin (-w sa = 1); whereas para-positronium has zero spin 
(-^ s-p = 0) and therefore does not undergo that phenomenon of level duplication. But it is 
important to note that the energy difference of both levels (emerging by duplication) is not 
of magnetic but rather of electric origin! Namely, for the ortho-spin so = 1, the eigenvalue 
problem for angular momentum admits (mostly) two different solutions belonging to the 
same values of total angular momentum ja and its 2;-component This then yields two 
different angular distributions of the corresponding charge densities (-w dimorphism, see 
Fig. IV. B on p. 1751) which thereby acquire different interaction energies of the electrostatic 
type. 
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Subsequently, these results will be elaborated in detail along the following arrangement. 



Fluid-Dynamic Character of RST 

As proponents of a fluid-dynamic description of the elementary matter we tend to the 
hypothesis that 'matter' (in its most general sense) is always spread out over some finite 
region of three-space and can never be concentrated in a truly point-like manner. The 
fiction of a point-like particle appears to us as a (more or less realistic) idealization; and 
consequently the finite-size effects of matter must be somehow simulated in the point- 
particle picture, preferably by saying that the state |\E'> of the point-particle refers to a 
"superposition" of states being due to really point-like positions |^>, see equation (jl.81) . 
By this construction it may seem now that the concept of wavefunction iIj{x) (=<x|\E'>) 
is the crucial point in the probabilistic approach and thus must be made responsible for 
its overwhelming success. 

But actually RST is also based essentially on the use of wavefunctions ^I/(x), namely 
in order to generate the physical densities of the considered system such as, e.g., the 
total current density or the energy-momentum density T^^ of the considered system, 
see equations (]II.29|) and (1II.35|) below. Thus, in contrast to the situation with the point- 



particle case, the RST wave function owns a well-defined meaning, namely to generate 
all the physical densities of the considered system (Sect.s II and III are taken over from 
the precedent paper [2] in order to elucidate the common basis of both para- and or- 
thopositronium -w ortho/para dichotomy). The RST handling with the wave functions 
\E'(x) is quite different from the probabilistic point-particle case when many-particle sys- 
tems are considered: here, the conventional theory relies on the tensor product of the 
one-particle Hilbert spaces and therefore sets up the dynamics f lI.13|) in the configuration 
space. However, RST adopts the Whitney sum of the one-particle fibre bundles 

^(x) = i)i{x) © i)2{x) © i)^{x) © . . . (1.16) 
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which then admits to set up the dynamics in the real four- dimensional space-time, namely 
in form of the Relativistic Schrddinger Equation for the wave function in combination 
with the (non-abehan) Maxwell equations for the bundle connection A^{x), see equations 
dlLTj) and dlLTSj) below. 

This formal difference of the dynamical equations of both approaches entails an impor- 
tant physical consequence: whereas (in the non-relativistic approximation) the probabilis- 
tic conventional approach takes the predetermined Coulomp potential as the basis for the 
interaction mechanism, the corresponding RST mechanism retains the interaction field as 
a truly dynamical object being equipped with its own field equation. For instance, for 
the positronium system (to be considered subsequently) the non-relativistic Hamiltonian 
H (ll.isp contains the non-dynamical Coulomb potential in a manifest way 

fp' e^ 

On the other hand, the corresponding RST interaction potential (^*^Ao(r^)), say, must 
obey the Poisson equation (IIV.19P below, where the source is given by the electric charge 
density (^^^kQ{r)^ generated by the RST wave function Through this arrangement the 
RST eigenvalue problem becomes a system of coupled differential equations, which (in the 
spherically symmetric approximation) consists of the Schrodinger equation f lIV.17|) and 



the Poisson equation fllV.lSp . Naturally, such a coupled eigenvalue system is much more 
difficult to solve than the simple conventional problem ( lI.lSp which for the positronium 
Hamiltonian ( II.17P admits to determine exactly the energy spectrum E^^^ in terms of the 



principal quantum number n: 

Ortho- Dimorphism 

In contrast to this simple situation in the conventional probabilistic theory, the solution 
of the corresponding RST eigenvalue problem can be worked out only approximately. The 
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difficulty refers here to the fact that the wave function becomes considerably anisotropic, 
according to the value of angular momentum ^^''jz', as an example see Fig. IV. A below. But 
since the wave function generates the charge density which then acts as the source of the 
interaction potential, the anisotropy of the wave function is transferred to the interaction 
potential via the Poisson equation so that we ultimately have to solve the (RST form 
of the) Schrodinger equation for an anisotropic potential. Clearly, this can be attained 
only in an approximative way; and most of the paper is concerned with setting up an 
adequate approximation procedure (Sect. VIII) for managing this anisotropy effect. The 
key point is here a suitable factorization of the wave function in a radial part and an 
angular-dependent part which then also becomes transferred to the charge density, cf. its 
product form firvTOjl - firOTj) below. 

As usual, the (approximate) product form of the wave function splits up the energy 
eigenvalue problem in two subproblems, namely (i) the radial problem which depends 
upon the specific physical system to be considered and which ultimately yields the energy 
eigenvalues, and (ii) the eigenvalue problem of angular momentum which is of quite 
general nature and thus is independent of the details of the considered system. Now 
it is just this latter problem, cf. (lIV.6ap - (]IV.6bp . which provides us with the origin of 



the ortho-dimorphism. Namely, for the same values of the angular-momentum quantum 
numbers {jo, ^^^jz} with (—jo) < '■^^jz < jo there do exist two solutions of the eigenvalue 
problem. This ambiguity then is transferred to the angular part ^^^ko^i!)) of the charge 
density, see Fig. IV. B below. 

Induced Energy Difference of Dimorphic Partners 

Since the charge density is the source of the electrostatic interaction potential, cf. 
the Poisson equation (111.43^ . the angular ambiguity of the charge density is immediately 
transferred to the interaction potential which thereby inherits the anisotropy of the charge 
density. More precisely, for any pair of angular- momentum quantum numbers {jo, ^''^jz} 
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the interaction potential ^''^Ao(f) acquires a corresponding angular dependence; and the 
elaboration of this interconnection necessitates an extensive search for the most rational 
representation of the angular dependency of the interaction potential (Sect. VII— VIII). 
There are various possibilities, but ultimately it turns out that the separative method 
provides the most pleasant representation of the wanted multipole expansion because 
here the energy contributions of the various multipole modes can be clearly separated. 
Consequently, one can subdivide the total anisotropy energy into the set of contributions 
of any mode (-w quadrupole energy, octupole energy, . . . ). 

But if once the anisotropy energy is determined, one can substitute this (together 
with the isotropy energy and the kinetic energy) in the RST energy functional and try to 
extremalize this by use of an appropriate trial ansatz. For our present purpose, we are sat- 
isfied with a very simple ansatz for the radial part Q{r) of the wave amplitude, cf. (IVI.lap - 
flVI.lbl) . with only two trial parameters (/3,i^). The value of the RST energy functional 
flVI.25p on the chosen trial configuration yields the corresponding energy function 
, z/) as a function of both trial parameters (3 and u, cf. (1VI.36|) . After minimaliza- 
tion of this function with respect to the first trial parameter /3 we are left with the problem 
of determining the minimal value of the reduced energy function E^^(z/) (1VIII.106P or, 
resp., the maximal value of the associated spectral function Sq^u) (IVIII.lOSp . This 
function contains both quantum numbers jo and ^^'^jz (^ ^''^rrii, ^^^1112) of angular momen- 
tum while the principal quantum number no is restricted to no = jo + 1 on account 
of our too simple ansatz fIVI.lap - flVI.lbp . Thus, one finally obtains the energy spectrum 
of ortho-positronium simply by looking for the maximal value of the spectral function 
u) (IVIILlOSp . namely by admitting all possible quantum numbers jo = 0, 1, 2, 3, . . . 



and ^''^jz, with —jo < ^^^jz ^ jo- This may be performed by means of some appropriate 
numerical program, and the most important results are the following (Fig. VII. A and 
table on p. I145p : 

The energy spectrum resembles the conventional one fll.181) but its (n^)-fold degeneracy 
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becomes eliminated. Such a result has already been found for the spectrum of para- 
positronium [2]; but for ortho-positronium there additionally occurs now a doubling of 
most of the energy levels. The origin of this strange effect traces back to the specific 
spin composition for ortho-positronium {sq = 1), where two different solutions of the 
eigenvalue problem for angular momentum are mostly possible for fixed values of quantum 
numbers jo and ^^^jz- The corresponding energy difference (due to level splitting) amounts 
to (roughly) 1% of the binding energy, see the table on p ll45j it is not caused by magnetism 
but is a purely electric effect due to the ambiguity of the electric charge distribution, see 

Fig.IV.B, p.m 
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II Positronium Eigenvalue Problem 



In order that the paper be sufficiently self-contained, it may appear useful to mention 
briefly some fundamental facts about RST. As its very notation says, the central idea is 
the Relativistic Schrodinger Equation which leads one in a rather straight-forward 

way to the Dirac equation ( ]II.13|) for few-particle (or many-particle) systems. In order to 
ultimately end up with a closed dynamical system for the fluid-dynamic quantum matter, 
one adds the (generally non-Abelian) Maxwell equations flII.15|) where this coupled system 
of matter and gauge field dynamics automatically entails certain conservation laws, such 
as those for charge flll.25p or energy-momentum flll.32p . This fundamental structure of 
RST is then subsequently specialized down to the non-relativistic positronium system, 
especially to its ortho-form, with the main interest aiming at its energy spectrum. 

1. Relativistic Schrodinger Equation 

A subset of problems within the general framework of RST concerns the (stationary) 
bound systems. The simplest of those systems is positronium which consists of two 
oppositely charged particles of the same rest mass (M). The physical behaviour of its 
matter subsystem is assumed here to obey the Relativistic Schrodinger Equation 

thcv^^ = n^^ (ii.i) 

where the two-particle wave function \E'(a;) is the direct sum of the two one-particle wave 
functions ^ipaix) (a = 1,2) 

^(x) = V^i(x) ©V^2(x) . (II.2) 

The gauge-covariant derivative V on the left-hand side of the basic wave equation flll.ip 
is defined in terms of the u(2)-valued gauge potential as usual 

V^^ = d,^ + A,^ , {11.3) 
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or rewritten in component form 

D^^,=d^iP,-iAliJi (II.4a) 
D^^2 = ^,^^J2-^Al^2. (n.4b) 

Here, the electromagnetic four-potentials (a = 1, 2) are the components of the original 
gauge potential Afj, with respect to some suitable basis (a = 1, ... 4) of the u(2)-algebra 

A^{x) = A-(a;) = + B^{x) x - B;{x) x • (11.5) 

Here, the electromagnetic generators (a = 1, 2) do commute 

[ri,r2] = (II.6) 

and the exchange potential -B^ is put to zero B^{x) = 0) because the two positronium 
constituents (i. e. electron and positron) do count as non-identical particles. Recall 
that the exchange effects, being mediated by the exchange potential B^{x), do occur 
exclusively for identical particles so that B^{x) is inactive for the positronium constituents 
(-^ B^{x) = 0). Thus the bundle connection Afj,{x) flll.Sp becomes reduced to its u(l) © 
u(l) projection 

A,ix) A;:(a;)r,. (II.7) 
2. Dirac Equation 

For Dirac particles, which are to be described by four-spinors ipaix), the Hamiltonian 
in the Relativistic Schrodinger Equation (lll.ip obeys the relation 

T^'n^ = Mc\ (H.8) 

where is the total velocity operator and thus is the direct sum of the Dirac matrices 

r^ = (-7^)©7^. (H.9) 
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The mass operator A4 specifies the two particle masses (a = 1, 2) 

M = iM''Ta (11.10) 

and is required to be Hermitian {A4 = Ai) and covariantly constant 

V^M = 0. (11.11) 

This requirement is trivially satisfied for particles of identical rest masses (M^ = = M) 
since for such a situation the mass operator becomes proportional to the identity operator 

M = M 1 . (11.12) 



Thus the result is that, by virtue of the relation (1II.8|) . the Relativistic Schrodinger Equa- 



tion f lll.l|) becomes the two-particle Dirac equation 

ihcT'^V^'^ = Mc"^^ , (11.13) 

or in component form 

ihc 7^ z^i = -Mc^ iJi (II. 14a) 

ihc Y Df, %l)2 = Mc^ ^2 , (II. 14b) 

where the gauge-covariant derivatives (D) of the single-particle wave functions ipaix) are 
given by equations (ITL4a])-(ITL4b]). 



For identical particles, the Dirac equations flIL14ap - (III.14bp would couple both par- 



ticles much more directly since the exchange potential is generated cooperatively by 
both particles and simultaneously does act back on any individual particle which then en- 
tails the phenomenon of self-coupling. However, for the present situation of non-identical 
particles the coupling is more indirect: any particle does generate a Dirac four-current 
ka^i{x) {a = 1, 2) which is the source of the four-potential A'^{x) (see below). And then this 
four-potential of the a-th particle acts on the wave-function ipb{x) of the other particle 
{b 7^ a) as shown by equations flII.14a|) - flII.14b|) in connection with the gauge-covariant 



derivatives D f llL4al) -f llL4b|l 
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3. Maxwell Equations 

The bundle connection A^{x) fill. 51) is itself a dynamical object of the theory (just as 
is the wave function \E'(x)) and therefore must be required to obey some field equation. 
This is the (generally non-Abelian) Maxwell equation 



V^F^,, = -A'Kia,X (11.15) 

he 

Here, the bundle curvature is defined in terms of the bundle connection as usual, 
i.e. 

J'^u = V^A - + [A, A] . (11.16) 

For the present situation of non-identical particles, the connection becomes reduced 
to its (Abelian) u(l) © u(l) projection, cf. ( III.7|) . which then also holds for its curvature 



Decomposing here both the curvature J-'^^, and current operator J'^j^ with respect to the 
chosen basis of commuting generators (a = 1, 2) 

J't^u F;^Ta (11.18a) 
Ju i3\ra , (11.18b) 



one obtains the Maxwell equations flII.15p in component form as 



V'^F;, = 47ra,j,V (11.19) 

Since for the present Abelian situation the curvature components F'^^ (111.18ap are linked 
to the connection components (111.7^ as usual in Maxwellian electrodynamics (in its 



Abelian form) 

= - V^Al , (11.20) 
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the Maxwell equations f lII.19|) for the field strengths become converted to the 



d'Alembert equations for the four-potentials 

□ /l;: = 47^asJ^ (11.21) 
provided the gauge potentials do obey the Lorentz gauge condition 

W^A"^ = . (11.22) 
4- Conservation Laws 

One of the most striking features in the description of physical systems is that both 
classical and quantum matter do obey certain conservation laws. For the presently con- 
sidered Relativistic Schrodinger Theory, as a fluid-dynamic theory, this means that there 
should exist certain local conservation laws, preferably concerning charge and energy- 
momentum. Moreover, these local laws should turn out as an immediate consequence of 
the basic dynamical equations, i. e. the Relativistic Schrodinger Equation f lll.ip and the 
Maxwell equations (III. 150 . 

In this regard, a very satisfying feature of the Relativistic Schrodinger Theory is now 
that such conservation laws are automatically implied by the dynamical equations them- 
selves. In order to elaborate this briefiy, consider first the conservation of total charge 
which as a local law reads 

V% = . (11.23) 

But such a continuity equation for the total four-current can easily be deduced from 
both the matter equation flll.ip and the gauge field equations flll.lSP : and this fact signals 
the internal consistency of the RST dynamics. First, consider the gauge field dynamics 
(III.ISP and observe here the identity 

V^'VT^^ = (11.24) 
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which holds in any flat space-time. Obviously, the combination of this identity with the 

Maxwell equations f lII.15|) yields the following continuity equation in operator form 

V^'J^ = . (11.25) 

Decomposing here the current operator in component form yields 

D^^r, ^ (11.26) 

which furthermore simplifies to 

= (11.27) 

(a = 1,2) 



under the Abelian reduction (lII.18ap -( 1II.18b|) . But when the individual Maxwell currents 



do obey such a continuity equation (]II.27|) . then the total current 



2 



J, = E -^M ("-28) 

a=l 

must also obey a continuity equation which is just the requirement (111.23p . 

On the other hand, we can start also from the matter dynamics (III.1|) and can define 
the total current by 

J,, = ^T^^ . (11.29) 

The divergence of this current is 

V'^j^ = (V^) + ^T^ (P^^) + ^ (VT^) ^ . (11.30) 

Here, one requires that the gauge- covariant derivative of the total velocity operator is 
covariantly constant 

V^'T^ = , (11.31) 

and furthermore one evokes the Relativistic Schrodinger equation together with the Hamil- 
tonian condition f lll.Sp which then ultimately yields again the desired continuity equation 
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(1II.23I) . Thus, the conservation of total charge is actually deducible from both subdy- 
namics of the whole RST system; and this fact supports the mutual compatibility of both 
subdynamics (i.e. the matter dynamics (lll.ip and the gauge field dynamics (III.ISP ). 

A further important conservation law does refer to the energy-momentum content of 
the considered physical system. Aiming again at a local law, one may think of a continuity 
equation of the following form 

(^>r^, = , (11.32) 

where ^'^^^^ is the total energy-momentum density, i. e. the sum of the Dirac matter part 
and the gauge field part (^^'^^^^): 

= (^>r,. + . (11.33) 



Clearly, the validity of the total law (111.32^ does not require an analogous law for the 



subdensities but merely requires the right balance of the energy-momentum exchange 
between the subsystems, i. e. 

(^>r^. = -v^ (^>r^. . (11.34) 

Indeed, the matter part has been identified as 

(^>r^, = ^ [^T^{V,^) - (P.^)r^^ + ^T,{V^^) - {V^^)T,'^] (11.35) 

and the gauge field part by 

= ^K^, {Fl,Fi ^ - lg,.Fl,F^^>^'j , (11.36) 

where K^p is the fibre metric in the associated Lie algebra bundle. The (local) conserva- 
tion law (1II.32P comes now actually about through the mutual annihilation ( III.34P of the 
sources of both energy-momentum densities, i. e. 

(^>r^. = - (^>r^. = he . (11.37) 
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Obviously the sources of the partial densities ^^^^u and '■'^•^T^ty are just the well-known 
Lorentz forces in non-Abelian form. 

It should now appear self-suggesting that the definition of the total energy (-Et) of 
an RST field configuration is to be based upon the time component '■^-Too of the energy- 
momentum density ^"^^^u, i- e. 

Et = J ci¥(^>roo(f) . (11.38) 

But since the total density ^'^^^u is the sum of a matter part and a gauge field part, cf. 
flll.33p . the total energy E-^ 0II.38P naturally breaks up in an analogous way 

= + , (11-39) 

with the self-evident definitions 

^D = y" rf¥(^>roo(f) (II.40a) 

Eg= I d^f^^'^ooi^. (n.40b) 



But clearly, such a preference of the time component Too among all the other components 
T^j, entails the selection of a special time axis for the space-time manifold. This then 
induces a similar space-time splitting of all the other objects in the theory, i. e. we have 
to consider now stationary field configurations which are generally thought to represent 
the basis of the energy spectra of the bound systems. 

5. Stationary Field Configurations 

In the present context, the notion of stationarity is coined with regard to the time- 
independence of the physical observables of the theory, i.e. the physical densities and 
the electromagnetic fields generated by them. In contrast to this, the wave functions 
do not count as observables and therefore are not required to be time- independent. But 
their time- dependence must be in such a way that the associated densities become truly 
t ime- independent . 
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Gauge-Field Subsystem 



The simplest space-time splitting refers to the four-potentials A""^, which for the sta- 
tionary states become time-independent and thus appear in the following form: 

A'^^ix) ^ {("Uo(f) ; -la(f)} (11.41) 
(a = 1,2). 

A similar arrangement does apply also to the Maxwell four-currents j"'^ 

f, ; -Ur)} , (11.42) 



so that the d'Alambert equations fill. 211) become split up into the Poisson equations, for 
both the scalar potentials ''"^Ao(r) 

A^'^^Aoif) = -A7raJ%{r) (11.43) 

and the three- vector potentials Aair) 

AAaif) = -AnaJaif) ■ (11.44) 

Recall here that the standard solutions of these equations are formally given by 

("Uo(rO = as /" d^r' ^^^^^Kl^ (II.45a) 

^a(r-) = as y d'f' ■ (n.45b) 

Similar arguments would apply also to the exchange potential 5^ = ^Bq, but since 
we are dealing with non-identical particles the exchange potential must be put to zero 

It is true, the particle interactions are organized here via the (electromagnetic and 
exchange) potentials which, according to the principle of minimal coupling, are entering 
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the covariant derivatives D^ipa of the wave functions ipa as shown by equations (]II.4aP 



(111.4bp . But nevertheless it is very instructive to ghmpse also at the field strengths F"" . 



Their space-time splitting is given by 

E,= {('^)E^} = {nJ (II.46a) 

Ha = {^''^H^] = \^-e'\F\'^ (II.46b) 

and thus the linear Maxwell equations 0II.19P do split up in three-vector form (a = 1,2) 
to the scalar equations for the electric fields 

V'K = ^T^a,''% (11.47) 

and to the curl equations for the magnetic fields 

Vx#, = 47rasJa. (11-48) 

There is a pleasant consistency check for these linearized (but still relativistic) field equa- 
tions in three-vector form; namely one may first link the field-strengths to the potentials 
in three- vector notation (cf. flll.20p for the corresponding relativistic link): 

E^if) = -V("Uo(f) (II.49a) 
i?,(f) = V X i„(f) (II.49b) 

and then one substitutes these three-vector field strengths into their source and curl 
equations flll.47p - flll.48p . In this way one actually recovers the Poisson equations flll.43p - 



flll.44p for the electromagnetic potentials ^""^Aq, Aa- 

Matter Subsystem 

Concerning now the stationary form of the matter dynamics, one resorts of course to 
the generally used factorization of the wave functions ipa{'f^,t) into a time and a space 
factor 



'ipa{,r,t) = exp 



.M,c' 



2 1 



V'a(r) . (11.50) 
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Here, the mass eigenvalues Ma {a = 1, 2) are the proper objects to be determined from the 
mass eigenvalue equations which we readily put forward now. For this purpose, observe 
first that the Dirac four-spinors ipa may be conceived as the direct sum of Pauli two-spinors 



Mr) = ^%+{r) © • (11.51) 



Consequently, the Dirac mass eigenvalue equations (to be deduced from the general Dirac 
equations (lII.14ap - (]II.14bp by means of the factorization ansatz (1II.50P ) are recast to their 
equivalent Pauli form for the two-spinors ^"•^Lp±{r) 

ta'V%±{f) + (2Uo(f) ■ %^{r) = ~ c ■ %^{r) (II.52a) 
ia-V%±{f) + «Ao(f) ■ %^{r) = -Ml±^c- (2V^(f) . (IL52b) 

(Observe here that we do neglect for the moment the magnetic effects by putting the 
three- vector potentials Aa{r) f lll.4ip to zero: Aa{r) =^ 0). The mass eigenvalue for the 
positron (with rest mass Mp) is denoted by Mi and for the electron (with rest mass Me) 
by Ma. 

Summarizing, the RST eigenvalue system consists of the mass eigenvalue equations 
(III.52ap - flII.52bp for the Pauli spinors ^"'^Lp±(f) in combination with the Poisson equations 



(III.43p . Since the magnetic effects are neglected, the Poisson equations (III.44P need not be 



considered here. But what is necessary in order to close the whole eigenvalue problem is 
the prescription for the link of the Pauli spinors ^"'^(p±{r) to the Maxwell charge densities 
^"'^jo{r), or more generally to the Maxwellian four-currents j""^ f lll.42p as the sources of 
the four-potentials A""^, cf. the d'Alembert equations (III.2ip . Surely, such a link between 
the wave functions ipa and the currents j"^ will have something to do with the Dirac 
four-currents ka^ which are usually defined by (a = 1, 2) 

ka,, = Ipalu'ipa ■ (11.53) 
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Indeed, a more profound scrutiny reveals the following link |3l] 

= = V^i 1u i^i (II. 54a) 

j'^ = -h^ = -^2 7u V'2 • (II. 54b) 

The change in sign of both Dirac currents reflects the positive and negative charge of both 
particles. In terms of the Pauli spinors ^"'^Lp±(f) (a = 1, 2) the space and time components 
of the Dirac currents read 

(r) = ("VV (r) (r) + ^"Vt (r) ^'^V- (r) (II.55a) 
(f ) = (^) ^ ^" V- (^) + VI (r) <? V- (0 • (II. 55b) 

Although the magnetic effects, which originate from the three-currents fca(^) via the mag- 
netic Poisson equations fill. 440 . are neglected in the present paper these current densities 
nevertheless play now an important part for identifying two essentially different kinds of 
positronium. 
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Ill Ortho/Para Dichotomy 



In the conventional theory, the manifestation of two principally different kinds of 
positronium is traced back to the two possibilities of combining the spins of the electron 
(e) and positron (p); if both spins Se and Sp add up to the total Spin S = 1 

So = Sp + Se= 1 (III.l) 

_ 1\ 
2 / ' 

one has ortho-positronium; and zero spin 

Sp = - Se = (III.2) 

yields para-positronium. This is the well known ortho/para dichotomy mentioned in any 
textbook on relativistic quantum mechanics. In contrast to this (generally valid) com- 
position rule for angular momenta, the ortho/para dichotomy in RST is based upon the 
(anti) parallelity of the Maxwellian three-currents jai'f^)- Here, it is assumed that both 
particles do occupy physically equivalent one-particle states ipa{r) (a = 1,2) in the sense 
that the Dirac currents (and therefore also the Maxwell currents) are either parallel or 
antiparallel. Thus we propose the following characterization of ortho-positronium [51] : 

fci(f) = -hir) = h{r) (III.3a) 

Ji(^ = Mr) = jb{r) = h{r) (111.3b) 

A^{r) = A^i^r) = A,{r) (III.3c) 

Hi{r} = H2{f) = Mr) , (III.3d) 
and analogously for para-positronium: 

ki{r) = hir) = kpif) (III.4a) 

Ji(r) = -J2(r) = Jp(r-) = kp{r) (III.4b) 

A(r) = -A2(f) = Ap{r) (III.4c) 

Hi{r) = -Mf) = Hp{^ . (III.4d) 
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The interesting point with such a subdivision of positronium into two classes is the fact 
that this subdivision is based upon the magnetic effects which however are neglected for 
the present paper; but despite this neglection the subdivision is of great relevance also for 
the presently considered electrostatic approximation! Namely, even within the framework 
of the latter approximation scheme, there do emerge different distributions of electrostatic 
charge (^"■'fco(r)) ^r ortho- and para-positronium with corresponding quantum numbers 
(-W ortho/para dichotomy); and moreover there does arise also a certain ambiguity of the 
electric charge distribution even within the subclass of the ortho-configurations due to 
the same quantum number (-w ortho- dimorphism) . In order to elaborate these effects it 
is necessary to first specify the general eigenvalue problem down to the subcases and then 
to look for the corresponding solutions. 

1. Mass Eigenvalue Equations 

The hypothesis of physically equivalent states for both positronium constituents entails 
that both mass eigenvalues Ma are actually identical, i. e. 

Ml = -M2 = -M, . (III.5) 

Next, one considers the Maxwellian charge densities ^""^joir) which must of course differ 
in sign for oppositely charged particles 

= (ni.6) 

and this must be true for both ortho- and para-positronium. On the other hand, the link 
(111.54ap - (111.54bp of the Maxwellian currents j° to the Dirac currents ka^ shows that the 



requirement (1III.6P demands the identity of the Dirac densities ^°'^ko{f), i.e. 

(%ir) ^ (%{r) , (III.7) 
or rewritten in terms of the Pauli spinors flll.55a|) 



%l{f) %+{f) + %l{f) = %i{r} + %Ur) %-{f) . (111.8) 
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However, the situation is different for the Dirac three-currents ka{r), since they may 
differ in sign, cf. ( 1III.3aP vs. (lIII.4a|) . Expressing the (anti) parallehty of both Dirac 



currents in terms of the Pauh spinors, cf. (]II.55bp . one requires 

(^VU^ + %l{r)a'^^^ip+{r) (III.9) 



where the upper/lower sign refers to the ortho/para case, resp. Both conditions OIILSp 
and fllll.Dp can be satisfied by putting 



(iV^(f) = %+{r) = ('/pV+(^) (III. 10a) 

(f) = i %_ if) = (''/^V- if) (III. 10b) 

where the upper/lower sign refers again to ortho/para-positronium, resp. 

It is true, the disposal flIII.10al) - flIII.10bl) satisfies both algebraic requirements fllll.Sp 



and fllll.Qp . but additionally there must be satisfied also a differential requirement: ac- 



tually, the spinor identifications f llll.lOap - fllll.lObp leave us with just one spinor field 



(i.e. ^^V±(^) for ortho-positronium and for para-positronium); and therefore 

both spinor equations (111.52ap - (111.52b|) must collapse without contradiction to only one 



spinor equation (either for '•''V±(^ or *'^V±(0)- In order to validate this requirement, 
we have to make a disposal also for the electrostatic gauge potentials *^"''ylo(r*). But this 
can easily be done by observing the link between the Dirac densities ^"'^^{f) and '■"''y4o(r) 
as it is implemented by the Poisson equations flll.43p . Indeed, this link entails that the 



potentials '^°''>AQ{f) must differ (or not) in sign when this is (or is not) the case also for 
the Maxwell densities ^"^jo(0- Therefore one concludes that for the positronium situation 
both electrostatic potentials must always differ in sign 

(lUo(r) = -('Uo(f) = (^/PUo(r) . (III.ll) 

But when this circumstance is duly respected, both mass eigenvalue equations ( ]II.52ap - 
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(]II.52bl) actually do collapse to a single one for ortho-positronium: 



ta-V ^%±{f) - (%(r) ■ ^%^ir) = c ■ %^ir) . (III. 12) 



The existence of one and the same eigenvalue equation (1III.12|) for both ortho-positronium 



constituents thus validates our original hypothesis that both the electron and the positron 
should occupy physically equivalent states. 

The relativistic pair (1111.12^ of Pauli equations has a single Schrodinger-like equation 
as its non-relativistic limit. Indeed, assuming that the "negative" Pauli-spinor ^^^(f^(f) 
is always considerably smaller than its "positive" companion ^^^{p^{r) one can solve the 
upper one of the equations flIII.12p for ^^'>ip^ (r) approximately in the following form: 

^V(r1^^<?-V(V(r1, (III. 13) 

and if this is substituted into the lower equation (1III.12I) one finally ends up with the 
following non-relativistic eigenvalue equation of the Pauli form 



2M 



A ^%4f) - (%(f) • = ■ ^\+{r) . (III.14) 



Here, the non-relativistic eigenvalue emerges as the difference of the rest mass M and 
the relativistic mass eigenvalue M^,, i.e. 

E, = {M,-M)c\ (III. 15) 

Subsequently, we will be satisfied with clarifying the phenomenon of the ortho-dimorphism 
in the non-relativistic version (1111.14^ of the original relativistic eigenvalue equation 0111.12^ . 



But the case of para-positronium is a little bit more complicated. To begin with 
the positron equation flll.52ap . this becomes transcribed by the identifications fllll.lOap - 
(jlll.lObp . lower case, to the following form (para-positronium): 

i a • V (^V±(^ - ^''^^o(r) ■ ^^Vt(^ = ^*^^ c ■ ^''VtI^ • (HI- 16) 
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Obviously, this positron equation is just of the same form as the joint positron/electron 
equation (1111.12^ for ortho-positronium. However, the electron equation (]II.52b|) of para- 
positronium becomes transcribed by the identifications (lIII.10ap - (lIII.10bl) to a somewhat 
different form: 

ia-V (pV±(^) + ^"^Aoif) ■ (^Vt(^ = - ^*^^ c ■ ^P^if^ir) . (III.17) 

Since the sign of the potential term (~ *^^^y4o(r)) is reversed here in comparison to the 
positron equation (1111.16^ . the latter electron equation (1III.17P is the ''^charge conjugated'^ 
form of the first equation (lIII.lGp . 

The charge conjugation is defined here by the following replacements: 



'(^_(f) , ^P^ifi^ir) (111.18a) 



(pUo(r) -(^Uo(f) (111.18b) 
M, -M, . (III. 18c) 

Indeed, one is easily convinced that the two forms of eigenvalue equations flIII.16p and 
f llll.lTp for para-positronium are transcribed to one another by these replacements flIII.18ap - 
f lIII.18cp . This means that any solution of the positron equation flIII.16p can be interpreted 
also to be a solution of the electron equation (]III.17P : however, the charge conjugation 
does not leave invariant the Poisson equations, cf. f lIII.2ip - (]III.22p below. Therefore, we 
prefer here the use of solutions with the same non-relativistic limit! Indeed, it follows 
from both equations (IIII.16I) and (IIII.17I) that the "negative" Pauli spinor (r) can be 
approximately expressed in terms of the "positive" spinor (P^(y9_|_(r) through 



ih 



y^-[r) ~ ±— — — cr 



■2Mc" ^™-19) 
where the upper case refers to (1III.16P and the lower case to (lIII.lTp . This result may 
then be substituted in either residual equation (IIII.16P and (jlll.lTp which in both cases 
yields the same Schrodinger-like equation for the "positive" spinor ^^''ip^{r): 

- ^ A (^'V+(r) - hc^^^Aoif) ■ (^V+(r) = ■ (PV+(f) . (III.20) 
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Thus it becomes again evident that the corresponding solutions of ( 1III.16I) and (]III.17|) 



do actually describe physically equivalent states. The notion of "physical equivalence" is 
meant here to refer in first line to the numerical identity of the energy being carried by 
anyone of the constituents of para-positronium, see below. 

2. Poisson Equations 

The mass eigenvalue equations do not yet represent a closed system and therefore can- 
not be solved before an equation for the interaction potential ^^/^''Ao(r) has been specified. 
On principle, this has already been done in form of equation (1II.43P so that we merely have 
to further specify that equation in agreement with the ortho/para dichotomy. Observing 
here the circumstance that the Maxwellian current of the first particle (a = 1, positron) 
agrees with the Dirac current, cf. flll.54a|l and flll.55a|l . the Poisson equation fill. 43 11 reads 



in terms of the Pauli spinors 

A (^/PUo(f) = -Ana, '^^'P^oi.r) (III.21) 

This Poisson equation closes the relativistic eigenvalue systems, both for ortho- 
positronium (1111.12^ and for para-positronium f lIII.16l) - flIII.17l) . For the non-relativistic 
limit, one merely suppresses the "negative" Pauli spinors (*/p)(y9_(r) so that the relativistic 
Poisson equations flIII.2ip simplify to 

A (^/^Uo(f) = -47ras ^''/^VU^ ^''^^V+(^ (III.22) 



which then closes both the non-relativistic eigenvalue equations ( 1III.14|) for ortho- 
positronium and (lIII.20p for para-positronium. 

Clearly, these coupled systems of eigenvalue and Poisson equations cannot be solved 
exactly (though exact solutions do surely exist), and consequently we have to resort to 
some adequate approximation procedure. But this suggests itself when we subsequently 
will establish the variational principle of minimal energy. 
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3. Non- Unique Spinor Fields 



Despite the fact that we originally subdivided the whole set of positronium configura- 
tions into two subclasses, i.e. ortho- and para-positronium (lIII.3ap - (lIII.4d|) . it may seem 



now that by the neglection of the magnetic forces we ended up with an eigenvalue problem, 
which does no longer offer any handle for sticking to that original subdivision into two 
peculiar subsets. Indeed, the adopted electrostatic approximation does admit exclusively 
an interaction force of the purely electric type (being described by the electric potential 
(''/p)Ao{r)), whereas the original ortho/para dichotomy was based upon the three- currents 
ka{r) as the curls of the magnetic fields, cf. (1II.48|) . As a result of this neglection of mag- 
netism, the Poisson equation (IIII.2ip . or (IIII.22p . resp., holds equally well for both the 
ortho-configurations (b) and the para-configurations (p). But also the mass eigenvalue 
equations, especially in their non-relativistic forms flIII.14p and flIII.20p are formally the 



same for ortho- and para-positronium. Does this mean that, through passing over to the 
electrostatic approximation, the difference between ortho- and para-positronium has gone 
lost? This is actually not the case because the difference of total angular momentum (in 
combination with the hypothesis of the physical equivalence of both constituent states) 
leaves its footprint also on the electrostatic approximation. 

The crucial point here refers to the fluid-dynamic character of RST, as opposed to the 
probabilistic character of the conventional quantum theory. This entails that in RST the 
angular momenta of the subsystems cannot be combined (to the total angular momentum 
of the whole system) in such a way as it is the case in the conventional theory {-^ addition 
theorem for angular momenta). More concretely: if we wish to insist on the viewpoint 
that the (observable) angular momentum of the considered two-particle system should 
emerge as the eigenvalue jz of the angular momentum operator J'z = + Sz, i. e. 

J, ^,/p{r) = ^'Mj, h ■ m„^{r) , (III.23) 

then the quantum number '-^^''^jz due to the whole two-particle system must be carried 
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already by any individual particle! Namely, the wave function \l/;,/p(r) refers here to 
the two-particle system as a whole and thus, according to the RST philosophy, is to 
be conceived as the direct (Whitney) sum of the one-particle constituent wave functions 
^,(f) (a = 1,2) 

^,/p(f) = (^/fVi(r) © (''/^(r-) . (III.24) 

According to this sum structure, the total angular momentum is also the sum of 
the individual angular momenta 

Jz = ^'^1 © ^'^J. . (III.25) 

Therefore the eigenvalue equations (1III.23P for the total angular momentum jf^ can be 
decomposed as follows 

Consequently both one-particle spinors ?/'i(r) and 4j2{r) must obey the same eigenvalue 
equation as the total wave function flIII.24p . i. e. 



(^/^Vi (r) = ^^/f^j, h ■ (^/^Vi (r) (III. 27a) 

(2)J, (^/PV2(r) = h ■ (*/?'V2(r) . (III.27b) 

Furthermore, the Dirac four-spinors ipai'f^ can also be conceived as the direct sum of 
Pauli two-spinors ^"'^Lp_^_(r) and i-^- 

Mr) = ^"V+(0 © ^"V-(^) • (111.28) 

Therefore the same eigenvalue equation must also hold for the individual Pauli spinors, 
especially after those identifications Olll.lOap - piI.lObp : 



j{+) Wp)ip^{r) = h (III. 29a) 

Wp)^_(f) = iWp)j^ (^/pV_(f) . (III.29b) 
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This means mathematically that the eigenvalue ^^^^''jz of the two-particle state '^b/p ( 1111.23^ 



becomes transferred to any individual Pauli component of this two-particle state! In 
physical terms, the bosonic or fermionic character of the two-particle state \E' becomes 
thus incorporated in any individual constituent of the two-particle system. 

But now it is clear that positronium as a whole carries bosonic properties (-w is 
integer- valued) ; and this must therefore hold also for any Pauli constituent '^"•^ip±{f) of 
both particles (a = 1,2), cf. flIII.29ap -f lIII.29bl) . On the other hand, it is well known 
that the Pauli spinors do form a half- integer representation of the rotation group SO (3). 
This means that one can select in any two-dimensional spinor space a certain spinor basis 
{C'^e™} "with the following eigenvalue properties: 

Per = 3{3 + l)n\T (Ill.SOa) 

Xc^^i"' = mher (iii.30b) 

C^l"^ = l{l + l)h\'l'^ (111.30c) 
|2 ^i, ™ = + i);^2 ^i, m _ (III.30d) 

Here the electron/positron spin is s = |; the orbital angular momentum is £ = 0, 1, 2, 3, . . . 



and thus the lowest possible value of j(= £ ± s) is j = 2 with £ = or £ = 1. Therefore 
we have two basis systems for j = \, namely \Q'^]Q' ^ \ and s Ci ; Ci^' ^ \. 

Fermionic States 

The existence of these two basis systems admits us to decompose now a fermionic 
state in the following way 



(4 



1 1 , s 11 



v+(ri = (''>;^+(r-) ■ cr + ^"^s+i^ ■ Cr~' (ni.3ia) 



('^V-(^) = {^''^-(r) ■ CP + ^''^5_(f) ■ Cf' ^} . (III.31b) 
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The action of the angular momentum operator on such a fermionic state is obviously 
= (4 ^">/^+(r1) ■ Cr + (4 ^"^'5+(f)) ■ Cr~' (ni.32a) 



4 + ^)(">^+(r) 



Co" 



2 



1_ 1 

■ 2 ' 2 



j(-) Wy,_(f) = (l, ('^>;^_(f)) ■ cr - ^ (4 (''^^^f)) • cf 



(III.32b) 



4 + ^]('^>;^-(r) 



1 1 

2 



1 _1 

2 ' 2 



The required results (IIII.29ap - flIII.29bp of the action of the angular momentum operator 
Jz on the Pauli spinors '^"•^ip±[f) (with a = 1,2 or a = h/p) are now deducible from the 
present equations f lIII.32ap - flIII.32bp by making the following arrangements: 

Lz (''^i {r)=izh- ('*>;^± (f) (III.33a) 
Lz (")5±(f) = (4 + l)h ■ (")5±(r-0 . (III.33b) 

Indeed, with these disposals the equations flIII.32ap - flIII.32bp adopt the required form 
of the eigenvalue equations f lIII.29ap - flIII.29bp with the eigenvalue of angular momentum 
being found as 



Jz 



+ 



1 



(III.34) 



Since the quantum number of orbital angular momentum is adopted as integer [i^ = 
0, ±1, ±2, ±3, . . .), we actually end up with half-integer quantum numbers ^"•'j^ flIII.34p 
for fermionic states! 

These fermionic states can obviously be realized by use of unique amplitude fields 
^"•^±(r) , ^'^^S±{f) and also unique spinor basis fields Co''^^, Ci^'^^ • Evidently, the latter 
fields work as the carriers of the spin, whereas the amplitude fields do contribute the 
orbital angular momentum. If this philosophy is tried also for the bosonic states we are 
forced to give up the uniqueness of the spinor basis! 
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Bosonic States 

Joining here the general conviction that bosonic states should have integer quantum 
numbers ^"^j, f lIII.29aj) - flIII.29bjl . i. e. ^"^j, = 0, ±1, ±2, ±3, . . ., it suggests itself to think 



that the amplitude fields should furthermore carry integer quantum numbers (f^) of orbital 
angular momentum; i. e. such equations as flIIL33ap - flIII.33bl) should persist also for the 
bosonic states. Thus the necessary modification must refer to the basis spinor fields (-^l 
f lIII.30a|) - (]III.30dp . More concretely, we think of four basis spinor fields ^g^'', ^j^'', 



which for para-positronium obey the following eigenvalue equations 

Ji+) ej-) = Ji+) = Ji-' = Ji-' d"^ = ; (III.35) 
and similarly for ortho-positronium one wishes to work with a spinor basis 
rio^\ Vo'K vi^^ of the following kind: 

(III.36a) 
(III.36b) 
(III.36c) 
(III.36d) 

This says that for the para-case (1III.35I) the spins Sa (= |) of both constituent particles 
(a = 1, 2) do combine to zero spin quantum number sp of the para-type (sp = Si — S2 = 0, 
Si = S2 = |); whereas for the ortho-case (]III.36ap - (lIII.36dp the individual spins combine 















Ji- 




u {+) 









to unity {so = si + S2 = 1). Thus both basis systems f lIII.351) and flIII.36ap - flIII.36dp 



carry integer spin and therefore may be used for the corresponding decomposition of the 
Pauli spinors ^°'^Lp±{f) due to a bound two-particle system. 

For para-positronium {s-p = 0) one has now in place of the fermionic situation f]III.31ap - 
f ]III.31bp the following decomposition 



(^V+(r) = (^>7^+(rl ■ + ^"^'5+(f) • e^-^ (III.37a) 

-^|(^>;^_(f)■e^^^ + (^')5_(f)■e^^} . (iii.37b) 
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Here, the action of the angular momentum operator looks now as follows 



(-) 





(L,(^>;^+(r-))■e^+) + (L,(^)5+(f))■e 



(III.38a) 
(IIL38b) 



Consequently in order to satisfy again the para-form (p) of the eigenvalue equations 
(llll.29al) - (IIII.29bp one puts here 



(IIL39a) 
(IIL39b) 



so that the quantum number ^^^jz 0III.29ap - piI.29bp is solely due to orbital angular mo- 
mentum: 



(p). 



Jz 



(III.40) 



(4 = 0,±1,±2,±3,...) 



For ortho-positronium {sq = 1), the situation is somewhat different. First, the de- 
composition of the ortho-spinors ^^V±(^) looks quite similar to the para-case flIII.37ap - 
flIII.37bD : 



V.+ (r-) = (^>;^+(r-)■r^W + ('')5^.(f)■r^(-) 



ibyj^_ 



(f) ■ r^l^-' + 



(III.41a) 
(III.41b) 



But the action of the angular momentum operator Ji^"* on these ortho-states looks now 
as follows 



j(+) (b) 



■ 4"-' + 



(III.42a) 



J(-)(V(r) = -^|^f4 + /l)(^>7^_ 



■ 4^' + 



(L, - (^)5_(f) ■ r^S^^I . (III.42b) 
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For satisfying here again the eigenvalue requirement ( 1III.29a|) - (lIII.29b|) in its ortho-form 
(6) one puts in a self-evident way 

L, (^>;^± =i,h- ^^>R± (III.43a) 
L, ^^'^S± = {i, + 2)h- (^)5± . (IIL43b) 

This arrangement yields then the following eigenvalue equations for the ortho-spinors 

J± W^i(f) = (i^ + l)h- (^V±(r) . (III.44) 

Thus one finds the quantum numbers of the ortho-system: 

^'y. = 4 + 1 (III.45) 
= 0,±1,±2,±3,...). 

For a realization of the required basis spinors one takes the fermionic ba- 
sis < Co ' ^ , Ci ' M as the point of departure and introduces a general spinor basis 

through 



Cr (ni.46a) 

J-) ^ ,Q^--2 (III.46b) 

^(+) _ g-*6<^ .^il (III.46c) 

c^i") =e*^-Cf~^ . (IIL46d) 
Here it is easy to see emerging the following relations concerning angular momentum 

jWa;W = -(6-0 h.J^^ (III.47a) 

jW4-)= n-J,^ (III.47b) 

jWc.W = -(?)-l);i.c.W (111.47c) 

Ji-) c.;-) = (^?. - h-J-\ (III.47d) 
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Obviously, this w-basis contains some free parameter S (i. e. the boson number) and if 
this is chosen as S = | we obtain the desired ^-basis (1III.35I) for para-positronium; and 
if S is chosen as 5 = — | we obtain the r^-basis (IIII.36ap - (]III.36d|) for ortho-positronium. 



For S = we get back the purely fermionic C-basis (lIII.30ap - (]III.30dp . 

With the choice of the w-basis flIII.46a|) - flIII.46d|) the loss of uniqueness becomes now 



evident: since the original ^-basis flIII.30a|) - flIII.30dl) is unique (i.e. the C"'] "^(t?, </>) consti- 
tute a "unique" spinor field on the 2-sphere S^), the other two basis systems ^ flIII.35p and 
T] ( 1III.36ap - (lIII.36dp are double-valued. More concretely, for both the ^- and the ?7-basis 



(^-w 5 = ±1) one finds by performing one revolution around the z-axis (0 < < 27r): 

eg) (0 + 2n) = e±- -ej J (0) = J (III.48a) 
r,it^ (0 + 2n) = e±- -r^g^ (0) = -r/g^ , (III.48b) 

and this says that we need two revolutions around the ^-axis (0 < < Att) in order 
to return to the original basis configurations. Since we adopt all the amplitude fields 

unique scalar fields, the double-valuedness of the para- and 
ortho-basis becomes transferred to the para-spinors ^^V±(r) flIII.37ap - flIII.37bl) and ortho- 
spinors '•^V±(r) f|III.41ap - flIII.41bp and from here ultimately to the Dirac spinors '^v^o 



^p(f) = (pV+(0 © ^^V-(^) (III.49a) 
^^(f) = (^V+(^ © (^V-(^) • (III.49b) 

Summarizing, the Dirac wave functions \l/p(r') and \l/c)(r*) for ortho- and para-positronium 
must in RST be double-valued in the following sense: 

^p(r, ^9, + 2n) = -^p(r, ^, 0) (111.50a) 
^o(r, ^9, + 27r) = -*o(r, ^, 0) . (III.50b) 
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Uniqueness of the Physical Densities 

Naturally, in a fluid- dynamic theory (such as the present RST) the proper observable 
objects are the physical densities, such as those of charge, current, energy, linear and 
angular momentum, etc. A plausible condition on these densities is surely given by the 
demand that these physical densities should be single-valued tensor fields. But, as we will 
readily demonstrate, this condition can be satisfied also by non-unique wave functions; 
and this fact allows us to actually deal with such non-unique wave functions, as given 
for example by the double- valued positronium states flIII.50ap - flIII.50bl) . Therefore one 



wishes to have some condition on the wave functions which on the one hand admits their 
non-uniqueness but on the other hand ensures the uniqueness of the associated physical 
densities! 

Now according to the present RST philosophy, the non-uniqueness of the (Dirac) 
wave functions is to be traced back to the spinor basis, whereas the amplitude fields are 
furthermore required to be unique. Therefore the non-uniqueness of the wave functions 
is measured by the boson number S, cf. flIII.46ap - flIII.46dl) : and thus the condition of 
uniqueness of the densities is to be retraced to some condition for the fixation of the 
boson number 6. Such a fixation may be attained now by considering specifically the 
Dirac density ^""^Icq (f) and the three-current kai^^ which read in terms of the Pauli spinors 
^°V±(^ as shown by equations flll.55ap - fjll.55bp . Decomposing these Pauli spinors with 
respect to the w-basis flIII.46ap - flIII.46dp lets then appear the (Dirac) charge densities 



flll.55ap in the following form: 

4:71 

Evidently, these charge densities are unique in any case and therefore do not yet provide 
an immediate handle for fixing the parameter S. 

This situation changes now when one considers also the Dirac currents ka flll.55bp . 
which by their very definitions are always real-valued objects: 
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(^\^ = 1- . ^^yR_ + ^^)*s+ ■ ('^^^^ - (">^_ • (">;^+ - ^""^h^ ■ (111.52a) 

^""^k, = |e2^(^-^)<^ ■ C(,) - e-2^(^-^)* ■ C(.)| (IIL52b) 

W/t^ = i-z. ) i-yji^ . (-}jz_ + {-}jz_ . i-yjz^ _ . _ . I (111.52c) 

But here a nearby restriction upon the parameter S suggests itself, namely through the 
plausible demand that the Dirac currents ka ( ]II.55bp . with their components being spec- 
ified by 0III.52ap - (IIII.52c| . must be unique (!) vector fields over three-space (albeit only 
apart from the origin r = and the z axis = 0, vr). Evidently this demand of uniqueness 
reads in terms of the spherical polar coordinates {r, i?, 0} 

fc„(r,^?,0 + 27r) = 4(r,^9,0), (III.53) 

and thus the values of S become restricted to the range 

S = ^(n+1) (III.54) 

{n = 0,±1,±2,±3,...) 

which then entails also (half- jinteger quantum numbers for the z component of the angular 
momentum flIII.47ap - flIII.47dp of the spinor basis: 



= ± - 1 j = ±^ . (III.55) 

Here it will suffice to admit for the spinor basis (of the two-particle systems) exclusively 
the values 5 = ±i; other values of S come into play for bound systems of more than two 
fermions. 

Notice here that this (half- jintegrity arises as a consequence of the demand of unique- 
ness with respect to certain physical densities (i.e. Dirac current), whereas the corre- 
sponding integral quantum numbers of conventional non-relativistic quantum mechanics 
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are mostly traced back in the textbooks to the uniqueness requirement for the wave func- 
tions themselves (not the densities). The lowest values of Sz f llll.SSp are Sz = ±i for 6 = 
and Sz = 0, ±1 for S = ±i. Thus for the first case (5 = 0) we have a fermionic basis and 
for the second case (5 = ±|) one deals with a bosonic basis. In this sense, a Dirac particle 
is said to occupy a fermionic quantum state if the boson number'^ 5 of its spinor basis 
is zero {h = 0), and a bosonic quantum state if the boson number S equals ±|. Observe 
that through this arrangement the fermionic or bosonic character of the quantum state 
of a Dirac particle is defined by reference to the corresponding spinor basis. 
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IV Ortho-Positronium (5 = — ^) 

After the general RST logic for the occurence of the ortho/para dichotomy is suffi- 
ciently displayed, we will now turn (for the remainder of the paper) to ortho-positronium. 
For this specific two-particle system there occurs a further ambiguity, i. e. the ^^ortho- 
dimorphism" , which is not present in the para-configuration. The effect of dimorphism 
consists in the circumstance that the angular momentum quantization admits the emer- 
gence of two different charge distributions ^^^kQ{r) which, however, are both due the same 
configuration of quantum numbers! Since these two electrostatic charge distributions are 
differing slightly, they carry a slightly different electrostatic interaction energy, and this 
causes a slightly different binding energy. In this way, ortho-positronium does appear in 
RST in form of electrostatic doublets. Observe here that, for our subsequent discussion 
of this ortho-dimorphism, we will be satisfied with the electrostatic approximation, where 
magnetic effects are neglected completely! 

1. Mass Eigenvalue Equations for the Amplitude Fields 

The point of departure for our quantization of angular momentum is the mass eigen- 
value equation (1111.12^ in Pauli form. Here, we will adopt for the moment the spher- 
ically symmetric approximation where the interaction potential ^''^Ao{f) appears to be 
spherically symmetric (i.e. ^^^Aq^t) =^ [^Uo(r), r = ||r||). For this situation, the ampli- 
tude fields (^>/^±, fllll.41al) - (IIII.41bD can be assumed to be the product of a purely 
angular-dependent factor and a purely radial factor which then entails the splitting of the 
original eigenvalue system into a purely angular and a purely radial problem. As usual, 
the solution of the radial problem yields then the energy spectrum (i. e. quantization of 
energy) . 

To begin with, one inserts the decomposition of the Pauli spinors ^''V±(^) (ini.41ap - 



(]III.41bp into the Pauli equations 0111.120 and thereby obtains the following system of 
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eigenvalue equations for the amplitude fields *^^^±(r, -(9, 0), ^''^S±{r,i9, 



+ e 



-2i(j> 



dr ~'~ r dcp 2r 



(IV. la) 



M + M, 

h 



2i(j) 



dr r d(j) 
1 d^^^R. cot 7? 



r d-d 



(IV. lb) 



9r r 
1 d 



90 
cot 



2r 



(90 



(IV.lc) 



M -M^ 

h 



+ e 



2j0 



dr 

1 9 (*>;^_ 



r 90 



2r 



cot '(9 



9^9 



(IV. Id) 



h 



2 " 90 

This may appear as a relatively complicated system but it can be simplified by imposing 
some plausible requirements on the desired solutions. Naturally, these requirements will 
refer to angular momentum, which usually serves to classify the solutions of the energy 
eigenvalue problems for the few-particle systems. Thus it is surely reasonable to demand 
that the eigensolutions of the present ortho-system (lIV.lap - QIV.ldp be classifyable by 



means of the associated eigenvalues "^^^j^ (IIII.45P of Jz{=f Jz © Jz )■ Recall here our 
hypothesis that the total spin of a bound system agrees with the individual spins of both 
Pauli two-spinors which build up the common basis for each of the Dirac four-spinors (in 
the present case sq = Si + S2 = I, see the discussion below 0111.530 ). 
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In this sense, we try now to obtain solutions for the present ortho-system ( 1IV.la|) 
(llV.ldp which are of the following form: 



(IV.2a) 



sin^9 



(IV. 2b) 



Indeed, by this ansatz the eigenvalue equations (]III.43ap - (]III.45p for angular momen- 
tum are actually satisfied; and furthermore the original mass eigenvalue system ( llV.lap - 
(llV.ldp becomes transcribed to the new amplitude fields ^^^R±{r,§),^^'>S±{r^§) and thus 
reappears now as follows: 



dr 



sm' 



+ 



sm 



t9cost9- (*)^4 



or r 



r sin 



r sin d 



cot d ■ , 



dr 
sin 

r 



in^-i? d ^''^S-{r,'d) sin -(9 cos z^^^. ^ j-^-)^ 



9r 



■('')^_(r,79) + (%(r,^9)-(% 



1 

r sin^ -(9 



r sin^-i? 



cot^?- 



r, {}) 
r, 7?) 



r, 7?) 
r, 



r, d) 
r, 7?) 



r, -(9) 
r, 



M + M, 

h 



M + M, 

h 



h 



M -M, 



(IV.3a) 



c- ('')^_(r, t9) 



(IV.3b) 



c ■ 



(IV.Sc) 



(IV.3d) 

-c-^^^S+{r,^) . 
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This is still a system too complicated in order that there could be hope to find exact 
solutions, even if one could consider the ortho-potential ^''^Ao(r,'d) being prescribed from 
the outside. However, this potential does actually couple back to the ortho-wave function 
ipbii^, via the Poisson equation (]III.2ip . i. e. in terms of the new amplitude fields ^''^R±{r, -d) 
and ('')5±(r,^9) f llv:2il) -f llv:2b|) 



A(%(r,^9)=-as 



1 2 



r sin 



(IV.4) 



+ 



sini? 



1 2 



Product Ansatz 

In view of these complications it should be obvious that it is necessary to resort to 
some approximative procedure in order to extract the physically relevant solutions of that 
intricately coupled system flIV.3ap - flIV.4p . For this purpose, the spherically symmetric ap- 



proximation suggests itself where the angular dependence of the gauge potential ^''Mo(r, •§) 
is neglected, i. e. we replace '^^^y4o(r, •§) by its spherically symmetric approximation [^]y4o(r). 
But observe here that this approximation assumption does not entail the spherical sym- 
metry of the new wave amplitudes (^)^±(r,^9) and ^^'^S±{r,d) fHv:2il) - flIV.2bjl ! However, 



the angular dependence of the wave amplitudes becomes now manageable, namely by 
applying the following product ansatz: 

('')^±(r,t9) = f7«(t9)-"i?±(r) (IV.5a) 
^'^~S±{r,^)=gs{^)-''S^{r). (IV.5b) 



Indeed, substituting this ansatz into the eigenvalue system 0IV.3ap - pV.3dp enables us 
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to separate this system into two first-order subsystems, namely the angular system 

- + 1) cot ^ ■ gR{^) = £o sin' ^ ■ gs{^) (IV.6a) 

+ + 1) cot ^ ■ = ^ ■ gnm (IV.6b) 

av ^ ' sm w 

and the radial system 

^^jMI _ . _ [^Uo(r) ■ "i?_(r) + ^ ■ "^+(r) (IV.Ta) 

c- 



dr r r 

M + M, 



h 



dr r r 

M + M, 



h 



^^^^ + ■ + • "i?+(r) - ^ ■ "5_(r) (IV.7c) 

dr r r 



M -M. 



'5_(r) + [%(r) ■ "5'+(r) + — ■ "i?_(r) (IV. 7d) 



dr r r 



M-M, 

h 



c-"5+(r) 



The meaning of the newly introduced constants io and io will readily be clarified, but 
first we will perform a further simplification of the present ortho-system 0IV.7a|) - (]IV.7dp . 



Indeed, it may appear somewhat strange that all four equations ( ]IV.7ap - 0IV.7b|) for the 



four variables "5"-!- are necessary in order to determine the one single mass eigenvalue 
M^,. Such a system surely appears to be overdetermined; and therefore the overdeter- 
mination must be eliminated, namely by identifying the amplitude fields in such a way 
that there remain two independent amplitude fields {Q±{r), say) which satisfy a coupled 
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first-order system of only two equations. Tlie wanted identification of tlie amplitude fields 
is the following: 

"i?+(r) = "5+(r) = n+(r) (IV.8a) 
"i?_(r) = "5_(r) = fi_(r) . (IV.8b) 



By this identification requirement, the four equations (]IV.7a|) - (lIV.7d|) become concen- 
trated to only two equations, namely 

iMl + .iX(r)- I'l4„(r) = c . f!_ (r) (IV.9a) 

dr 2r a 



dn (r) 2— (£0 — £0) ^ _ _ jw[ 

+ ^ ^ ■ n^{r) + [%(r) ■ n^{r) = ^^-^c ^]+(r) . (IV.9b) 

dr 2r n 

Moreover, the identification requirement flIV.8ap - flIV.8bp yields also the constraint 

^'^Jz = l (io + io) . (IV.IO) 



2 

Of course, it is self- suggesting that those newly introduced constants io and io 
flIV.6ap - flIV.6bp must have to do something with the quantum number {jo, say) of the to- 



tal angular momentum of the ortho-states. In order to preliminarily reveal that supposed 
interrelationship we pass over to the non-relativistic approximation of the latter eigen- 
value system fjIV.9ap - flIV.9bp for the amplitude fields Q±{r). This non-relativistic limit 
may be attained by approximately solving the first-equation f lIV.9aP for Q_{r), yielding 

and substituting this into the second equation f ]IV.9bp . The result is a decoupled second- 
order equation for fi+(r) (being simply denoted by fi(r)) which then looks as follows: 

2 



d^n{r) idn{r)\ K 



+ ^ I = I ■ m (IV.12) 



2M \ dr^ r dr I 2Mr'^ \ 2 

- he [%(r) ■ Q{r) = (M, - M) c^ ■ Q{r) . 
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Obviously this looks like a Sclirodinger-like eigenvalue equation for the determination of 
the (non-relativistic) energy eigenvalue E^, 



= (M, - M) 



(IV.13) 



where the effect of angular momentum is expressed by the term containing io and io- 
In view of this strong analogy it seems reasonable to define the quantum number {Jq) of 
the total (i. e. spin plus orbital) angular momentum for the solution of the Schrodinger 
equation pV.12p by 

(IV.14) 



Jo 



Thus we have now two equations, namely fllV.lOp and flIV.14p . which relate the quan- 
tum numbers of angular momentum ^^^j^ and jo to the constants £o and io] and this then 
admits to express the latter in terms of the first ones (although in an ambiguous way) as 



'Jz =F Jo 



(IV. 15a) 
(IV. 15b) 



But obviously, the product of both constants is unambiguous: 



eo-io = ^h!-jh 



(IV. 16) 



Anticipating here the fact that the magnitude jo of angular momentum is never smaller 
than its z-component ^^^jz (i. e. jo > |^'''!72|)) we infer that both constants io and io must 
always be of opposite sign {io ■ < 0). 

Thus, one formally obtains from f lIV.12|) a true Schrodinger equation 



2M 



r dr 



2Mr2 



(IV. 17) 



E, ■ n{r) . 



60 



However, one should keep in mind that the potential f^^Ao(r) is not fixed from the outside 
(as in the conventional Schrodinger theory) but must be considered the solution of the 
non-relativistic Poisson equation (lIV.4p . i.e. 

A[%(r) = -a, (IV.18) 
r 

One guesses from this fact that the interaction potential [''ly4o(r) will depend on the ortho- 
state Cl{r) itself so that the eigenvalue problem appears highly intricate. Different states 
fi(r) will generate different potentials ^''^Ao{r) according to f lIV.18|) so that the energy 



eigenvalues E^, due to different states will be found to be different, too elimination of 
the conventional degeneracy (II.ISP ): and thus the RST energy spectrum of positronium 
must be expected to be much more intricate than predicted by the conventional theory 
f lI.18|) : already in the non-relativistic domain and in the electrostatic approximation. 

2. Angular Momentum Quantization in RST 

It is rather evident that the off-separated angular system (lIV.6ap - (]IV.6b|) represents 
an extra problem which may be tackled independently of the proper energy eigenvalue 
problem (lIV.lTp -f llV.lSp . On the other hand, the eigenvalues -E* of the latter problem 
will clearly be influenced by the solutions of the angular problem (IIV.6ap - (IIV.6bp . namely 



via the quantum number Jq of angular momentum f llV.Mp and in general also via the 
interaction potential ^^^AQ(r) since this potential "feels" the angular behaviour of the 
amplitude flelds via the Poisson equation (lIV.4p . 

Anisotropy of the Gauge Potential 

This anisotropic influence can be expressed even more distinctly by observing our 
product ansatz flIV.5ap - flIV.5bp together with the identiflcations f lIV.8ap -f lIV.8bp which 



then recasts that former Poisson equation f lIV.4p to the more concise form 



A(%(r, = -Ana, (^)fco(r, = -Ana, ■ ^%(r) , (IV. 19) 
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i. e. the charge density ^^^o (^) piI.Sip becomes now also factorized {^^^ko(r, -i?) = ^^^ko('}!}) 
^''^ko{r)) with the angular factor ^^^kQ{'d) being specified by 



An [ sm v 



(IV.20) 



and the radial factor ^''^ko^r) by 



r 

This factorization effect will subsequently allow us to develop a suitable approxima- 
tion procedure in order to manage the anisotropy of the interaction potential ^*''Ao(r, -(9). 
Namely, the formal solution of the Poisson equation flIV.19|) reads 

(''Uo(r, ^)=a. [ d'f' ^^l^^^^ (IV.22) 



where the charge normalization condition 



d^f^^'^koir,^) = 1 (IV.23) 



allows now also a factorized form: 



J dn = / ^ + sin^ ■ = 1 (IV.24a) 

j drr^ ^%(r) = j drr |(]^(r) + (^)| = i . (IV.24b) 

This, however, helps us now to get approximate results for the gauge potential ^''^Ao(r, -(9) 
from its integral representation (lIV.22p . For instance, the roughest approximation consists 
in neglecting completely the angular dependence of the radial integral in (IIV.22P (by 
replacing it through its angular average) so that the angular normalization condition 
flIV.24aP can be immediately applied which then yields 

oo _ _ 

An J r \\r — r \\ 
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whose non-relativistic approximation reads 



"U„(.r) = fl/^^ (,V.26) 

Air r r — r' 



and thus is the solution of the "spherically symmetric" Poisson equation (1IV.18|) . Conse- 



quently, if we wish to go beyond this spherically symmetric approximation for the inter- 
action potential (^Uo(r,^9) (Hv:22D we are forced to first work out the solutions Qni'd) and 
gs{'d) of the original angular eigenvalue problem (IIV.6ap -f lIV.6b|) . 



Solving the Eigenvalue Problem of Angular Momentum 

The close interrelationship between angular momentum of both particles and anisotropy 
of their interaction potential ^^^AQ{r^ -d) flIV.22p necessitates now to explicitly solve the cou- 
pled angular system flIV.6ap - flIV.6bp . For this purpose, one first decouples that system 
by differentiating once more which yields 

,2 ■ 1 \ 2 ^ 

■ gR{d) = (IV.27a) 



gs{d) = (IV.27b) 



where ^'^'^jz is defined by equation flIII.45p and jo by f]IV.14p . 



Once we have arrived at two decoupled equations one can try to separately solve each 
of both by some transformation of variables which renders the problem more manageable. 
Thus, putting 

a; = sini9 (IV.28a) 
gR{^) Gr{x) (IV.28b) 
gsi^) Gs{x) , (IV.28C) 
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the decoupled second-order equations (]IV.27a|) - (]IV.27b| ) adopt the following form 



(1-x^) 



d"^ Gr{x 



2\ dGui^x 



X — 



X 



dx 



(IV.29a) 



• Gr(x) 



x^ 



_ d^Gsjx) ^ 2_^\ d_Gs{x,^ 
^ ' dx^ \x I dx 



It is true, the coupled first-order system flIV.6ap - flIV.6bl) does not allow that both angular 
functions Qnifi) and gs{^) can simultaneously be written in the form claimed by equations 
flIV.28ap - flIV.28cP : but the fact, that each one of these two equations flIV.6aP and flIV.6bp 



can emerge in the form flIV.29ap - flIV.29bp . means that always one of both equations 



flIV.29aP or flIV.29bp may be realized. Thus, our strategy of elaborating the desired 
solutions works as follows: first, select one of both equations flIV.29ap - f1IV.29bp and 
determine its solution as a function of x f lIV.28ap . Then determine the associated angular 



function gn^'d), or gs{^), resp., according to the relations flIV.28bl) - flIV.28cP : and finally 
determine the residual angular function gs{'&), or gni'^), resp., from the corresponding 
other angular equation flIV.6ap . or flIV.6bp . resp. 

The very form of the second-order equations flIV.29al) - flIV.29bl) suggests to try a power 
series 

G'r(x) = J]pn-x" (IV.SOa) 

n 

= -x" . (IV.30b) 

n 

If this is substituted in the second-order equations flIV.29ap - flIV.29bp . one gets by the 
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standard methods the following recurrence formulae 



Pn+2 



(n + 2)(n-l)-(Wj; + l)(Wj,-2) 



■ Pn 



n(n + 2)-(jo + l)(jo-l) 



(IV.Sla) 



(IV.31b) 



(n + 2)(n + 3)-Wj,(Wj, + l) 

In order that these solutions be unique, we next have to fix the lowest order coefficients 
of both series. For the first series flIV.30ap we obtain for the lowest power rimin 




and for the second series fllV.SObI) 



^'^Jz > 
^'^Jz < 

("^J. > 



(IV.32) 



(IV.33) 



On the other hand, for — )■ oo the coefficients pn+2 and (Tn+2 would tend to pn and 
an, resp.; and therefore both sums OIV.30ap - flIV.30bl) cannot adopt definite values if they 
do not stop at some finite nmax- However, such a maximal value nmax of the power n can 
immediately be read off from the recurrence formulae (]IV.31al) - (]IV.31bp 



rir, 



jo + 1, (HESiD 

Jo - 1 , flIV.31bp 



(IV.34) 



These results enable us to closer specify now our original ansatze (lIV.30ap - (]IV.30bl) : 
namely 



Gr{x) 



(IV.35) 



E Pn-X'- 



< 
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and 



io-i 



Gs{x) = { 



„=-((fc)j3+i) 



^'^Jz > 



< . 



(IV.36) 



From these results one concludes that the range of quantum numbers jo and ^^''jz is the 
following: 



-(jo - 1) < ^'^3z < Jo 

-Jo < ^'^J. < {jo - 1) 



for Gr{x) fll\r35]) 
for Gs{x) (ll\r36|) 



with jo being an integer {jo = 0,1,2,3,...). 

As a brief concrete demonstration of the general method, we consider the situation 
with jo = 2 ^ —2 < ^''^jz < +2. Starting here with ^''^jz = 0, we first conclude from the 
arrangements (llV.35p - (llV.36p that the first case flIV.351) does apply where the function 
Gr{x) must then look as follows 

3 

Gr{x) = J2p^.x^ = p,x + P3X^ , (IV.37) 



n=l 



or by means of the recurrence formula (lIV.31ap 

Gr{x) = Pi(x — 2x^) 



(IV.38) 



with the constraint pi to be fixed (below) by the normalization condition. Next, since the 
counterpart Gs{x) of Gr{x) cannot exist for jo = 2, ^^^jz = we have to determine the 
second angular function gs{'&) directly from the first-order equation (lIV.6ap with gnij}) 
being deduced from Gr{x) flIV.38p as 



gj^{^d) = pisin^9 (1 -2sm^^) . 



This yields 



9sm 



4p 



- cos '(9 



(IV.39) 



(IV.40) 



f^o 
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Now that both angular functions are fixed up to the normahzation constant pi, we can 
pin the latter down by the normalization condition (lIV.24aP which adopts the following 
form for the present situation 

^lp;/«{[i-2si„-.]= + (i-)^ 

n V 



2 

,2 



sin ^cos "(9 



This ascribes to pi the value 



(1 + 4 



Here, one wishes to express the constant io through Jq and ^''^j^, see equation (jIV.lSbp 
which, however, offers two possibilities, namely either 



^^o = ^'bz+Jo +2 (IV.43a) 
^o = ^'^J.-Jo -2 (IV.43b) 



or 



^'^Jz-io -2 (IV.44a) 
'^'bz + Jo +2. (IV.44b) 



But these two possibilities differ merely in sign and therefore do not generate essentially 
different solutions of the eigenvalue problem under consideration (in contrast to the other 
cases with 7^ 0, see the table below). Indeed, both situations flIV.43ap - (IIV.43bp and 



(IIV.44ap - (IIV.44bp yield the same normalization constant pi f lIV.42p . i.e. 



Pi = ±W^, (IV.45) 



where the ambiguity in sign is immaterial. Thus the solutions gR{'&) f lIV.39p and gs{'&) 
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(1IV.40|) adopt their final form as 



sin ^9 (1 - 2 sin^ ^) (IV.46a) 



gg{^) = -2 cos^9 (IV.46b) 
a.b). 



The peculiar point with this result is that the essentially ambiguous situation (llV.lSaP 



f llV.lSbp does not entail the existence of two essentially different solutions for all allowed 
values of '■^^j^ but only for some ^''^j^ 7^ (see the table below). Indeed, that essential 
ambiguity of some solutions of the angular eigenvalue problem (lIV.6ap - (]IV.6b|) is just the 



core of the ortho- dimorphism, i.e. the dimorphism of the ortho-positronium configura- 
tions. Naturally, the most urgent question arising in this context must refer to the energy 
difference of such a dimorphic pair, which will readily be worked out. 
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Ortho- Dimorphism for jo = 2 



The table (precedent page) displays the possible configurations of the pair (ioi ^^o cor- 
responding to the quantum numbers = —2,-1,0,1,2 due to jo = 2. The associ- 
ated angular functions gni^), gsi''^) as solutions of the first-order system (lIV.6ap -( ]IV.6bp 



emerge ambiguously for ^^^jz = ±1 but are unique for ^''^j^ = 0, ±2. The case with ^^^jz = — 1 
is the simplest one and has already been treated in ref. [31] (in that paper the present 
quantum number Jq has there been denoted by io). 

The table of angular functions gji^-d), gs{^) (for jo = 2) obviously displays some 
regularities (or symmetries, resp.) as far as the permutation of both functions gni^'d) and 
gs{'&) is concerned. Surely, this must have to do something with the isotropy of three- 
space, especially with respect to the replacement '■^-'j^ —^^^jz- Indeed, if this replacement 
is associated with the following replacements (irrespective of an irrelevant change of sign) 

gj^(^) =^ gs(^) . sin^ ?9 (IV.47a) 
9s{^) ^ -^4^, (IV.47b) 



(IV.47c) 



then both eigenvalue equations f lIV.6aP -( ]IV.6bp become merely interchanged and the an- 



gular density ^^\i^{{f) is left invariant. 

In order to get a broader basis for our inductive reasoning, we extend the calculations 
to jo = 3, see the table below. Combining both tables for jo = 2 and jo = 3 we 
conclude that the angular functions gni'd) and gsii!^) are always unique for ^''^jz = and 
for ^^''jz = ±jo for arbitrary jo (= 1, 2, 3, 4 . . .). This is proven quite generally in App.A. 
For the other values of ^^^jz (i.e. < \^^^jz\ < jo) the angular functions gni'd), gs{^) are 
generally expected to be ambiguous. However, the extent of this ambiguity deserves an 
extra discussion where perhaps the energy difference of the ambiguous configurations may 
be taken as a quantitative measure of just that ambiguity. 
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Ortho- Dimorphism for jo = 3 

The possible configurations due to jc, = 3 are ^^^j, = 0, ±1, ±2, ±3 and display again the corresponding ambi- 
guities (for ^''^jz = ±1,±2) and uniqueness (for ^^^jz = and = ±jo)- Thus, one concludes that for general jo 
{— 1, 2, 3, 4, . . .) one always has 2jo + 1 — 3 = 2(je) — 1) dimorphic pairs and 3 solitary angular distributions. 



The somewhat hidden symmetries of the pair of functions Qr^^) and gs{i^) due to 
the interchange of the positive and negative 2;-direction = -v^ -i? = vr) becomes 
now more obvious by looking at the angular densities ^'^^kQ^'d) flIV.20p . The subsequent 
table presents a collection of all these angular objects due to jo = 1,2,3. From here it 
becomes immediately obvious that the ambiguous pair gR^i!)), gsii!^) due to a given ^''^jz 
(with < \^''^jz\ < jo) does generate the same couple of angular densities ^^^ko{'d) as does 
the pair gji{'d), gsi'd) due to —^^^jz- Consequently, the electric charge density ^''^ko^f) does 
not single out a preferential direction on the z-axis (contrary to the current density kb{r)). 
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Angular Density Wfco(^^) (IYTM) for the States Due to jo = 1,2,3 
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Quite intuitively, one would suppose that the shape of the charge density ^^^ko{i!}) is 
influenced by both the azimuthal quantum number and the total number jo- However, 
the above table says that the angular densities ^''^kQ{'d) are the same for 
and jo = 3, ^''^j^ = and also for jo = 1, ^''^jz = 0. Indeed, one can show that the 
angular charge density ^^^ko^i!)) for ^^^j^ = is the same for all values of jo, see equation 



(jA.Sp of App.A. Thus one becomes tempted to attribute the deformations of the angular 
charge distributions ^''^kQ{'d) to the action of "centrifugal forces", which are linked to 
the azimuthal quantum number ^''^jz rather than to the total number jo! In this sense, 
one expects that for ^''^jz = the charge distribution is concentrated in the vicinity of 
the rotational axis {'& = 0,7i); whereas for maximal ^'^^jz ^''^jz = ijo) the charge 
distribution ^^^k^ becomes rolled out to a disc-like shape. Indeed, the occurrence of this 
effect is demonstrated by the Fig. IV. A below. Moreover it should be considered a matter 
of course that such a deformation of the charge clouds will entail some change of the 
electrostatic interaction which then implies that the jz- degeneracy becomes eliminated 
(see below). This elimination of degeneracy, however, has nothing to do with the ortho- 
dimorphism which occurs for dimorphic partners belonging to the same pair of quantum 
numbers jo-, ^'^''jz] see Fig. IV. B below for a simple demonstration. 
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Fig. IV.A Rolling out of the Rotating Charge Clouds 



For small rotational quantum number ^^^jz (solid curve: ^''^jz = 0) the charge density 
Wfco(r, ^) = ^%o{^) ■ ^%o{r) flIV.20l - flIV.21j) is concentrated close to the rotation axis (i.e. 
z-axis; '(9 = 0). The angular pre-factor ^^^kQ{i!}) due to ^''^jz = is the same for all values 
of jo(= 0,1, 2, 3,...), i.e. 

sin^? ■ ^%(^9) = (sohd line) 

see App.A, equation flA.Sp . However, for extremal value of ^^^jz (i.e. ^''^jz = ijo) the 
charge cloud becomes rolled out to a disc-like shape through the centrifugal forces; dotted 
curve: ^''^jz = jo = 3 with 

sin^9- (^)fco(^9) = -sin^^? , 
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see the table on p. [7TJ Such an inertial deformation of the rotating charge clouds is to 
considered the origin of the elimination of the j^- degeneracy which itself occurs only 
the spherically symmetric approximation. 
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Fig. IV.B Dimorphism of Angular Density ^^^ko('d) 

for jo = 3, (^)7, = ±2 

Even for the same rotational quantum numbers jo and ^^^j^ tlie angular density ^^^kQ{'&) 
(lIV.20p occurs in two different forms, provided the quantum number ^''^j^ is restricted by 
the condition < \^''^jz\ < jo, see App.A. For the present situation {jo = 3, ^^^jz = ±2), 
the two forms of angular density are given by (see table on p. ITTj) 

sin^ , solid curve 

I (I sin^ ''^ ~ I sin^ i?) , dotted curve . 
The different shape of the electric charge clouds entails different electrostatic energy and 
therefore also different binding energies ( "ortho-dimorphism"). 
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V Anisotropy of the Gauge Potential *^'^)Ao(r 

Besides the matter fields, the second basic building block of RST is constituted by 
the gauge fields. Since the latter fields are coupled to the matter fields via the Poisson 
equations, the anisotropy of the matter fields is transferred to the gauge fields just by 
these Poisson equations. For the present electrostatic approximation, the gauge fields enter 
the theory only in form of their time component '•^''Ao(r) due to the corresponding four- 
potential ^''^A^. The field equations for the gauge fields are the (generally non-Abelian) 
Maxwell equations, which however for the present electrostatic situation reduce to the 
ordinary Poisson equation flIV.4|) . Here, the charge density ^^^ko{r) factorizes as shown 
by equation flIV.19|) : two examples of this dimorphism of the angular factor ^''^ko{'d) are 
given by the above tables at the end of Subsect.IV.2. Thus it should be evident that 
exact solutions of the Poisson equation (1IV.19|) can not easily be found. However, this 



should not be considered a bad situation since there are approximation procedures at hand 
which are able to mediate a rough picture of how the exact solution (^^Ao(r*) will look like. 
The point of departure for this picture is the ^^spherically symmetric approximation", 
which has already been briefly mentioned in connection with the equations flIV.181) and 
flIV:25|) -f Hv:26|) . 



1. Spherical Symmetry as a First Approximation 

The spherically symmetric approximation of the gauge potential has been based upon 
the assumption that this potential is S0(3) symmetric, i.e. we put 

(^Uo(f) [%o(r) (V.l) 
(r = ||rl|) . 



The corresponding Poisson equation fllV.lSp for the approximative potential ^^^Ao{r) is 
obtained through extremalizing the electro-static part ^^^[n] of the full energy functional 
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E[n] which of course must consist of both the gauge field part and the (Dirac) 

matter part see ref. [2] 

Ep] = ^""tp] + ^""^[u] . (V.2) 

Here, the matter part '■^-'Ep] will be readily specified and discussed, but in the present 
context we are rather interested in the electric contribution to the gauge part (^terni 

because the extremalization of the latter part already yields the desired electrostatic field 
equations. 

Indeed, the wanted electrostatic part ^^-1E[q] is composed of two contributions 

('=l[o]=4'^ + Ag)-Ng). (V.3) 

The Langrangean multiplier adopts the value A^^ = —2 and the first contribution (-E^"*) 
is the electrostatic gauge field energy 

4^) = J d'f[E,{r)'E,{r)) (V.4) 

where the electrostatic field strength Eb{r) of the ortho-configurations is given as usual 

by 

E,{f) = -V (%(f) (V.5) 

(this simple form of Eh{r) does exclusively apply to non-identical particles, see ref. [3^). 
For our presently considered spherical symmetric approximation (IV. the general form 
Eb{r) fIV.Sp does reduce to 

E,ir) -V [^Uo(r) = -e. , (V.6) 

(e) 



and consequently the electrostatic energy E^' (IV.4p becomes simplified to 



.(e) _ p[e] _ ^ f ^^^2 fd ^'%iry ' 
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A similar conclusion does also hold for the ^^Poisson constraint" N^^ which generally is 
defined as the difference of the gauge field energy E^^^ and its ^^mass equivalent" M^'^^c^ [2]: 

AT^^) = - m(^)c2 , (V.8) 

where the mass equivalent itself is defined through [Zitat] 

M^^^c^ ^-hc [ £f (*Uo(f) ■ % (f) . (V.9) 



For the spherically symmetric approximation, the gauge potential is assumed to be SO (3) 
symmetric, cf. fIV.ip . and the charge density ^^^ko (r) is of the product form described by 
f lIV.19p -f lIV.2ip . Therefore the mass equivalent fIV.Qp factorizes as follows 



oo 

M^^^c^ -he j dn j dr [^Uo(r) ■ ^%(r) 



(V.IO) 







But here we can resort to the angular normalization condition flIV.24ap for ^^^ko^i!)) and 
also to the non-relativistic version of the radial charge density ^^^ko{r) flIV.2ip 



i%(r) = ^ (V.ll) 



so that we end up with the final form of the ( "isotropic" ) mass equivalent 

oo 

M(")c2 ^ Ml'^lc^ = -he j drr [%(r) • ^?{r) . (V.12) 







With all these arrangements, the "spherically symmetric" form of the electrostatic 
gauge field functional is of the shape 

(^l[o]^[*p] = 41 + Ag)-Ng (V.13) 

with being given by f lV.Tp and 

= E^'l - Ml'^lc^ . (V.14) 
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Now that the specific SO (3) symmetric form of all individual contributions is clarified, 
one can carry through the variational process for ''^EfQ] (IV.lSp and thereby obtains the 
well-known Poisson equation (IIV.ISP for ^''^Aq^t). The standard formal solution is here 
given by equation (1IV.26P which we will readily exploit for the study of the corresponding 
anisotropy corrections. 

2. Anisotropy Corrections 

If one disregards for the moment the spherically symmetric approximation (IV.ip . the 
general Poisson equation is found as the extremal equation due to the general gauge field 
functional '■^•fepj. Furthermore, the use of the product ansatz flIV.5ap - flIV.5bp recasted 
the general Poisson equation flIV.4p to the more specific form fllV.lQp with the formal 
solution 

<-U„(r..) = a./.v'"*°;,?-!;':°<'-''. (V.15) 



Clearly, this is the anisotropic generalization (up to a factor of An) of the former result 
f lIV.26p for the spherically symmetric approximation. 

It is true, the latter approximation disregards completely the anisotropic features of 
the exact solution of the eigenvalue problem, but nevertheless it plays a certain role when 
considering now the anisotropic potential ^^^AQ{r,'d) f lV.15p . Namely, since it appears 
rather difficult to exactly calculate the integral of that formal solution f lV.15p we are 
forced to think about some approximative evaluation. For this purpose, one may resort 
to an expansion of the anisotropic potential ^''^Ao^r^'d) around the spherically symmetric 
potential [^Uo(r) flTOej) . i. e. an expansion with respect to the magnitude of anisotropy. 
In this sense, one may look upon the spherically symmetric approximation t^]y4o(r) as the 
leading term of such an expansion, namely such that it appears as the angular average of 
the more general anisotropic solution ^^^Aq^t,!!)) f lV.15p . i.e. explicitly 

[%(r) = I ^^oir,^) (V.16) 
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In order to build up now such an expansion with respect to the magnitude of anisotropy, 
one starts with the following expansion of the denominator in the integral representation 
dVTSD of W^o(r,'^) [2] 



(^^2.^^/2^1/2 (1^2.^^/2^3/2 \^ J 2 (^^2 ^/2 j 5/2 

5 r^r'^ i.A3 35 r^r'^ 

^2 (^2 + ^/2)7/2 [r-r ) + y "^^;7^;^;;iji7^ (^^-^ j +••• 
^ ^' 1 1 -I I 1 1 -'I I 1^ 

However, this expansion is not directly introduced in the integral of equation flV.lSp but 
we first define the properly anisotropic part ^^^A^^{r,'d) of ^^Mo(r, (IV.lSp through 

^''U""(r, ^) = ^*Uo(r, ^) - [%(r) , (V.18) 

so that the spherically symmetric approximation [''ly4o(r) actually appears as the expected 
leading term of the proposed expansion: 

{''Uo(r, ^) = [%(r) + ^''U""(r, ^) . (V.19) 

As a consistency requirement, the former averaging postulate (IV.16P entails now that the 
average value of the anisotropic part ■^^^^'^(r, ■(9) must vanish 

f ^ {b}A-^(r,^) = . (V.20) 
J 47r 

But once the anisotropic part of the gauge potential is precisely defined now, its 
power series expansion can be easily carried through. Here, the point of departure is the 
corresponding integral representation of ^^^A^'^{r,i!}), which is obtained by substituting 
both integral representations pV.26p and f lV.lSp into f lV.lSp yielding 

^''U-(r,^9) = as f ('^%(^') - ^] ■ ^%(rO , (V.21) 
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and if the expansion (IV.lTp is inserted herein one finds emerging a series expansion of the 
following product form: 

WA""(r,7?) = • W^'(r) + • W^"(r) + ^^'"(i?) ■ ^^'"(r) (V.22) 

Here, the angular factors appear first in integral form, and the lowest orders may easily 
be calculated exactly: 

= j dQ! - i-^ = (V.23a) 

= j dQ! (^%o{:d') - ^) = (V.23b) 

{dQ! # dd'd(P' sin^9') . 

The importance of these lowest-order results lies in the circumstance that they are valid 
for any angular distribution ^^^Q^-d)] namely the first result f lV.23aP is a consequence of 
the separate angular normalization condition ( 1IV.24aP and the second result (]V.23bp is 



an implication of the skew symmetry of the integrand with respect to the refiection of the 
two-sphere at the equatorial two-plane {d = ^). Moreover by the latter argument, all the 
angular factors of odd order are trivial, i. e. 



j dQ! (^^\{^') -^{f- ^')'"^' = 



(V.24) 



(n = 0,l,2,3,...) 



Consequently, the first non-trivial angular factor is '"^''^'"(■(9), followed by 

(''Ui'V) (^%(^^') --^{f- ^')' ' quadrupole (V.25a) 

{n^V(^) j - ^) . octupole 133] . (V.25b) 

However, it should be evident that the precise value of these integrals will depend on 
the actual density ^^^A;o(^); for instance, any member of a dichotomic pair of ortho-states 



(due to the same quantum number jo) may have its own ang ular factor 

etc. 

being different from that one due to its jo-partner. But a common property of all angular 
factors ■"^''^'"('(9) etc. is that their average on the two-sphere is zero 



0, 



(V.26) 



An J Att 

which in turn implies that the anisotropic correction ■^^^^"(r, t?) flV.2ip has also zero 
average over the two-sphere, cf. flV.20p . 

For the sake of brevity, we are satisfied for the moment with a treatment of only the 
lowest-order case flV.25aj) for jo = 2,3,4 and = ±1. But even for this restriction to 
only a few values of jo we can observe the ortho- dimorphism, see the subsequent table 
which transcribes the densities ^''^ko{'&) to the corresponding angular factors ^^^'"(-(9) 
according to the prescription flV.25al) . 



jo 


Jz 


Wfcnf^^) fllV.20) 


WA'"f^^) flV.25ap 


2 


±1 


sin i9 


_A [cos2^-i] 


1^ • (9 - 8 sin^ ^) 


3 
16 


COS^ d-Q 


3 


±1 


^sin?9 (l - fsin^ ^) 


3 
16 


cos^ ^ — I 




sint? (1 - f sin^ ?9 + | sin^t?) 


9 

32 


cos^ ^ — I 




4 


±1 


^sin^ (1 -4sin2^9 + ^sin"^9) 


9 
32 


cos^ — I 




i| sin t9 (1 - 4 sin^ ^ + f sin" t9 - | sin^ ^) 


51 
160 


cos^ — I 





This table demonstrates some striking features of the quadrupole correction term 
First, the quadrupole corrections ( ]V.25aP for both ortho-densities 

^''^ko^'d) due to jo = 2 differ merely in sign. This circumstance has its consequences as far 
as the energy (being concentrated in the anisotropic field configuration) is concerned, see 
below for the discussion of the corresponding degeneracy. But in order to see here some- 
what more clearly the origin of this strange effect, which partly (i.e. from the energetic 
point of view) suppresses the ortho-dimorphism just elaborated, one reconsiders the orig- 
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inal definition (]V.25a|) of ^^^A^^^i!}) and observes also tlie separate angular normalization 
(lIV.24a[) of the density ^''^ko{'d), together with the specific angular structure of the scalar 
product of unit vectors r, r ': 

^r-r'j = sin^ ^9 sin^ -(9' cos^(0 — 0') + cos^ "i? cos^ "i?' (V.27) 
+ 2 sin {} sin {}' cos -t? cos {}' cos(0 — (j)') . 

Indeed, these simple mathematical elements are sufficient in order to deduce, with a little 
bit algebra, the following general form of ^^^A^^^d): 

(1 - ,3 ■ K'rW ] cos'iy - ^ 

4 



cos^ d — 
3 



(V.28) 



where the integral K^^ly is defined by 

K^^l^ = j dn'cos^^' = 27ry" d^' sin ^' cos^ ^' ■ . (V.29) 

This result (1V.28P explains why both ortho-configurations due to jo = 2 and ^''^jz = ±1 
share the same angular dependence, apart from sign. The value of this integral K^^ly 
f lV.29p on both densities ^^^ko{'d) due to jo = 2 (see the table of densities ^^^ko^i!)) on p. [82] 
yields now 

{■J : first density , ^^^kr,('d) = sin'i? 
^ ' (V.30) 
^ ; second density , ^^^kQ{i!}) = sm-& ■ (9 — 8sin^'i9) . 

Inserting now these two possible values of the integral K^^l} general relation flV.28p 

yields just the two quadrupole corrections '^^^^"'(^9) for jo = 2 with their different signs, 
as displayed by the table on p. |82j As we shall readily see, this difference in sign of 
both versions is not sufficient in order to equip both corresponding solutions 

of the eigenvalue problem with different energies EJp . Thus the result will be that the 
lowest approximation order (with respect to the anisotropy corrections) cannot induce 
an electrostatic energy difference for the dimorphic pair of states due to no = 3, jo = 
2, ^''^jz = =tl- (See also some further comments on this in App.D). 

83 



However, as is clearly demonstrated by the lower half of that table on p. [821 this 
result does hold only for the non-relativistic state due to jo = 2 of ortho-positronium. 
For the excited ortho-states {jo = 3 and jo = 4), the angular densities ^^^kQ{'d) and 
therefore also the quadrupole corrections ^''^A^^\'d) of both dimorphic partners are seen 
to be distinctly different from each other (App.D) which then will entail also different 
electrostatic interaction energies. Of course, this energetic degeneracy of the dimorphic 
ortho-states {jo = 2, ^^^jz = ±1) does refer exclusively to the first anisotropy correction 
as expressed by ^^^A^^^i!}) (]V.25a|) ! The next higher approximation order ^^^A^{i)) ( ]V.25b|) 



will surely eliminate this first-order degeneracy. Indeed, this expectation can easily be 
verified by inserting both ortho-densities ^'^^kQ{'d) due to jo = 2 (see table on p. |82|) into 
the general prescription f lV.25bp for the quadrupole correction ^''^A^{'d) in order to yield 



the following results: 



^%(^) = ^ ^ WAV) = ^{^-3cos2^-icos'^^| (V.Sla) 



{%(^?) = ^.(sin2t9 + 9cos2^9) 

^ {b}A'^{^) = J-^\-7 + 30cos^^-15cos^^} (V.31b) 
256 

Uo = 2, ("^j. = ±1) . 

Obviously, this second approximation order ^''^A^{'d) displays no common angular depen- 
dence for the dimorphic pair as is the case with the first approximation order ^''^y4"'('(9) 
(]V.28p - (]V.30l) . but nevertheless the average of all multipole corrections ^''^A'^{'d) 
over the two-sphere is zero, as required by equation (lV.26p . However, in contrast to their 
trivial mean values, the contribution of these anisotropy corrections to the binding energy 
is readily seen to be non-trivial. 
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VI Energy of Matter and Gauge Fields 



It should be a matter of course that different interaction potentials ^^^Aq^t,!!}) for 
both positronium constituents will entail different binding energies E^-*^. This is surely a 
reasonable assumption and may be expected to hold also for the approximate solutions of 
the eigenvalue problem to be discussed subsequently. The corresponding approximation 
method is the same as has been applied in the precedent paper [2], where the RST energy 
functional Ep] is extremalized on a suitably selected set of trial configurations. 

The trial amplitudes Q{r) are postulated to be of adequate form and include some 
variational parameters (here /3, z/): 

Q{r) = Q,r''e-^' (Vl.la) 

Obviously, after separating off the angular dependence (see flIV.5ap - flIV.5bp ). the selected 
class of trial amplitudes Cl{r) is here of the spherically symmetric form which itself is 
surely an approximation to the exact solutions; but for the purposes of the present paper 
we will be satisfied with such a spherically symmetric form of the wave amplitudes Q{r). 

In contrast to this, the interaction potential ^''^Aq^t) is admitted to be anisotropic; 
however, it is assumed to be the sum of a spherically symmetric part ^^^Ao(r) plus an 
anisotropic perturbation ^''^A^'^{r,'d), see equation f lV.lOp . Furthermore, these two con- 
stituents of the interaction potential ^^Mo(r, t?) are assumed now to be kinematically in- 
dependent, where this independence of the anisotropic part ■^^^^^"(r, t?) refers exclusively 
to the radial auxiliary potentials ^''^A^^\r),^^^A^(r) etc., cf. the series expansion flV.22p : 
whereas their angular counterparts ^^^A^^\'d),^^^A^('d) etc. are fixed through the integral 
relations flV.25ap . f lV.25bp etc. Thus our selected trial configurations consist of a set 
of spherically symmetric fields |fi(r); ^^^Aq^t), ^^^A^^\r), ^^^A^(r), . . .| where only the trial 
amplitudes Q{r) are of a pre-specified form (cf. fIVI.lap - flVI.lbl) ) with undetermined varia- 
tional parameters /3 and u. In contrast to this, the spherically symmetric potential ^^^AQ{r) 
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and the radial auxiliary potentials ^^^A^^\r), ^''^A'^(r) etc. are first left unspecified and must 
afterwards be determined from their extremal equations, i. e. from extremalization of the 
energy functional E[q] to be understood as the anisotropic generalization of the spherically 
symmetric version (IV. 2p . The extremalization of E[n] with respect to the trial amplitude 
Q{r) occurs then through extremalization with respect to the variational parameters /3 
and u which themselves are entering also the potentials ^^^AQ{r)]^^^A^^\r),^^^A^{r) etc., 
namely via the solutions of the extremal equations for those potentials whose source is 
determined by just the square of the amplitude Q{r). 

Thus, summarizing the approximation procedure, we will end up with a certain energy 
function E''^'^J'(/3, u) of the two variational parameters /3 and u; and the stationary points 
of this function, being determined by the requirements 

9E^'v}(/3,z/) 



df3 

aEi'v>(/3,z/) 



(VI.2a) 
, (VI.2b) 



yield the equilibrium values Z?*,//,, so that the wanted energies appear as the values 
E^'v>(/3,,i/,) of that energy function E^'^J'(/3, u) on these special values (3^:, i/*. 

According to the sum structure of the ortho-funtional E[q] (IV. 2^ . its value En(/3, z^) on 



our selected trial configurations will appear as the sum of the matter energy z/) 
and the anisotropic gauge field energy E^\/3, v) plus the constraint terms: 

Ef,(/3, v) = 2Ekin(/3, 1^) + 1^) + 2A, ■ Nf,(/3, v) + Ag^- Ng^(/3, v) (VI.3) 

= Ef^>(/3, v) + 2A, ■ N^(/3, v) + Ag^- N^^^(/3, v) . 

Whether or not the Poisson constraint term (IV. Sp in its anisotropic form i.e. 

= E^^ - M^'^c^ , (VI.4) 

will contribute in a non-trivial way must be discussed separately (see below). In contrast to 
this, the other constraint term refers to the normalization condition for the non-relativistic 
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amplitude Q{r) 

I drr n''{r) - 1 (VI.5) 



and can always be put to zero by simply admitting only those trial amplitudes Q{r) which 
are normalized to unity (-^ Nf^ = 0). However, from formal reasons it is necessary to first 
include the constraint terms in the energy functional E[q] so that both the matter equation 
(IIV.17|1 and the gauge field equation (IIV.18|) do appear as the extremal equations due to 



that ortho-functional E[q] flV.2|) whose matter part appears now as 

= 2Ekin + 2A, ■ (VI.6) 

and similarly its gauge field part as 

^%n] = E^^ + Ag) ■ Ng^ . (VI.7) 

1. Energy of the Matter Fields 

First, we may turn to the kinetic energy which for the two-particle system is twice the 
one-particle energy E]^i^ 



oo ( / ~ \ 2 

^2 

2/ 



2m/ +^ 



^kin = 7777 / ^^M -7^ + ^ ■ ^ (0 ) ■ (VI.8) 







Here, the centrifugal term (~ j^,) is induced by the former product ansatz flIV.5al) - flIV.5b[) 



and evidently does account for the total angular momentum of the considered ortho-state. 
Substituting herein the selected form fIVI.lap - flVI.lbp of the trial amplitude fl{r) yields 
for the one-particle kinetic energy Ei,i^ fIVI.Sp 

^ki„(/3, i^) = ^{2aBl3Y ■ (VI.9a) 

see equation (VI. 154b) of ref. [3l]. This is formally the same kinetic energy function 
as was used for para-positronium with merely the quantum number 1^ (= jp) for para- 
positronium being replaced now by jq for ortho-positronium. 
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2. Energy of the Gauge Fields 

Next, the gauge field part '■'^■fepj is to be specified in terms of the three gauge field 
constituents ^^^AQ{r),^''^A^^\r),^''^A^ (r) which we treat as independent constituents of the 
gauge field subsystem. Since, on account of the specific form of the Poisson constraint 
(]VI.4|) . the desired gauge field part ^'^'^[n] (1VI.7|) consists essentially of the gauge field 
energy E^^ and its mass equivalent M^^J'c^ it is sufficient to explicitly display these two 
objects in terms of the newly introduced gauge fields [^ly4o(r), ■>^^^y4"'(r), ■>^''^y4^(r). 

Turning here first to the electrostatic gauge field energy eI^\ the spherically symmetric 
relation (IV. 7p is to be generalized for the presently considered anisotropic situation to 

or r ov 

so that the gauge field energy (lV.4p reappears as 

4.) ^ 4.. = /■ I ('^"'^f''>) ' + 1 ('2}^^] J , (vi.li; 



But here we have to observe the splitting of the anisotropic potential ^''^Aq (r, ■(9) into the 
sum of the isotropic part [^Uo(r) and the anisotropic perturbation i^^A^^^r, cf. equation 
(IV.19P ; and this must then induce a similar splitting of the gauge field energy E^^ (IVLTT|) 
into two parts, i. e. we put 

Ei^^ = E§+Eif. (VI.12) 

Here, the "isotropic" part E^ has already been specified in terms of t^]y4o(r) by equation 
(IV. 7p and the "anisotropic" energy Ein is found to be of the following form 



AttOs J I V dr 7 '"^ V '^'^ 
(See the precedent paper [1] for the reason why there is no mixed term involving both 
parts [%(r) and ^''U^^(r, 

For further specifying down this anisotropic part E^'n^ (I VI. 13 P one observes the de- 



composition (IV.22P of the anisotropic correction ^''^A^^{r,i!}) into the sum of products of 



an angular and radial auxiliary potential; here we retain only the first two terms of the 
series expansion (1V.22I) : 



(VI.14) 



This arrangement lets then reappear the "anisotropic" gauge field energy m 
the following final form 



M!> = -- 



as 



dr r 
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em 
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III/ 



+ (^Vi 



III 



V ■ 



(VI. 15) 



eiv 



dr 



dm'^ir) 
dr 



IV- 



The angular coefficients cn, /n = III, IV, V) do arise here by carrying out the corre- 
sponding angular integrations in the original energy integral ( IVI.13P which yields 

2 . ,^ .111/ nx N 2 
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em 



dn 

Air 



W^iii(^) 



III 



dQ fdm^^\^) 
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eiv 
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d^ 



(b) 



ev 



dQ 



(b) f 

Jv 



dn (dm^i^d) 

Ait 



(VI.16a) 
) 

(VI.16b) 
(VI. 16c) 



displaySince the angular functions ■^^^^'"(■(9), ^^^A^(i3) are known and depend on both quan- 
tum numbers jo and ^''^j^ but depend also on the dimorphic type, the angular coefficients 
eN, /n (N = III, IV, V) are functionals of the considered ambiguous quantum state (jo, ^'^^jz) 
but are otherwise to be treated as fixed numerical constants. We will readily see that 
through this arrangement the principle of minimal energy becomes reduced to a purely 
radial variational problem. The numerical values of the constants are displayed by the 
subsequent table for the quadrupole approximation (~ '^^^^'"('(9)). 
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jo 


Jz 




(ivr 


16a|) 

(b) f„, 
Jill 


flVI.22a) 
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±1 


3 

16 


cos^ ^ — I 
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3 
160 


1 
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3 
16 


cos^ — I 




1 
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3 
160 


1 

64 


3 


±1 


3 

16 


cos^ — I 




1 

320 


3 
160 


1 

64 


9 

32 


cos^ ^ — 1 




9 

1280 


27 
640 


9 
256 


4 


±1 


9 
32 


cos^ ^ - 1 




9 

1280 


27 
640 


9 
256 


51 
160 


cos^ — I 




1 /17\2 
5 \80J 


1 /SIX 2 
30 Uo/ 


\80J 



The striking feature of these numerical values refers to the fact that the ratios — and 



^ are all the same (see App.A) 



"' 6 (VI.17a) 



= 5. (VI. 17b) 



Wmiii 



The meaning of this will readily become obvious when considering below the quadrupole 
equation for the radial potential correction ^^M"'(r). 

3. Mass Equivalent 

Similar arguments as for the anisotropy energy flVI.lSp do also apply to the mass 
equivalent M^^-'c^ (1V.9|) . Indeed, inserting therein the series expansion (1V.22P (up to the 



V-th order) in combination with the factorized density ^''^ko (r) (IIV.20|) -( IIV.21|) and the 
separate angular normalization (lIV.24aP lets emerge the (non-relativistic) mass equivalent 
in a form analogous to the "anisotropic" gauge field energy E^^ (1VI.12P : 

M(")c2 ^ M^">c2 = Ml^lc^ + M±c^ . (VI.18) 

This says that the general mass equivalent splits up into the sum of the spherically sym- 
metric part MI'^^c^ f|V.12p and the anisotropic correction M'^-anC^, where the latter part 
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appears as the sum of the 1 1 l-rd and V-th order, i. e. 

M^f^c^ = Millie' + M^v^c' . (VI.19) 
This yields finally the following decomposition 

M^^^c^ = M^^lc^ + M^l^ijc^ + M^^;}c^ (VI.20) 

where the spherically symmetric part Mt'^'c^ has already been specified by equation f lV.12p 
and the anisotropy corrections are given by 



I^^>c2 = -he ^''^rrnu ■ j drr ^^U"'(r) ■ n\r) (VI.21a) 



oo 

I^^>c2 = -he (*)mv ■ [ drr ^''U^(r) ■ n^{r) , (VI.21b) 





with the angular constants ^''^m\\\ and '^^-'my being defined through 

Wmiii = J dn ■ ^%(^) (VI.22a) 

('')mv = [ dn • ^%('^) . (VI.22b) 



The angular coefficient ^^^mm (lVI.22a|) is obviously due to the formerly considered quadrupole 
approximation and some values have already been displayed in the table on p. [901 

The last point concerning the gauge field subsystem refers to the "anisotropic" Poisson 
constraint flVLil) . Since both constituents E^^ and M^'^J'c^ have been split up into 
sums of the isotropic and anisotropic parts, cf. flVI.12p and (IVI.lSjl . the constraint itself 
must undergo a similar splitting: 

= Ng + , (VI.23) 

where the isotropic part Nq is of course built up by the isotropic parts of the gauge field 
energy and its mass equivalent, cf. (IV. 14^ . and analogously for the anisotropic part N^^n 



^ = Eif - M^t^' . (VI.24) 
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Summarizing, we ultimately arrive at a splitting of the whole energy functional E^^] 
(1V.2|) not only in two but rather three constituents which correspond in a self-evident way 
to the matter part and the isotropic and anisotropic gauge field parts: 



E 



2Ekin + 2A, ■ 



+ 



+ 



(VI.25) 



(^1 

(^1 



These three subsystems will now be adopted as three kinematically independent (but 
coupled) degrees of freedom so that we can deduce the corresponding field equations by 
extremalization of this functional E^j^] flVI.25p on the space of selected trial configurations 
n(r); [''Uo(r), WA"'(r), WAV(r)' 



4- Principle of Minimal Energy 

Once the energy functional for the RST configurations has been set up, cf. flVI.25p . one 
can now look for the corresponding extremal equations. For this purpose, one may derive 
benefit from the fact that the various terms of E[n] do not contain all four independent 
quantities '^'Aq, ■>^^M"'(r), ■>^''M'^(r)|. Therefore the variational process for any of these 
field quantities does comprehend only a certain subset of all those terms constituting Ej^j . 
For instance, the extremalization of E[q] with respect to the amplitude field Cl{r) does 
disregard the electrostatic gauge field energy E^^ fIVI.llI) and thus takes into account 
only the matter part '■^•5E[f7] flVI.6p and the mass equivalent M^^^c^ flVI.20p . Thus the 
extremalization {6q, say) with respect to the amplitude Cl{r) looks as follows 



[n] 



:[^ic2 + 



(VL26) 



and this yields for Q{r) the Schrodinger-like equation (provided we put Ag^= —2): 



2M 



cPVL(r) ^ 1 dVL{r) 



r dr 



+ 



2Mr 



(VI.27) 



(b) 



mm 



■ W^"'(r) + ('')mv ■ ^^U^(r)} ■ f2(r) 
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Obviously, this is the anisotropic generahzation of the former eigenvalue equation (1IV.17|) 



which thus is seen to refer exclusively to the spherically symmetric approximation (IV. ip . 



If the anisotropic corrections ^''^'"(r, ■(9) and ^''^A^(r,'d) are taken into account, their ra- 
dial parts ^^^A^^\r),^^^A^{r) do obviously also couple to the spherically symmetric trial 
amplitude (l{r), where ^''^m and '^^Wy flVI.22al) - flVI.22bp play the role of coupling con- 
stants. 

Clearly, we are not able to solve the amplitude equation (IVI.27P exactly, but instead 
we will resort to the previously used trial functions Q{r) (IVI.lap - llVI.lbj) . Thus, denoting 
the value of the matter functional ( ]VI.25|) on this normalized trial amplitude by 



2i?kin(/3, z/), one readily is led back to the former results (]VI.9ap -( 1VI.9b|) . 

Isotropic Energy Function E['^](/3,z/) 

Next, concerning the extremal equation for the spherically symmetric approximation 
[''l74o(r), one observes that this isotropic potential does enter the functional ]E[n] flVI.25p 
only via its second constituent f^Ejf^] fIV.lSp . Therefore the "isotropic" extremalization 
process {Sq, say) becomes active only on this "isotropic" constituent ^'^^[a\ 

5oE[n] =5ot*[n] =0 (VI.28) 

so that of course the "isotropic" Poisson equation (IIV.ISP is recovered here. But the 
point with our present approximation procedure is now that we do not try some (more 
or less reasonable) potential [^Uo(r), but rather we prefer to solve exactly that isotropic 
Poisson equation, albeit for the situation where the amplitude Q{r) is given by our trial 
ansatz flVI.lap - flVI.lbl) . For this purpose, one recasts the Poisson equation fllV.lSp in 
dimensionless form 

"^^^^ = -iWT^ = " r(2. + 2) , ' ^^'-^'^ 
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see equation (A.l) of the precedent paper P; here the solution can easily be found for 
integer 2z/ as 

This can be shown to be a special case of the general situation with arbitrary real values of 
the variational parameter v (> — |) where the desired solution then adopts the following 
form: 

~a^{^) = (i-e-y-y -2/2-+" ) . (VI.31) 



But once the isotropic potential f^^Ao(r) (or a,^{y), resp.) is known, one can readily 
calculate the associated "isotropic" energy (\V.7\i and its mass equivalent M'^^c^ (IV. 12^ . 
In the dimensionless form, both objects appear as [21 |3l] 

°° -22 

Eg = -e' (2/3) J dyy' (^^) = (2aB (3) ■ e,.,{u) (VI.32a) 





oo 

1 r p2 



= -e' (2/3) / dyy ~a,{y) 9.t{y) = {2a^ (5) ■ £pot(z/) (VI.32b) 







with the "isotropic" potential function Spoti^) being given by [H [M] 



— (^-—■f: 

\ n=0 



n r(4z/ + 2 + ra) 

2z/ + 1 P 24'^+2 ^ 2^ ' r(2z/ + 2) •r(2t/ + 2 + n) 



Observe here the fact that the electrostatic gauge field energy flVI.32aj) equals its 
mass equivalent M'^'lc^ (lVI.32bp so that the "isotropic" Poisson constraint (1V.14P is zero 



N^Si/3,u)^0, (VI.34) 

and thus this constraint term disappears from the energy function ]En(/3, u) due to (lVI.25p . 
The reason for this is that for the isotropic potential [''Uo(r) (or ajy(y), resp.) we use an 
exact solution of the corresponding Poisson equation, e. g. in the dimensionless form 
(IVI.30p and flVI.3ip : and for this situation it can easily be shown through integration by 
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parts and use of the Poisson equation (]VI.29[) that both integrals in (lVI.32al) - (]VI.32b| ) 
are actually identical: 

^pot(^) = ldyy' (^^) = / dyy ~aM ^liv) ■ (VI.35) 



At this point of the discussion it may be instructive to discontinue for a moment in 
order to get some survey of what has been attained up to now. Evidently, with both 
the matter part ^^^[q] and the "isotropic" gauge field part ''^E^qj being at hand now, 
one has at one's disposal the isotropic part of the energy functional (1VI.25p . Indeed, 
with the results (lVI.9a|) -f lVI.9b|) for the kinetic energy i?kin and with the "isotropic" gauge 
field energy (lVI.32aP one possesses knowledge of the corresponding "isotropic" energy 
function E['^1(/3, z/), i.e. 



El'^] i^) = 2Ekin(/3, ly) + 4'^ ly) (VI.36) 
{(2aB/3)^ ■ ekinl^^) - (2aB/3) ■ epot{j^)} ■ 



Ob 

This result may now be exploited in order to obtain the (non-relativistic) energy spec- 
trum of ortho-positronium in the isotropic approximation. To this end, one subjects 
this "isotropic" function E['^](/3, u) (IVI.36P to the first one (lVI.2ap of both extremalizing 
conditions in order to find the corresponding equilibrium value as 

2a.^ = . (VI.37) 

Substituting this back into the original energy function E['^1(/3, u) yields its reduced form 
as 

E(iv)(/3, u) ^ Et(z/) = --^ ■ Soiu) , (VI.38) 

4aB 

with the "isotropic" spectral function ^^(z/) being given by 

Saiiy) = • (VI-39) 
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Finally in the last step, one looks for the extremal points of this function 

dSo{i^) 



(VI.40) 



in order to find the equilibrium value (z/*) of the second variational parameter (z/). 

In this way, one gets the energy spectrum |e^|?^(z/*)| = jlEo'l of ortho-positronium 
in the spherically symmetric (i. e. "isotropic" ) approximation. Of course, this spectrum is 
numerically the same as for para-positronium, if calculated on the same level of approx- 
imation, see the table on p. 11631 Observe however that there is an important difference 
between both energy spectra. Namely, the presently considered ortho-levels are dimor- 
phic doublets, due to the occurence of different charge distributions ^^^ko{r, '&) for the same 
quantum numbers jz,jo (= — 1). This dimorphism of the ortho-states will induce a 
certain energy difference between the members of such a doublet which is readily to be 
discussed now in greater detail (the degree of degeneracy of the ortho-levels will be studied 
in an other paper). 

Anisotropic Energy Function u) 

Naturally, the "isotropic" energy function E['^1(/3, z/) (1VI.36P must undergo certain 
corrections if the "anisotropic" constituent '•^"•fep flVI.251) 

^^"^i^[n] = i5^i? + Ag)-Nt> (VI.41) 

of the total functional Ep] is taken into account. The anisotropy corrections ^^^'"(r) and 
{-^^J^y [r) of the interaction potential ^^^y4o(r, t?) are explicitly present in this anisotropic 
part ( ]VI.4ip so that the extremalization (denoted by and (5v, resp.) of the total 



functional E[q] with respect to those correction potentials ^^^'"(r) and ^^^A^{r) becomes 
active only for that anisotropic part, i.e. 

5111 E[n] = 5,11 (^'^>E[n]=0 (VI.42a) 
5vE[Q] = 5v^'"i:p] = 0. (VI.42b) 
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Carrying here through both extremahzation processes yields the following coupled Poisson- 
like equations for the anisotropy corrections ■"^''^'"(r), ■f^^M^(r): 

■ A. i^U'"(r) + (%,v ■ A. WA^(r) - (^)/,„ ■ - (^)/,v • (VI.43a) 

= ^miiia, 



r 



(^)ev • A, ^^U^(r) + Weiv ■ A, ^^U'"(r) - (^)/v ■ - ^^Viv ■ (VL43b) 

= — ^ ''my as • . 

r 

If it were possible to find exact solutions of the present coupled system (IVI.27P plus 
flVI.43ap - flVI.43bp . then the "anisotropic" Poisson constraint term flVI.24p would be zero 

^ . (VI.44) 

But clearly, it is very hard to get these (surely existing) exact solutions; but even if we 
manage to get exact solutions only of both coupled Poisson equations flVI.43ap - flVI.43bp 



alone (i.e. for arbitrary source fi(r), e.g. f lVI.lap -f lVI.lbp ). then the "anisotropic" Pois- 
son identity flVI.44p will nevertheless hold. This may easily be verified by multiplying 
through the coupled equations f lVI.43ap -f 1VI.43b"|) with {''^^'"(r) and ^^^A^(r), resp., and 



integrating over whole three-space under regard of E^^ flVI.lSp and its mass equivalent 
M^anC^ flVI.19D . For such a situation where both Poisson identities flVI.34p and flVI.44p 



are valid, the anisotropic part '■'^'^■fejQ] flVI.4ip gets rid of its constraint and thus becomes 
reduced to its proper physical part E^^ , i. e. 

^''^n] Ei:} . (VI.45) 

But the same conclusion does hold also for the isotropic part, i. e. the "isotropic" 
Poisson constraint term Nq f lV.14p vanishes whenever an exact solution of the "isotropic" 



Poisson equation fllV.lSp is used, see equation flVI.34p . Taking this altogether we get rid 
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of all the constraint terms in the original energy functional E[n] (1VI.25P which then also 
becomes reduced to its proper physical part E'^^'^J': 

tin] E^'^> = 2Ekin + 4'^ + ^in^ . (VI.46) 

The value of this reduced functional on the set of our trial configurations, parametrized 
by the variational parameters {/3, u}, yields again some energy function E,^^'^^{(3, u) whose 
stationary points according to (IVI.2ap - flVI.2bl) determine the energy spectrum of ortho- 
positronium in the next higher orders beyond the spherically symmetric approximation. 
Naturally, one expects that there will emerge some energy difference between the cor- 
responding ortho- and para-states and also between the dimorphic partners of ortho- 
positronium. With the precedent preparations, these questions will now be discussed in 
some numerical detail. 
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VII Energy Difference of Dimorphic Partners 



The Poisson-like equations flVI.43ap - flVI.43bl) provide us now with the possibihty to 
study the energy difference of both dimorphic partners. In order to keep our demon- 
stration of this effect as simple as possible it may be sufficient here to study only the 
case ^''^jz = ±1 for jo = 2, 3, 4. Indeed the results, obtained so far, say that the di- 
morphism cannot occur for ^''^j^ = and ^''^j^ = ±jo, ^jo] and this excludes a priori its 
occurrence for jo = 0, 1. 

It should be evident that such an energy difference is in the first line due to the 
corresponding anisotropy corrections ^^M"'(r) and ^^^A^ [r] which stand for the first and 
second order of our anisotropy approximation, cf. the perturbation expansion flV.22p . 
In contrast to this, the kinetic energy E'km is assumed to be still of the simple form 
flVI.9ap - flVI.9bl) . Thus one may raise the question whether perhaps the expected energy 
difference emerges already in the first approximation order (~ ^'^^A^^\r))l If yes, is then 
the magnitude of the second order-correction (~ ^^^A^ {r)) considerably smaller than the 
first-order correction (~ ^^^A^\r))l Obviously, these questions can now be settled by 
first solving one (or both) of those Poisson-like equations flVI.43ap - flVI.43bp for given 
source VL{r) (e.g. fIVI.lap - flVI.lbp ) and then minimalizing the resulting energy function 
E'>^''^J'(/3, z/) with respect to both variational parameters (3 and u. This will yield the 
binding energy of ortho-positronium for '^^'>jz = ±1 due to jo = 2,3,4, quite similarly 
to the discussion of the isotropic situation (cf. flVI.36p - flVI.40p ). The result will be 
that, for jo = 2, ^^^j^ = ±1, the energetic degeneracy of both dimorphic configurations 
survives the first anisotropy correction (~ ■>^^M"'(r)); but in the second approximation 
order (~ ^^^A^(r)) the degeneracy becomes eliminated, albeit only by 0, 04 [eV] (i.e. 0, 6 % 
of the binding energy), see the table on p. 11451 
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1. Quadrupole Corrections (~ ^''^A"'(r)) ; jo = 2,3,4; ^^^jz = ±1 

For a first estimate of the anisotropy corrections one may accept for the moment the 
hypothesis that the second-order corrections (~ ^^^A^{r)) will turn out as being much 
smaller than their first-order counterparts (~ ■>^^M"'(r)) so that in the lowest-order ap- 
proximation we can neglect ^^^A^{r) altogether and thus arrive at the following truncated 
form of equation ( IVI.43aP : 



■ A. i^U'"(r) - (^)/,„ . = -as (^)m,„ ■ ^ . (VII.l) 

Here, the constants ^''-'em, /iii) '•^''"^iii are defined by the equations (lVI.16aP and (lVI.22aP : 
and this clearly demonstrates their dependency on the angular density ^^^fco('i9). However, 
since the latter object is of dimorphic nature, see the two possibilities for ^^^^id) and 
the associated '^''^'"(■(9) in the table on p. [H21 one expects to get two different sets of 
constants {^''•'em, '•^-'/iii) ^''"'^iii}) which then would entail two different equations of the 
present form f lVII.ip . But such a dimorphic equation f lVII.ip would imply the occurence 
of two different solutions for the anisotropy correction ■>^^M"'(r); and this, in turn, would 
yield two different anisotropy energies 

Furthermore, since the anisotropy energy is an essential part of the total energy Ef^>, cf. 
flVI.46p . one would also end up with two different energy functions WS^'^^ [13 ^ u) ^ namely 
as the values of that functional E^'^J' flVI.46p on both trial configurations ^^^A^ir)^ 
^*^74"'(r)}. But these two energy functions would immediately include the existence of a 
certain energy difference between both dimorphic partners! 

Actually, in the considered quadrupole approximation there is no energy difference be- 
tween both dimorphic states due to jq = 2, ^^^j^ = ±1 ! The crucial point here is that, de- 
spite the dimorphic nature of ^''^ko(-t}) and ■'^^^y4"'('i9), all three constants ^''^em, /m, ^^^w^iii 
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turn out to be identical, no matter whether they are calculated by means of the first trial 
configuration (-w first line of the table on p. [82|) or the second one (-w second line of that 
table). For both cases one finds the same result: 



""e,„ 



l§i'''^"'Wr-m (V...3a) 

'"/,„= /^("T'^'V'- (vn.3b) 



y dfi ^^U"V) ■ = ^ • (VII.3c) 



Of course, this identity of '-^■'eiii and f\\\ for both dimorphic configurations due to j'o = 2 
comes not as a surprise, because the quadrupole correction '^^^'"('(9) differs only in sign 
for both dimorphic partners; but on the other hand, the equality of '-^■'min (IVII.ScP for 
both dimorphic partners is surely somewhat amazing since the angular density ^^\i^{{f) 
is distinctly different for both cases (recall the table on p. [5^ . However, the situation is 
different for jo = 3 since in this case all three constants *-^^eiii, '■^-'/iii; '■^''^iii are different for 
both dimorphic partners; and thus it is merely the ratio fl VI . 1 7ap - (I VI . 1 7b|) what remains 



identical. Therefore one expects that the dimorphic partners will have different binding 
energy exclusively from jt'o = 3 on (in the quadrupole approximation). 

Quadrupole Equation 



But in any case, with the first-order coefficients (lVII.3a|) -( IVII.3cP being fixed in such 



a way, the equation (IVII.ip becomes specified (for all cases due to jo = 2,3,4) to the 
same quadrupole equation" (recall the ratios (IVI.17aP - IVI.17bp 



r2 —^r- . (Vn-4) 



where an example for the trial amplitude Q{r) has already been specified by equations 
fIVI.lap -f lVI.lbp . This equation has been studied extensively in connection with para- 
positronium in the precedent paper (see ref. PJ). It is possible to find the exact solution 
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of this quadrupole equation (1VII.4|) for the case where the trial amphtude Q{r) is given 



by equations (IVI.lap - llVI.lbp . For integer values of the variational parameter 2z/ the 



appropriate solutions read in dimensionless units {y = 2(3r; = {2(3as) ^■■'^^^'"(r)) 

{ 2u+3 n-6 ,^ 
n=6 ■ m=0 ^ J 

where for u = 1 the second term (in the bracket) must be put to zero; and the solution 
^^^Ai\y) for u = 1 is given by 

^ I ^ ^' J 

Indeed, the case for u = 1 plays an important part for this set of solutions because 
it fixes the boundary conditions. Namely, observe that the homogeneous (dimensionless) 
form of the quadrupole equation flVII.4p 

(a.-^) ^'U!j'(i/) = (VII.7) 

admits two simple solutions 

W^!j'(y) y'ory-' (VII.8) 
from which we conclude the limit behaviour 

lim ^^Mjj'(?/) = const ■ (VII.9a) 
lin, i^U-(v) = P-^ + ^'f + . (VII,9b) 

The second limit {y — )■ oo) is immediately obvious from the combination of both equations 
fIVII.Sp and flVII.6P : and for verifying the first limit {y — ?■ 0) it may suffice to look at the 
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lowest-order cases = 1, |, 2, |: 



^'^X/2(y) = Y ■ ^^U^y) - Y^y^e-^y^ + lOy + 66) 



6(^'-6^'+14^'-48^'---) 
1 

60 





(VII. 10a) 


) 


(Vll.lOb) 




(VII. 10c) 




(VII. lOd) 



Once the quadrupole solutions are exactly known one may substitute them into the 
first-order form of the anisotropy energy E^^^ flVII.2p which then must equal its mass 
equivalent M^[^|jc^ flVI.21ap . i.e. the value of the Poisson constraint N^^n (I VI. 24 p on our 
selected trial configuration |fi(r), ■i^*M"'(r)| is again zero 



Nt>(/3,^)^Nf,^ =i?,^^-Mf,^V = 0. 



(VII.ll) 



Here, the anisotropy energy E^^^ flVII.2p adopts the following form 



/rfWA"'(r 



dr 



(VII.12) 



em — i2^aB) ■ e';,Uu) 
as 



=F I dyy 



and equals its mass equivalent Mfi^'^c^ (]VI.21a[) which may be rewritten as 

oo 



<'>m,„ — (2aB/!).f,«((/) 

ttB 



(VII. 13) 
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so that the claimed Poisson identity (IVII.lip can here exphcitly be vahdated by reference 
to the solution (]VII.5|) of the quadrupole equation ( 1VII.4P in the form 



^'^e,„-e;^LM= (V,M-/iSf(z/). (VII.14) 

Properly speaking, the present final form of the energy functional Ejj^], i.e. E,^^'^^ 
(lVI.46p . would require now the use of the anisotropy energy Ej^^^ (IVII.12P : but since the 
calculation of this object is somewhat tedious one may resort to its mass equivalent Mfil-'^c^ 



(IVII.13P which is known to be numerically identical to the desired eI^^ because of the 
vanishing of the Poisson constraint n[|'|-'^ f lVII.lip . A closer inspection of the mass equiv- 



alent f lVII.lSp reveals the following structure for the involved mass-equivalent function 

/if?(^) = j^^2u+2) I y"^^' -^'^^'iy) (VII.15) 

r(2z/ + 4)-r(2z/-l) , 1 r(4z/ + 3 + n) 



r(2z/ + 2)2 r(2z/ + 2)2 ^^"^ ' 2*'^+^+' 



where the coefficients Pnit^) are given by 

, , r(2z/ + 4) T(2u- 1) ,,,,, , 

1 (2z/ + + n) 1 (2z/ + n) 

2. Energy Function of the Anisotropic Configurations 

Since the total gauge field energy E^^ flVI.121) appears as the sum of the isotropic 
and the anisotropic parts, with the isotropic part E^ being given by equations flVI.32ap - 
flVI.33p and the anisotropic part Ein by equation flVII.12p . one ultimately obtains the 
gauge field energy in the following form: 

2 

4'^ = -— (2aB/3) {^pot(^/) + ^'^miii ■ /if,f (VII.17) 



as 

^-^{2asP)-e\i{v) 



as 
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Obviously, the dimorphism enters the total gauge field energy E^^ via the mass-equivalent 
parameter '^^-'mm which in general has been shown to adopt different values for the di- 
morphic partners, see the table on p. [901 The important point with this result is that, 
also after the inclusion of the anisotropy correction, the total gauge field energy E^^ still 
persists in form of a linear function of the first variational parameter /3. On the other 
hand, the kinetic energy E]^[^ has already been found to be a quadratic term of /3, cf. 
f lVI.9aj) -f lVL9b|l : and therefore the total energy function E^'^}(/3,z/) emerges as the sum 



E^'^>(/3, u) = — |(2as/3)' ■ e^^) - (2aB/3) ■ 4o1(^^)| • (VII.18) 

Or L J 



of a linear and a quadratic term of (3, i. e. 

,2 

This "anisotropic" energy function is of the same form as its "isotropic" precursor 
E['^1(/3,z/) flVI.36p and therefore the extremalization process runs in a quite similar way, 
cf. equations flVI.36p - flVI.40p . i.e. we ultimately have to look for the minimal values of 



the spectral function 6*0 (z/) being defined through 

s}^\^) ^ \ f ; , (VII.19) 

where {j} symbolizes the pair configuration of quantum numbers {jo, ^''''jz}- Here ^^t (z^) 
is to be taken over from equation flVII.lTp 

45(^) = ^pot(^) + (')m,M ■ ^^(i.) , (VII.20) 



and £kin('^) is specified by equation f ]VI.9bp . Denoting the maximalizing parameter for 



the spectral function S^\u) f lVII.19p by 



dS^^\u) 



dv 



the corresponding energy spectrum "K^I^ is then obtained again from the local maxima of 



(VII.21) 
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E 



o 



4a/ 



6,8029... [eV] 



(VII.22) 



The results for ^''^jz = ±1, jo = 1, 2, 3, 4 are collected in the following table: 



Jo =^ 





1 


2 


3 


4 


^% 





±1 


±1 


±1 


±1 


Wmiii =^ 


1 


1 


1 


1 

64 


9 
256 


16 


M 


M 


9 
256 


289 
6400 


[eV] 


-7.5378 . . . 


-1,575. . . 


-0.6811... 


-0.3808... 


-0.2501... 


-0.3900... 


-0.2533... 



As expected, the ortho-dimorphism does not occur for quantum-numbers jo = 
and jo = 1- The lowest-order possibility for its occurrence is jo = 2, ^'^^jz = ±1. But 
now the important point with this result is that the calculated energy E^'^^ (IV11.22P IS 
(in this first-order approximation to anisotropy) numerically the same for both dimorphic 
configurations with their rather different charge distributions ^''^ko{'d), see the upper half 
{jo = 2) of the table on p. |82j Clearly, such an amazing result calls for an explanation. 
But here the first observation is that both angular anisotropy corrections '^^^^"'('(9) differ 
only in sign, cf. the table on p. |H21 and therefore the coefficients ^^^m\\\ f lV11.3ap and '-^■'/iii 
f lVII.3bp must be numerically identical for the dimorphic pair with jo = 2. Furthermore, 
the coefficient ^''^m\\\ (IVII.3cp turns also out to be the same for both configurations. This, 
however, entails that the first-order equation (IVII.ip must be identical for both dimorphic 
configurations, see its common form equation (IVI1.4p . But for this situation it is clear that 
the corresponding solution ^''^"'('(9) of that "quadrupole equation" is also the same, since 
for both cases one takes the same trial amplitude Cl{r) fIVI.lap - flVl.lbp . Furthermore, the 
latter argument does also apply to the spherically symmetric approximation ^^^Aq^t) (IV.16P 
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so that now all three fields |f2(r), [^]74o(r), ^^M"'(r)| are identical for both dimorphic 
configurations {jo = 2), albeit only in the lowest order of the anisotropy approximation. 

But this is sufficient in order that all three energy contributions to E^'^> flVI.46jl 
are also identical; and thus both configurations {2,±1} carry the same energy Kq' 
(IVH.22p . In physical terms, this says that the degeneracy of the dimorphic partners 
(occuring in the spherically symmetric approximation) is not eliminated by the present 
first-order anisotropy corrections (see Fig. VII. A below)! See also App.D for entering 
into particulars. 

3. Comparison of Ortho- and Para-Levels 

The reference basis of the comparison of the various RST levels is the conventional 
level system Eq'' (]I.18|) . solid line. The various RST approximations are the spherically 
symmetric approximation e)^"^ (■ ■ ■ ) and the quadrupole approximations for para- 
positronium Ep -'^ (xxx) and ortho-positronium E^ -'^ (•••). 

All three RST approximations E^, eJ,^'\ eJ/'^ are based on the partial-extremalization 



process which relies on the maximalization of the spectral function Sq^^^, cf. flVII.22p 



kotW+ (''^)m,M-/if,f(z/) 
= . . .^ ■ (Vn.23) 



1 

For the spherically symmetric approximation E^^ [IJ, one puts jo.p = 0; and concern- 
ing ('''P^miii for the quadrupole approximation of para-positronium one reads off the mass- 
equivalent parameter ^P^m\\\ from the table on p J163j and finally for ortho-positronium one 
takes '^''•'miii from that table on p. 11631 The conclusions of these numerical results for the 
spherically symmetric and quadrupole approximations are now the following (Fig. VII. A): 
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Quadrupole Approximation 



-0.23909 
-.2501 ^ 

-0.2533 



n=5 



E(4,l} 
O 

E{4,1} 

p 



O 




0.3736 
-.3808 H 
-.3900 



7[4] 



o 

P 



EI3J} 
O 



-.42518 



n=4 



.(4) 



Fig. VILA Comparison of Ortho- and Para-Levels for Principal 

Quantum Numbers n = 4 and n = 5 
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i) RST and Conventional Predictions 

All three RST approximations ■ mj^Kp^ (xxx), E^^ (• • •) yield nu- 

merical predictions above the conventional E^^j fll.isp ( — )! This may be under- 
stood as a consequence of the fact that the RST energy functional E^q] flV1.25p 
has a lower bound which is in the immediate neighborhood of the conventional 

(n) ~ 

levels E^cy Therefore the value of this RST functional Ejj^] on a (non-optimal) 
trial amplitude f2(r), such as our variational ansatz (lVI.lal) - (]VI.lbp . will yield pre- 

(n) 

dictions which always surmount the corresponding conventional prediction E(^ 
(i.e. e!j!'\eIP ,E^^^ \ > Efq^). Thus one may expect that a more clever selection 
of the trial amplitude Q{r) will shift all three kinds of RST predictions further to- 

(n) (n) 

wards the conventional E^ . Of course, one would like to know how close to -E^j 
the RST predictions will turn out for the (hitherto unknown) exact solution of the 
RST eigenvalue problem. 

However, one qualitative feature of the RST spectrum should be mentioned in this 
context: namely, a closer inspection of the larger principal quantum numbers (n > 
10, say) yields the result that the RST predictions E^^ become smaller than their 
conventional counterparts E^'^ (i.e. E^^ < E^^ for n > 10) even for the non-optimal 
trial amplitude Q{r) (VI. 1). Therefore the true RST energy eigenvalues should also 

(n) 

be smaller than their conventional counterparts E^ for sufficiently large principal 
number n (or no,np resp.). 

ii) Elimination of the j^-Degeneracy 

The second observation concerns the elimination of the j^-degeneracy which has 
been discussed in great detail in the precedent paper [2]. This means that the 
C^jo.p + l)-fold degeneracy of the conventional and the spherically symmetric RST 
levels (due to any fixed jo,p) becomes eliminated by the anisotropy of the elec- 
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trostatic interaction potential. This anisotropy may be thought to come about by 
the action of the centrifugal forces on the rotating charge clouds ( ^'''^^jz 7^ 0). We 
assume that this degeneracy elimination for para-positronium [2] does occur in an 
analogous way also for ortho-positronium, but we do not study this phenomenon in 
the present paper which is rather concerned with the effect of dimorphism (see also 
App.D). 

) Dimorphism of Ortho-Positronium 

The effect of dimorphism acts physically in a similar way as do the centrifugal 
forces in the case of the j^- degeneracy; namely by anisotropic deformation of the 
electric charge distributions, see Fig. IV. A. However the difference between both 
deforming mechanisms is that the ortho-dimorphism is an intrinsic geometric feature 
of the eigenvalue problem (IIV.6a|) - (IIV.6b|) itself; i.e. an indirect consequence of the 



presence of spin. If there is no angular momentum {^^^jz = 0, App.A), then there 
is no dimorphism (see the tables on p. 11451 and p. I152p . as one may expect also 
intuitively. But a somewhat counter-intuitive effect occurs for maximal angular 
momentum ^''^jz = ^jo-, where there is also no dimorphism (see App.A). A hint 
at the possibility that the present quadrupole approximation may not display the 
true occurrence of the ortho-dimorphism comes from the quantum number jo = 2. 
Here it is possible to satisfy the condition < I'-^^j^l < jo for the occurrence of 
dimorphism, i.e. by choosing ^''^j^ = ±1; but nevertheless the dimorphism does not 
exist (see the table on p. I145p . at least within the framework of the quadrupole 
approximation. Therefore the next step must consist in inspecting the octupole 
approximation (Sect. VIII), especially with respect to the dimorphism. 

) Ortho- and Para-Levels 

Finally, it is also interesting to notice the specific way in which the para- 
levels Ep"'-'^ become splitted into the ortho-levels E^-'^ and E^-*^ by the dimorphism. 
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The first observation is here (see table on p. 11451 and App.D) that the ortho-levels 
do agree with the para-levels for those states which are not subjected to the dimor- 
phism: E^ -'^ = This is exemplified here for jo = and jo = 1 (see the table on 
p. 11631) . but presumably it holds for all values of jo- But even if the dimorphic split- 
ting is active (here: jo = 3,4; ^^^j^ = ±1), one of the two split levels (i.e. E^"*^) does 
furthermore agree with the corresponding para-level K^P (Fig. VII. A on p. 11081) . 
The other split level of the ortho-configuration (i.e. E^ -*^) is shifted upwards off from 
the residual ortho/para-level. This enhances the deviation from the corresponding 
conventional prediction Eq\ It remains to be checked whether this result survives 
the consideration of the next higher (i.e. octupole) approximation. 



Ill 



VIII Multipole Solutions 

Naturally, one cannot be quite sure that all the essential features of the ortho-spectrum 
are already accounted for by the quadrupole approximation. The uncertainty refers here 
mainly to three questions: 

i) is the elimination of the '^^^j^-degeneracy a true effect or an artefact of the quadrupole 
approximation, 

ii) if it is a true effect, is then the order of magnitude of the level splitting width 
correctly accounted for already by the quadrupole approximation, or do the higher- 
order multipole results exceed/reduce the magnitude of the quadrupole predictions, 

iii) is the "incidental" missing of the ortho-dimorphism ioi jo = 2; ^^^jz = ±1 an artefact 
of the quadrupole approximation and thus disappears again when passing over to 
the multipole approximations of higher order? 

Subsequently, we will try to answer these questions by a first preliminary inspection of 
the multipole solutions of higher order. 

1. Alternative Multipole Expansion 

The perturbation expansion flV.22p of the anisotropic part ^^^A^^{r,'d) of the origi- 
nal gauge potential ^^^AQ{r,'d) fIV.lQp may be further elaborated by observing the spe- 
cific form of the two lowest-order contributions (table on p. 1821) and ^^^A'^('d) 
flV.31aD - flV.31bjl . Indeed, the latter two pertubations flV.31aj) - flV.31bD can be recast to 
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the following form 



(V.31a): 



(V.Slb): 



sin '(9 



(VIII. la) 



105 
256 



cos 



35 
512 



cos d — - 



sin-i? 



(sin^ ?9 + 9 cos^ 



(VIII. lb) 



105 
128 



cos 'd 



105 
256 



cos 'd 



Thus, one is lead to suppose here that the original anisotropic expansion (]V.22|) yields 
an angular expansion in terms of the specific functions 'C^'^^^{'d) being defined by (/ = 



1,2,3,4, 



C^2'}(^9) =cos''^ 



1 



2/ + 1 



(VIII.2) 



i.e. the original anisotropy expansion flV.22p would then look as follows 

oo 

^''U""(r,79) = ^C^^^\d) ■ ^^U2/(r) . 
1=1 

Properties of the New Basis Functions 



(VIII.3) 



The new basis functions C^'^''^{'d) have some interesting properties which we will later 
on exploit for our purposes of attaining better approximations for the wanted ortho- 
spectrum. First, observe here that the angular average of all the basis functions is zero: 



An 



(VIII.4) 



so that the anisotropic part i''}A^'^{r, f lVIII.3p of the gauge potential actually does obey 
the requirement f lV.20p of vanishing angular average. Next, split up the Laplacean A^^^i,,/, 
into the three contributions due to the spherical polar cooedinates r, {}, 



^r,'&,6 — Ar -\ 5- H o Arf 

y.2 r^sm ^9 



(VIII.5) 
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with 



92 2d 

92 ^ d 

902 ' 



(VIII.6a) 
(VIIL6b) 
(VIII.6c) 



and then find that the longitudinal part A^ (]VIII.6bp of the Laplacean A^ ,9 acting on 
the basis functions C'^-^'j' (-(9) as follows 



A^ C^2i}(^) _ 2/ {21 - 1) • C^2«-2}(^) _ 2/ {21 + 1) C^2i}(^) 



(VIII.7) 



i. e. the set of basis functions is closed under the action of the longitudinal Laplacean A^. 
From this we conclude that the Poisson equation for the anisotropic part ■f^M'^°(r, ^9) 



AWA^°(r,^9) 



-47ra<, 



47r 



(VIII.8) 



will transcribe to a system of coupled Poisson-like equations for the radial potential cor- 
rections ^^^A2i{r) (IVIII.Sp . see below. In order to deduce that system by means of the 
principle of minimal energy it is necessary to define a scalar product {Ci^rn} for the vector 
space spanned by the basis functions. A nearby choice for this is: 



(VIII.9) 



2m + 2n + l {2m + l){2n + 1) ' 

For instance, for a treatment including terms up to the octupole approximation ^^^A^{'{}) 
f lVIII. la|) -( 1VIII. l"b|) one has to work with the following (symmetric) coefficient matrix 



{Q,m} 



A. _8_ 

45 105 

8 16 

105 225 , 




49 



105 
2 



105 245 
2 4 



(VIII. 10) 
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The quadrupole approximation (~ ''^''^A" '(■(?)) is a subcase hereof where only one single 
matrix element Ci^i = ^ is needed. For the higher-multipole approximations one corre- 
spondingly extends the dimension of the matrix {Q m} according to (IVIll.Qp . 

We will readily see the relevancy of these coefficient matrices {Ci^m} which obviously do 
represent a generalization of the former coefficients 

occuring in connection with the original perturbation expansion flV.22p . The analogous 



generalization of the former coefficients ^''^f\\\, ^^^fw, '■^■'/v, • • • flVI.16ap -f lVl.16cp would now 
obviously be a certain matrix {-Fi,m} which is to be introduced by means of the derivatives 
of the new basis functions C^'^^^{'d), i.e. we define: 

The two-dimensional case < 2, m < 2) is again needed for the octupole approximation 
and is found to look as follows: 

^ ( ! 1 ) • 

\35 63/ 

The case of the higher-multipole approximations should be self-evident now. 

Finally, one has to inquire after the higher-multipole generalization of the mass- 
equivalent parameters '•^Vnin, ^^Wy, . . . flVI.22ap -f lVl.22bp . For this purpose, one restarts 
from the general definition of the mass equivalent M'^^'^c^ flV.Op and inserts therein the 
new anisotropy expansion fIVlll.Sp . i.e. we then obtain the "anisotropic" mass equivalent 



I^anC^ in the following form 

OO 



H = -he V ^^^mi / drr Q\r) ^^U2/(r) = V M^^fc^ (VI11.13) 

Jo 



1=1 1=1 

where the new mass-equivalent parameters ^^^mi are given by 



= dn C^2'}(^) . (V111.14) 
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In contrast to the precedent matrices Ci^m and Fi^m the new parameters ^'^^mi depend now 
on the angular distribution ^^^ko^i!}) and are thus not merely numerical objects of general 
character. 

The new definition (1VIII.14P of the mass-equivalent parameters looks very similar to 
the old one (IVI.22ap - flVI.22bp : and indeed, both definitions are almost (but not really) 



identical. Indeed, there is a certain (albeit technical) difference between the original 
perturbation formalism (IV.22P and the new one (IVIII.3P if one stops at a finite per- 
turbation order; although both formalisms would yield the same result if the infinite 
number of perturbation contributions could be taken into account. In order to see this 
technical difference somewhat more clearly, consider (e.g.) the peculiar situation with 
jo = 2, ^''^jz = ±1 (see the table on p. |H2])- Here, the new prescription (1VIII.14P yields 
for both dimorphic densities ^^^/co(''^) the following results for the mass-equivalent param- 
eter ^''^mi up to the octupole approximation 



i%(.) ^ ^ ^ i '^'"^ ^ ' (VIII.15a) 

(9 / 8 . . \ / ^'^"^1 = 1 

\ / I = ^ ^ octupole 

{jo = 2, ("^J. = ±1) . 

Thus the new quadrupole parameters are found to differ in sign ({f-lmi = and are 

some three times larger than their counterpart ^^Vn\\\ of the original formalism 

Wmi = T- ■ ('W,„ = ^^.- = ^-, (VIII.16) 
1^3 ^ 3 64 ^12 ' ^ ' 

see the table on p. 11631 Thus, on the level of the quadrupole approximation there is actu- 
ally no essential difference of both formalisms because the objects are of the same order of 
perturbation and differ merely by some unsubstantial numerical factor. It remains to be 



I i^}m2 = — Jt , octupole 



= sin^A ^ I = " ■ (VIII.15b) 
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shown (see below) that such a difference of the intermediate objects of both approaches 
does not generate any difference of the final predictions of physical relevance. 



Comparison of the Expansion Coefficients 

However, the analogous comparison of the octupole contribution of both formalisms 
reveals now some (albeit purely formal) distinctions. For the sake of simplicity, we restrict 
ourselves again to the specific state in question (jo = 2, ^''^jz = ±1) and observe that the 
octupole potential term f lVIII.lajl -f lVIII.lbjl emerges as a decomposition with 

respect to the new basis C^'^^^'d) flV111.2p : 

2 

WA^(^) = a'^^^ C^2Z}(^) (VIII.17) 
1=1 

Such a finite series expansion of the original angular modes ^^^A^^\'d), ^''^A'^{'d), ^''^A^^^('&), 
... in terms of the new modes C^'^''^ (i!)) seems to occur for any perturbation order. The 
expansion coefficients functionals of the angular density ^^^kQ('&), 

for instance for the present octupole cases ( IVIII.lap - llVIII.lbp 



105 

_ ) 256 ' 

105 
128 ' 

35 

_ ) 512 ' 

105 
"256 ' 



first density (IVllLlah 

(Vlll.lSa) 

second density ( IVIII.lbp 



first density fIVIII.laP 

(VIII. 18b) 

second density (IVIILlbl) . 



Naturally, in view of the definitions of both mass-equivalent parameters *^^Wv flVI.22bp 



and ^^^mi fIVlIl.Mp the expansion fIVIII.lTp can now be transcribed also to these param- 
eters: 

2 

(^)mv = ' (VIII. 19) 
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i.e. concretely for the present demonstration (]VIII.15a|) - (]VIII.15b[) 

{jo = 2, ("b". = ±1) 

I -^-(-T^) I -^-(-Jt) first density 

^'^my={ ^ +{ ^ ^ 1 4096' (VIII.20) 



105 J_ _105 J_ 49 

128 ' 12 I 256 ' 20 I 1024 ' 



second density 



It should become now obvious from the numerical comparison of both approaches ('^''Vrtiii 

■^^^^777,2) that their corresponding perturbation orders are of equivalent nu- 
merical magnitude (~ 10^^), at least concerning the present lowest orders (/ = 1, 2). This 
suggests that none of both competitors is numerically more advantageous, as far as the 
rapid convergence of their series expansion is concerned. However, the basis elements 
C^2'}(^9) f |Viii.2|) of the second expansion f lVIII.Sp are more favourable from the purely 



technical viewpoint because, in contrast to their counterparts ^^^A^^^i!)) flV.25aP etc. due 
to the original expansion flV.22p . they are independent of the angular density ^''^ko^'d) and 
therefore the second approach appears to be more manageable for dealing with the higher 
orders of the multipole expansion (see below). 

2. Multipole Equations 

Once a new basis system of some promise is at hand now, one would like to see 
how the former quadrupole and octupole equations (lVI.43ap -f lVI.43bp transcribe to the 



corrsponding new radial potentials ^^^A2i{r) (IVIII.3p . The desired Poisson-like equations 



for ^*^A2/(r) may be deduced by inserting the new expansion f lVIII.3p of the anisotropic 
part ^^^A^^{r,i!}) of the gauge potential ^^^AQ{r,'d) fIV.lQp in the "anisotropic" gauge-field 
energy E^^ ( ]VI.13p and also in the "anisotropic" Poisson constraint N^^n (1VI.24p . Since 



the latter two contributes to the energy functional E[q] (]VI.25p are the sole ones which 



contain the anisotropic part ^^^A'^^(r,'d) it is sufficient to consider the anisotropic part 
''^'^'^']E[n] of the energy functional alone, i. e. 

^-l[^]=Ei:> + Ag)-Nt>, (VIII.21) 
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namely in order to deduce hereof the Poisson-hke equations for the new radial potentials 
^''^A2i{r) by means of the principle of minimal energy 



(5/^^">E[f,] = . (VIII.22) 

Here, Si is the functional derivative with respect to any potential ^^^A2i{r) (/ = 1, 2, 3, 4, . . .) 
where the derivative of the mass-equivalent M^an (1VIII.13P is immediately evident: 

6e M^'i = dr r^ ^'Ua^r) j ■ |-/ic \ ■ (VIII.23) 

Anisotropy Energy 

n{c} 



The case of the anisotropy energy is a little bit more complicated. It is 



the sum of a radial derivative term Ein^ and a longitudinal term Ein\ i- e. 

Ei^ = + Ei:^ , (V111.24) 
with the radial contribution being given by 

and its longitudinal counterpart by 



Inserting herein the new series expansion (]VIII.3|) of the anisotropic part ^^^A^^(r, yields 
for the radial energy Ei:^ f lVIII.25|) 



l,m\=l 

and analogously for the longitudinal part Ein^ (IVIII.26P 

m2i{r)-m2m{r) 

l,m\=l 



El£^ . f; /,,. ■f'^'-'-) . (VIIL28) 

, , , J r 
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The functional derivatives of both energy contributions are easily seen to look as follows 



5/Ean(r) = 




(VIII.29a) 
(VIII.29b) 



But with all three functional derivatives being at hand now, one can substitute them 
in the anisotropic part of the principle of minimal energy flVIII.22D which then reads 
explicitly 

5, {->E[o] = (l + Ag)) • _ Ag) . 6, (mWc^) (VIII.30) 

with the Lagrangean multiplier A^"* = — 2. 



General Multipole Equation 

This then yields the following Poisson-like equations for the radial correction potentials 
^''^A2i{r) (^^ multipole equations") 

A. i^U.(r) - f: i^,, ■ ^^^^ = -a. ■ ^ , (VIII.31) 

k=l 

where the transformed mass-equivalent parameter ^''^m'l is given by 

oo 

Wm; = 5^ C,;,^ ■ Wm, (VIII.32) 

k=l 

and the transformed matrix {-F/fc} by 

oo 

For instance, for the octupole approximation where the matrix {C;~^} is given by equation 
flVIII.10|) and {Fi,m} by flVIII.12D one finds for {F/ J flVIII.33|l 



r , . . 6 -12 , 

J = • (VIII.34) 

^0 20 
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Observe also that all the higher-order multipole potentials ^^^A2i{r) become coupled 
by flVIII.3ip to the spherically- symmetric amplitude Q{r)\ 



Quadrupole Approximation 

As a brief demonstration of the alternative perturbation expansion (1VIII.3P one may re- 
sort again to the quadrupole approximation. Here, some results have already been worked 
out within the framework of the original perturbation expansion (IV.22p . see Sect. VII; 
and it may be instructive now to see how the same results can be obtained also by use 
of the alternative expansion f lVIII.Sp . Once we thus have gained some confidence in this 
alternative approach, one may then tackle the octupole approximation (or even the higher- 
multipole approximations) which will then appear in a more pleasant shape than in the 
original perturbation formalism. 

The quadrupole approximation is the lowest order of perturbation for which the system 
of multipole equations flVIII.31|) becomes reduced to one single equation (/ = m = 1). 
Thus the multipole matrices {C;,™} (IVIII.10|) and {-F«,m} f!VIII.12|) become reduced to 



their components Ci^i = ^ and Fi i = hence the transformed matrix simplifies to 



Furthermore, the mass-equivalent parameter ^''^m[ ( 1VIII.32P becomes 

^ V = C^l ■ '''rn, = I . (^1) = t'^ , (VIII.36) 



cf. (lVIII.15ci|) -f lVIII.15bp . As a result of all these quadrupole specifications, the general 
multipole equation (IVIII.3ip becomes cut down to the quadrupole equation: 



16 ° r 

where the upper/lower sign refers to the first/second density ^%(^9) flVIII.15aD / flVIII.15bj) 



Already from this form of the quadrupole equation it is possible to conclude that both 

121 



quadrupole correction potentials ^^^A2{r) (for jo = 2, ^^^jz = ±1) due to botli densities 
(lVIII.15a|) and (]VIII.15b|) will differ only in sign, just as was the case for the correspond- 
ing quadrupole potentials ^''^A"'(r), see the table on p. [821 This difference in sign has 
already been shown to represent the origin of the missing of the dimorphism for the states 
with jo = 2, ^'''^jz ± 1, so that this effect must occur also in the quadrupole approximation 
of the present alternative approach. Furthermore, the quadrupole equations of both ap- 
proaches, cf. flVII.4p vs. flVIII.37p . are almost of the same shape and differ merely in the 
coupling parameters, cf. (]VI.17bl) and (1VIII.36P : 



flVI.lTM : 



'mm 



'em 



5 ^ (IV111.36P : 



{bl 



+ 



15 
16 

15 
16 



, first density fIVlll.lSap 
, second density (IVIII.lSbp 



(VIII.38) 

This similarity of both quadrupole equations then entails that their solutions must be 
proportional, i. e. 



'mm 



'em 



{b}^ 



III/ 



(VIII.39) 



Concerning now the energy being stored in this quadrupole mode, one specialises 
the general mass equivalent M^anC^ f lVIII.13p of the alternative approach down to the 



quadrupole appoximation (/ = 1), i.e. 



M^!>c2 



= -he ^^>mi / dr ^"^A.ir) 



(VIII.40) 



or by use of the link flVIII.39p to the original quadrupole potential ^''M"'(r) 

fi2(r) 



—he 



{b}^^ . {b}^'^ 



(6). 



'mm 



'em 



dr r2WA"'(r) 



(VIII.41) 



Here, the pre-factor is found by means of flVIII.361) and fIVI.lTbp 



{b}jn^ . {b}^'^ 



45 
4 



1 ^2 



(b), 



'mm 



'em 



1 

64 



(VIII.42) 
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so that the quadrupole mass equivalent in the alternative approach M2'^^c^ (lVIII.4ip 
actuall equals the original form M|^|'f^c^, cf. (jVL22a]): 

Mfii^c^ = ml'^c^ . (VIII.43) 

For the specific form of Mi^i^^c^ see equation flVII.lSI) . 

Of course, this identity of the mass equivalents must have its counterpart for the 
quadrupole energies, i. e. one concludes 

E},f = , (VIII.44) 

where the original form of E^^^ is given by equation flVII.12p . with ^^^em = (table 
on p. I90|) . and its alternative form E2^^ is the specialization of E^'^ (]VIII.24p to the 
quadrupole case (/ = 1), i.e. 

oo r 1 2 

^ 41... ^ Et^ ^-^J dr [c., (^^^) (VIII.45) 





2 



Inserting here again the required matrix elements from the matrices Ci^m (IVIII.lOp and 



Fi^m (IVIII.12P in combination with the tranformation of the quadrupole potentials f lVIII.39p 
does actually validate the claimed equality f lVIII.44p . For the precise form of E^^^ see 
equation flVII.12p . 

This complete equivalence of the original and the alternative perturbation approach 
(within the framework of the quadrupole approximation) should provide now sufficient 
confidence for tackling the next higher order of the perturbation expansion, i. e. the 
octupole approximation, just by means of the alternative technique. 
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3. Octupole Approximation 

Surely, it is reasonable to suppose that a certain improvement of the quadrupole results 
(Fig. VII. A) can be attained by taking into account also the octupole character of the 
charge distribution ^^^kQ{r, -(9), i.e. we consider now the multipole equations flVIII.3ip with 
restriction of the multipole order to l,m\ < 2. The required octupole matrices {Ci^m}, 
{^/m}' {^i,rn} and {-F/m} have already been specified by equations ( ]VIII.10|) . (]VIII.12p 



and (IVIII.34p . Thus, we obtain from the general multipole system ( IVIII.3ip for / = 1 the 

modified quadrupole equation 

A. .'U.(.) - 6 + 12 = .'W; . ^ (VIIL46) 



and for / = 2 the pure octupole equation 

A, ^'^A,{r) - 20 ^^^^ = -a, ^''W^ • ^ . (VIII.47) 

Any individual multipole equation of order I owns its specific coupling constant 
(1VIII.32P : especially for the present octupole approximation (/ < 2) we have for the 
quadrupole coupling constant in equation ( ]VIII.46p 



{b}m\ = C^l ■ ^^>mi + C^l ■ ^''>m2 =1 ' ' " (VIII.48 



105 
256 ' 



first density f lVlll.lSap 



525 
128 ' 



second density (IVIII.lSbp 



and similarly for the octupole coupling constant in equation (1VIII.47P 



, , first density fIVIII.lBap 
{H^'^ = C^l . Wm, + . {^im2 ^ ^ (VIII.49) 



III , second density f IVIII.lSb]) 



Observe here that, on account of the particular form of the matrix {-F/^} (IVIII.34p . 
the octupole equation (IVIII.47P expresses the decoupling from the quadrupole influence 
whereas the quadrupole equation (IVIII.46P still contains the octupole potential ^^^A^^r). 



Surely, it would be more aesthetic if one could deduce also a pure quadrupole equation 
which formally does contain only the proper quadrupole terms (see below). 
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Octupole Energy 



Concerning the energy being stored in the quadrupole and octupole modes, one may 
first turn to the mass equivalent M^anC^ (IV111.13I) which is the sum of the individual 
multipole contributions, without any interference term: 



^:>c2 ^ mI%c^ = Ml'^h^ + M\'h^ , (VIII.50) 
with the quadrupole contribution being given by 

oo _ 

M^'^^c^ = -he ^^>mi ■ [drr^ m^ir) 



(VIII.51^ 







and analogously for the octupole contribution 

m\'^c^ = -he ^'>m2 ■ J drr^ ^''U^ir) . (VIII.52) 







(For the mass-equivalent parameters ^*>mi and ^^>m2 see equations flVIII■15aj) - ^VIII.15bj) )■ 



Unfortunately, such a pleasant separation of the various multipole contributions does 
not occur for the anisotropy energy E^^ (]VI.13|) . For its radial part Ein^ ( ]VIII.27|) one 
finds 







^15 c/WA2(r) rfWA4(r) ^ fd^''Ui{r) 



and similarly for the longitudinal part -Ean ( 1VIII.28|) 



7 dr dr \ dr 



oo 

Eit^ E^ll = ^ {l^ i^'^Mr)Y + I ^'^Mr) ■ ^'^A,{r) (VIII.54) 
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But now that all contributions to the "anisotropic" energy are explicitly known one can 
of course carry through the variational principle of minimal energy (1VIII.22P for I = 1,2 



and would then regain the quadrupole-octupole system (lVIII.46p - (lVIII.47p . Though this 
result provides sufficient confidence in the consistency of the general multipole approach 
flVIII.31| ) it is nevertheless very instructive to consider also a slight modification thereof. 

Conjugate Potentials 



Reconsidering for a moment the radial anisotropy energy ^in^ flVIII.27D . one is tempted 
to introduce the conjugate potentials ^^^A2i{r) through 



oo 



^'^«(^) = Yl ^'^^2™(r) , (VIII.55) 

m=l 

or conversely 

oo 

WA2z(r) = ^C,;^WAL(r). (VIII.56) 



m=l 



In terms of these conjugate potentials, the former multipole equation (lVIII.3ip reappears 
now in the following form 



oo 



A. i^U',(r) -J2f;:, '-^^ = -as ^ . (VIII.57) 



fc=i ^ 



The matrix {F;^} emerging here is in general asymmetric and is defined through 

oo 

771=1 

For instance, for the octupole approximation (i.e. m,l,k\ = 1,2) this matrix looks as 
follows 

{F;:,}=i ^ I (VIII.59) 

y-12 20 J 

and thus is the transpose of the former matrix {-^/fc} flVIII.34p . Furthermore, the octupole 



approximation {l,k\ < 2) of the general multipole system flVIII.57p consists now of a pure 
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quadrupole equation (/ = 1) 

A, ^^U'2(r) - 6 ^^^^ = -as ^'^m, ■ ^ (VIII.60) 

and a modified octupole equation (/ = 2) 

A W^Kr) - 20 ^^^^ + 12 ^ = -a, ^'^m, ■ ^ , (VIII.61) 

which is just the other way round if compared to the original case (lVIII.46p - (IVIII.4"7]) . 

Of course, the conjugate multipole system (IVIII.57P can also be deduced from the 
principle of minimal energy 

5[ (^^"^Ejo]) = (VIII.62) 

where 5[ denotes the functional derivative of ^^"^^Epj with respect to the conjugate poten- 
tials ^^^A'2i{r) f lVIII.SSP : and "^^'^"Ej^j is the former functional ''^^'^Ejq] expressed in terms of 
the conjugate potentials ^'^^A'2i{r)^ cf. f lVI.25p . In order to validate this claim, one first 
transcribes the mass equivalent M^anC^ f lVIII.iap to its conjugate form: 

M^:>c2 ^ 'M^:>c2 = -he ^ ■ I drr'' ^'U'2,(r) . (VIII.63) 



oo 







Next, the radial energy Ei:^ flV111.27p reappears as 

oo 

he ^ r . odm'jr) dmi, 

l,k\=l 

and finally the longitudinal energy E^t^ (IVIII.28P adopts the following shape 

^i,k I '■ 

/,fc|=i 

with the symmetric matrix {-^/fc} being given by 



oo °° 

Eit' 'Eit^ = -^Y.n'Jdr '"'^-^^yf^'^'^^^^^ (VIIL65) 



K'k= ^n} ■ (VIII.66) 

m,n\=l 
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Now one can carry through again the variational principle (1VIII.62I) (see (]VIII.30| ) for 
a more explicit form of this) in order to find the corresponding extremal equations just 
coinciding with the result (IVIII.57p . as expected. 

For a concrete demonstration of the conjugate formalism, one resorts again to the 
octupole approximation {j,k,l,m\ < 2). Here, the mass equivalent 'M^anC^ flVIII.63j) 
appears again as the sum of the quadrupole and the octupole contributions 

'M^^^c^ ^ 'mI%c^ = 'ml'^c^ + 'm\'^c^ . (VIII.67) 
The conjugate quadrupole part is given here by 

^ {b}^,^ J ^2 WA'^^r) ■ (VIII.68) 





'^2 



and the corresponding octupole part by 

oo _ 

'Mf = -he ^'Wa / dr ^%(r) ■ . (VIII.69) 







Next, the radial energy 'Ean flVIII.64p becomes specified down to 



oo r, 



^^^rfWA'2(r) d^^%{r) 245 /d^''U'^{r] 
dr dr 4 \ dr 

And finally, the longitudinal energy 'Eit^ f lVIII.65p is found as 

-^an ^ -^2+4 ~ 

45 he 2 924 (WA'2(r))' - 2100 ^'^A'^ir) ■ ^^U^(r) + 1225 (^^^(r))' 



16 as Jo ' ' r^ 

(VIII.71) 



For the desired deduction of the octupole system (IVIII.60p - flVIII.6ip by the extremal- 



ization of the "anisotropic" energy functional ''^'''^Ejq] flVI.25p we have to write down now 
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that functional in the octupole approximation again in terms of the conjugate potentials 
WA'2,(r), i.e. 



2+4 -^2+4 1 -^2+4 ^ -^2+4 ^*-^2 ^'^M 



By means of the standard variational technique, it is then easy to see that both extremal 
equations 

5'^ ^2+^>E'[f^] = (VIII.73a) 
^4 ^^^^^^'m = (VIII. 73b) 



actually do lead us back to just the claimed octupole system ( IVIII.60p -f lVIII.611) ! 

Hybrid Method 

There exists also a kind of hybrid formulation of the octupole system which then 
contains both kinds of potentials ^^^A2i{r) and ^^^A2i{r). Indeed, it is easy to see that the 
radial energy E^} flVIII.271) can also be written in the following form 

oo 

(vni.T4) 



=1 







and similarly for the longitudinal contribution £'an (]VIII.28| ) 

oo 



4„^^ Eil^ ^-^j:nJdr r 

' k,l\=l 



(VIII.75) 



where the matrix {-F/'^} has already been specified by equation fIVIII.SSp . Concerning the 
mass equivalent, one takes half the sum of its two forms given by ( 1VIII.13|) and ( 1VIII.63|) . 
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I.e. 



tW^2 , mW^2 



he 



oo ^ 02 



hy 



1 = 1 



(VIII.76) 

The corresponding hybrid form of the energy functional (1VI.25P is 



(VIII.77) 



and its extremahzation, cf. (1VIII.30P and (IVIII. 73a|) - ( IVIII. 73bp . with respect to both 
potentials ^''^A2i{r) and ^''^A2i{r) leads us again back to both former multipole systems 
( IVIII.46|) -f lVlTL47l) and f lVIII.60|) -f lVlTr6T|) . 



* * * 



At this point, it seems worth while to pause for a moment in order to briefly survey the 
results of the octupole approximation obtained so far. Perhaps the most interesting feature 
is that one can work with two alternative formalisms, namely the one being based upon 
the alternative multipole potentials ^''^A2i{r), which were introduced in connection with 
the new basis system, cf. (IVIII. 3p . or one may work with the conjugate potentials ^'^^A2i{r) 
(IVIII. 55p . Clearly, both formalisms are completely equivalent and one may prefer the first 
or the second one. Indeed, both approaches lead us to an octupole system of equations, 
i.e. in the first case f lVIII.46p -f rVIII.47p and in the second case flVIII.60p - flVIII.6ip . Each 
one of both systems consists of a quadrupole equation, ie. (1 VI 1 1. 46 P and (IVIII. 60p . and of 
an octupole equation, i.e. flVIII.47p and flVIII.6ip . Furthermore, both approaches admit 
the existence of a variational principle and thus are equipped with an energy functional. 
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^^^^J'Ein] and '^^^^'^-''E'j^j (]VIII.72p . whose extremalization yields just the corresponding 
octupole systems. 

But a difference (albeit only in a technical respect) exists concerning the specific 
shape of the multipole equations: whereas in the first case (1VIII.46P - (IVIII.47|) the sec- 
ond equation flVIII.47p for the octupole potential ^''^A^^r) appears decoupled from the 
quadrupole potential ^^^A2{r), the first equation flVIII.461) for the quadrupole potential 
(r) contains a coupling to the octupole potential ^^^y44(r), i.e. the octupole potential 
(r) appears to be independent of the quadrupole potential, but the converse is not 
true. This quadrupole-octupole interrelationship is just the opposite of the second case 
f lVlll.GOp - dVlll.eip : here the quadrupole potential ^^U'a (r) seems to be independent of 
its octupole counterpart ^''^A'^(r), whereas the latter is coupled back to the quadrupole 
potential ^^U'^{r), cf. f lVIII.Gip . 

Separative Method 

This coupling between the quadrupole {^''^A2{r), ^'^^A'2{r)) and octupole {^''^Ai{r), 
^^^A'^{r)) potentials cannot be removed by resorting to the hybrid method since the lat- 
ter does embrace both coupled octupole systems, ie. (]VIII.46P - (lVIII.47p and ( 1VIII.60P - 
(]VIII.6ip . However, it would be very desirable to have a completely decoupled system 



because in this case it would become possible to associate well-defined anisotropy energies 
to both the quadrupole mode (/ = 1) and to the octupole mode (/ = 2). Indeed, only 
for such a situation it would be legitimate to talk about a "quadrupole energy" and an 
"octupole energy" . But such a mode separation of the anisotropy energy is not possible up 
to now because the anisotropy energy Ein = -^2+4 + -^2+4 of three methods do contain 
a mixed term consisting of both the quadrupole potential ^^^A2{r) /^^^A2{r) and the oc- 
tupole potential ^''Ui{r) /^^W^{r), cf. equations (IVIII.53D - (IVIII.54p . flVIII.70p - (IVIII.71D . 
and (IV111.74l) - (IV111.75p . 

Fortunately, there exists a possibility for the complete separation of the quadrupole 
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and octupole modes, namely by introducing suitable new potentials ^''^A2{r) and ^''^A'l(r): 

Wy4'2'(r) = ^^U'2(r) (VIII.78a) 
^'U'iir) = WA4(r) . (VIII. 78b) 

This says that the new (separative) quadrupole potential ^^^A2{r) is identified with the 
conjugate potential ^''^A2{r) (IVIII.55P and the new (separative) octupole potential ^''^A'l{r) 
is identified with the alternative octupole potential ^^^A^{r) being formerly introduced 
by the transition fIVIII.Sp to the new basis system IclJ'-^^l. Since the new separative 
potentials are composed of both the alternative and the conjugate potentials, the present 
separative method is obviously the true hybrid method. Clearly, the separative potentials 
may be expressed also solely in terms of the alternative potentials ^''^A2i{r) 

WA'2'(r) = ^ • i^U2(r) + ^ ■ WA4(r) (VIII.79a) 
WA;'(r) = ^'U^{r) (VIII.79b) 

or soleley in terms of the conjugate potentials ^''^A2i{r): 

{''U'2'(r) = ^''U'2(r) (Vlll.SOa) 
^'^A'iir) = . WA'2(r) + 1^ . i%(r) . (Vlll.SOb) 

Conversely, the alternative potentials ^^^A2i{r) read in terms of the separative potentials 
^^U'^iir) as 

WA2(r) = ^ ■ WA'2'(r) - ^ ^^U:[(r) (Vlll.Sla) 
^'^A^r) = ^'^A'iir) , (Vlll.Slb) 

or similarly, the conjugate potentials ^''^A'2i{r) do reappear in terms of the separative 
potentials ^^U'^i (r) as 

{''U'2(r) = ^^U'2'(r) (VIII.82a) 
i^U;(r) = ^ . i''U'2'(r) + ^ . i^U:(r) . (VIII.82b) 
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By means of these results, one can now eliminate the alternative (or conjugate) po- 
tentials in favour of the separative potentials so that the corresponding systems, i. e. the 
"alternative" system flV111.46l) - (IV111.47p or the "conjugate" system flV111.60p - (IV111.6ip . 
are both recast, in terms of the separative potentials ^''^A2{r) and ^''^A'^{r) to the following 
form: 

A, ^'^A'^{r) - 6 ^^^^ = -^^Wi ■ a. ^ (VIIL83a) 
A, ^^U:(r) - 20 ^^^^ = -Wrn'^ ■ «s ^ • (VIIL83b) 



For the values of the mass-equivalent parameters and ^^Wg see equations flVIII.16P 
and (IVIII.49p . and observe also their relationship 



Wm'^ = 1^ (^{b}^, - ^ {^}m,^ . (VIII.84) 

Thus the separative method takes over both the pure quadrupole equation flVIII.60p 
and the pure octupole equation (]VIII.47p . inclusive the corresponding mass-equivalent 



parameters. Evidently, our attempt (lVIII.78a|) -( IVIII.78bp of decoupling the octupole 
systems has thus been successful; and consequently we can now look for the separative 
form of the energy functional so that the present octupole system flVIII.83ap - flVIII.83bp 
represents just the corresponding extremal equations thereof. 

Separative Energy Functional 

Prom the purely formal viewpoint, the desired "separative" energy functional (^^'^^E'^j, 
say) must of course be of the same shape as its "conjugate" and "hybrid" predeces- 
sors ^^+^^]Ejf^] flVIII.72p and ^''y^E[n] flVIII.77p : i.e. one expects again the splitting of the 



anisotropy energy in a radial and a longitudinal part: 

"Ei^ "e\:1 = "e\:1 + "e\'} , (VIII.85) 
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just as was the case with all the predecessors, e. g. (1VIII.24P for the alternative choice of 
basis system. Thus the "separative" functional ■"^^''^JEj'f^j will appear again as (cf. ■'^^"'^^^]E|j^j 
(IV111.72P and ^^y^E^n] flV111.77p ) 



{2,%// ^ 



2,4 • 



(VIII.86) 



But the pleasant property of this functional is now that it decays in a sum of subfunc- 
tionals, i. e. the quadrupole and octupole contributions: 



(VIII.87) 



where the quadrupole contribution is given by 



(VIII.88) 



and the octupole contribution by 



(VIII.89) 



Here, the individual energy contributions are obtained by simply inserting those in- 
verse transformation relations flVlll.81ap - flVIII.81bp or flVlll.82ap - flVIII.82bp in the corre- 
sponding "conjugated" or "separative" results; this yields for the quadrupole contributions 
flVIII.88p 



"^{^} _ _45 nc 



4 tts 



dr r 



dr 





oo 

-^Wm^.^c J drr^ ^'^A'^ir) 



QHr) 



(viii.goa) 
(viii.gob) 
(viii.goc) 
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and similarly for the octupole contributions (1VIII.89P 



"EP = -—— ■ — drr^ { (Vlll.gia) 



11025 a, J V dr 





oo 



"i^r = ■ f / 2° (^^^^'4'(o)' (viii.9ib) 



• ^''Wa y dr ^'^A'Hr) . (VIII.91c) 



11025 







Such a pleasant decay of the energy functional (as being given by equation flVIII.871) ) 
into pure quadrupole and octupole contributions ^^lEj'^j and ^^^Efj^p cf. flV111.88l) - flV111.89p . 



does not occur for the other approaches, i.e. "alternative", "conjugate", and "hybrid". 
And the consequence is that any subfunctional is sufficient in order to deduce the cor- 
responding multipole equation. Indeed, it is easy to see that the extremalization of the 
quadrupole subfunctional ^^^Ej'j^j (1VIII.88P yields just the quadrupole equation (lVIII.83ap . 
and analogously the extremalization of the octupole subfunctional ^^^Ej'^j flVIII.891) yields 
the octupole equation flVIII.83bl) . This neat separation of the quadrupole and octupole 
effects can now be used in order to clarify the question of octupole dimorphism for those 
states with quantum numbers uq = Jo + 1 = 3, ^^^j^ = ±1 which in the quadrupole 
approximation misses this effect of dimorphism (see the discussion of this effect on p. [82]) . 

4- Magnitude of Octupole Splitting 

For a closer inspection of the critical energy level {no = 3,jo = 2,Wj, = ±l}it IS very 
instructive to first consider the quadrupole part flVIII.83aP of the octupole approximation 
f lVIII.83ap -f lVIII.83bp ! 
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Quadrupole Pari 



Observe here, that this is not identical to the former quadrupole approximation 
flVIII.37p . being due to the "alternative" approach, since the coupling constant is now 
^^>mi flVlll.15ap -f lVIII.15bp in place of ^^^m[ flVlll.361) . But a common property is the 



fact that the quadrupole ambiguity shows here up again in form of the sign ambiguity of 
the mass-equivalent parameter ^^Vni flVIII. 1 5a]) - flVIII. 15bp . just as is the case with 
flVIII.36P : and this ambiguity transcribes then also quite correspondingly to both solutions 
^^^A2{r) and ^^^A2{r). The consequence of this is that, also within the framework of the 
separative method, there can occur no dimorphic energy splitting of those states specified 
by the quantum numbers {no = 3, Jo = 2, ^''^j^ = ±l}. The reason is again that the cor- 
responding anisotropy energy "Ef ^ "E^"'^ + "E^^^^ flVlll.90ap - flVIII.90bD is quadratic 
with respect to the quadrupole potential ^''^A2{r), and is therefore insensitive with respect 
to a change of sign of ^''^A2{r)\ 

A further common property of the former quadrupole approximation flVIII.37p and the 
present quadrupole part flVIII.83al) is the fact that both approaches do yield the same 
quadrupole energy, i.e. we have not only the energy identity f lVIII.44p but also 

E^;> = "E^"^ (VIII.92) 

where E^^^ has been defined in terms of the original potential correction {^^^'"(r) by equa- 
tion flVII.12p and the present energy "E^"^ (= "E^ ^ + "^2^''^) by the equations flVIII.90ap - 



flVIII.90bp . In order to validate this claim, simply observe that when i^iA"'(r) is the so- 
lution of the former original quadrupole equation flVII.4p then the solution of the present 
quadrupole part f]VIII.83a[) is 

{''U'Vr) = - ^^>mi ■ W^"'(r) . (VIII.93) 
5 

On the other hand, the quadrupole part f lVIII.83a|) of the octupole approximation is 
uniquely linked to the quadrupole energy "E^''^ (= "e!{^ + "E^^A flVIII.90aD - flVIII.90bp 
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so that we merely need to substitute therein the solution ^''^A2{r) (]VIII.93| ) in order to 
validate the claimed identification (]VIII.92p . 

Octupole Energy 

Finally, let us mention also the fact that the decoupling of the quadrupole and octupole 
potentials ^^^A2{r) and ^''^A'l(r) by the separative method, as expressed by (lVIII.83aP - 
flV111.83bp entails the separate validity of the octupole identity 

"eJ^^ = "m\'^c^ , (VIII.94) 



where the octupole energy "e\^^ ^= "e\^^ + "E^^^^ is given by equations flVlll.QlaP - 



IVIlI.91b|) and the corresponding mass equivalent "Mf ^c^ by f lV111.91cl) . If we are satisfied 



here again with the former trial amplitude Cl{r) fIVI.lap - flVl.lbp . both objects flVII1.94p 
become functions of the trial parameters /3 and u, i. e. 

{2/3aB) ( ' ■ 4'^ (J^) (VIII.95a) 
{2(3aB) ( ' ■ fif^ (u) (VIIL95b) 

') are explicitly determined in App.E. Thus, the 
octupole identity ( ]VIII.94|) reads in coefficient form 

4^^(z/)=/if^(z/). (VIII.96) 







e^ 






as 








"IV 


li'^^c^ = 


e^ 






aB 




[u) and 


W/ 



The latter results ( 1VIII.95a|) - (]VIII.95b|) can now be used in order to convince oneself 



that the absence of the energy splitting becomes eliminated by the octupole approxima- 
tion: namely, the mass-equivalent "M4''-'^c^ (]VIII.95bp shows that this object is quadratic 



with respect to the mass-equivalent parameter ^''W2, i.e. 

"mfc'' ~ (^'Wa)' , (VIII.97) 
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because jj^l^^u) is independent of 'f^Vrig, see equation (1E.27P of App.E. On the other hand, 
it has aheady been shown through the equations (1VIII.49P that ^''Wg must adopt two 
different values which differ by a factor of six due to the existence of two different angular 
densities ^''^/co('*9), see equations fIVIII. 15ap - flVIII. 15bl) . Therefore the octupole correction 



energy flV111.95al) - flVIII.95bp must lower the gauge field energy of the considered ortho- 
state (with no = Jo + l = 3, ^''^j^ = ±1) in two ways and thus creates the dimorphic energy 
splitting for this state! The corresponding two (energy lowering) octupole-corrections will 
differ by a factor of six, see below. 

Schro dinger Equation with Octupole Interaction 

Properly speaking, the Schrodinger-like equation for the amplitude field Cl{r) has 
already been deduced from the energy functional E[n] flVI.251) and yielded the equation 
(jVI.27p . However, that deduction had been based on the original decomposition flV.22p 



of the anisotropic gauge field ■t''M^"(r, so that the corresponding Schrodinger equation 
(1VI.27P refers to the potential corrections ■>^^^A"'(r) and ^^^A^(r), etc. However, in the 
meantime it has turned out to be more advantageous to decompose the anisotropic gauge 
field ^''^A^^{r,'d) as shown by equation (1VIII.3P where the alternative potentials ^''^A2i{r) 
have been further transcribed to the conjugated gauge fields ^''^A2i{r) and ultimately to the 
separative potentials ^^^A'^{r) and ^^^A'l{r) fl VIII. 78aD - fl VIII. 78bD . But since the separative 
method provided us with a clear decoupling of the quadrupole and octupole interactions, 
it will surely be helpful to write down the Schrodinger equation for the amplitude field 
also in terms of the separative potentials ^^^A2{r) and ^^^A'l(r). 

For this purpose, we have to go back to the original energy functional K^q] (]VI.25p 
and have to recast the anisotropic gauge field part ^^°j]E[q] in terms of the separative 
potentials ^^^A2{r) and ^^^A'l{r), whereas the matter part and the "isotropic" gauge 

field part I^fep] remains unchanged. The latter fact then implies also that the Poisson 
equation flIV.181) for the isotropic potential ^''^Ao{r) remains unchanged. Thus, it is only 
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the Schrodinger equation (]VI.27|) which must get a new shape. 



The latter is to be deduced from the energy functional (1VI.25|) by extremalization with 
respect to the amplitude field (l{r), which is contained only in the matter part ^^5E[q] (=f 
2 Ekin + 2 A, ■ Nn) and in all the mass equivalents M^c^ d\rT2|) . "mJ^^c^ (IVlll.QOcp . and 
"M4'^-'^ flVIII.91c|) . Consequently, denoting the functional derivative with respect to the 
amplitude field Q{r) by 6q, the corresponding extremal equation reads in abstract form 

= Sn %n] = Sn {^""^n]) - Ag^^f, (m^c^) (VIII.98) 
-X'g>6n ("mJ^^c^) - Ag^^f, ("Mf c^ 

or concretely 

(iji + \iz]^ir) + -^.jlm (VIIL99) 



2M ydr"^ r drj ' ' 2Mr'^ 
-he jl^Uolr) + ^ Wmi ■ ^^U'2'(r) + ^ ^^W^ ■ WAl'(r)| r2(r) 

= -A. ■ fi(r) , 

where the Lagrangean parameter A^ plays the part of the Schrodinger energy eigenvalue 
-E* (= —As). Neglecting the octupole interactions {^^^^[{r) =^ 0) and observing 

^ ^^>mi ■ WA'2'(r) ^ ("^miii ■ WA"'(r) (Vlll.lOO) 

lets the present octupole Schrodinger equation (IVIII.99P becoming identical to the former 
quadrupole Schrodinger equation flVI.27p . 

Variational Procedure 

Unfortunately, we are not smart enough in order to find the exact solution of the 
eigenvalue problem flVIII.991) which must be complemented by the monopole equation 
fllV.lSp for the spherically symmetric approximation t^Uo(r), by the quadrupole equa- 
tion flV111.83al) for ^^^A'/^{r) and finally by the octupole equation flVIII.83bD for ^^^^'(r). 
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Therefore we must be satisfied here with the extremahzation of the energy functional 
IE[a] f lVI.25p on the subspace of the trial amphtudes Q{r) (]VI.la|) - (lVI.lb|) and its descen- 
dants [%(r) a^{y) dVL30|) ). ^''^A'^ir) ^^UJj'd/) (lylLS]) ). and ^^^'(r) ^^^4(1/) 
(IE.26P ). Each of these spherically-symmetric multipole potentials carries its own energy 
content, i.e. 



flVI.32ajl : 
flVlLT2|) : 
flVIII.89|) : 



[%(r) 



{e} 
III 



(2/3aij) ■ £:potl 

^2 



(VIII. 101a) 



"E. 



— (2/3aB) ■ ^'^mii, ■ /if,f (z/) (Vlll.lOlb) 

ttB 



(VIII.lOlc) 



(Observe here that the quadrupole energy E^f = "E\'^ flVIII.101b|) appears to be pro 



{e} 



portional to the mass-equivalent parameter ^^\n\\\ whereas the octupole energy "El 
(IVIII.lOlcP is quadratic with respect to the corresponding parameter ^^Vn'g- The rea- 



son is that the universality of the ratio 



'mm 



(6), 



'em 



5 (]VI.17bp admits to absorb this into the 



quadrupole ptential ■f''M"'(r), cf. the quadrupole equation flVII.4p : see also the discussion 
in Appendix A. 

The total gauge field energy E^^ f lVI.12p in the octupole approximation is now due 
to the sum of all three multipole components (IVIII.lOlap - llVIII.lOlcp . i. e. 

^2 



E, 



{e} 
R 



aB 



(2/3as)-4o1(^) 



(VIII.102) 



with the total potential function e^t (z^) being defined through 
4ot(^) = eM^) + (^W,M ■ /if,f (z.) + (i^W,)' ■ /if (z.) 



(VIII.103) 



Recall here also that the functions epot(z^) (1VI.33I) . fJ'l^\^{i^) (IVII.lSp . and ^"^^ (u) (IE. 271) are 
independent of the quantum numbers no, jo, ^''^jz^- Furthermore, the kinetic energy Ey^^^ 
f lVL8|) is still given by flV1.9ap - dVL9b|) : and the three poisson identities f lVI.34p . (IVII.lip . 



and flVIII.94p are satisfied since we always use the exact solutions of the corresponding 
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multipole equations. Therefore these constraints can be omitted in the energy functional 
(1VI.25P together with the normahzation condition (IVI.SP : = 0, provided that 
we restrict ourselves to the use of normalized trial amplitudes Q{r). This cuts down the 
energy functional ]E[q] (1VI.25P to its physical part 

eI^^^ = 2 Ekin + 4' + e};;^ + , (VIII. 104) 

and the value of this reduced functional on the subspace of RST fields {^(r); [^Uo(r), 
^^^A2{r), ^'^^A'Kr)} due to our trial ansatz fIVI.laP - flVI.lbp yields then an energy function 
&'^^{I3, v) which depends on the trial parameters /3 and v. 

E{iv}(^^ ^) = ^ i2paBf ■ 5kin(^) - — mas) ■ e\i{v) , (VIII.105) 

where £:kin('^) is given by fIVI.Qbp and ^tot (z^) by flVIII.1031) . 

The quantum number jo (= no — 1) is contained explicitly in the kinetic function 
£kin(i^) (lVI.9b[) and implicitly in ^tot (z^) (1VIII.103P via ^''Vn\\\ and ■f^Vri2, whereas the az- 
imuthal quantum number ^'^^j^ is (implicitly) contained only in the potential function 
^tot l'^) (1VIII-103P . But after having fixed both quantum numbers jo (= no — 1) and 
^^^jz with —jo < ^^^jz < Jo, we can obtain the non-relativistic energy spectrum of ortho- 
positronium by simply looking for the minimal values of the energy function ¥.^^"^^(13, v) 
(|VIII.105p with respect to both trial parameters /3 and v. By means of partial extremal- 
ization (i.e. with respect to /3), there remains to minimalize the reduced energy function 
E^Q^iy) with respect to the residual parameter v 

e'^\u) = . S'Shu) , (VIII.106) 

i.e. the desired energy value E^'^, associated to the quantum numbers {j} = {jo-, ^''^j^}, 
is given by the maximal value of the spectral function Sq^u): 

Eg^ = -^^gLx- (VIII.107) 
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Here, the spectral function Sq^u) contains the relevant quantum numbers {jc^^'^Jz} in 



the following way 



'tot y 



1 _L J^O 
2u+l \ 4^ 2u 



1 



(VIII. 108) 



The simple potential function £:pot(^) due to the spherically symmetric approximation 
is given by equation f lVI.33P : its octupole analogue e^^^u) ^= /i4°^(z/)j is also explicitly 
prepared in App.E, cf. f lE.27P : and thus it is only the quadrupole mass-equivalent function 
/i2°^(z/) which necessitates a brief comment. However, the quadrupole approximation 
has already been studied thoroughly in the preceding sections; and if we identify the 



corresponding mass equivalent m|-|'|-'^c^ (1VII.13P with its separative counterpart "M2''^c 



W.2 



fIVIII.QOcI) one is immediately led to the identification of the mass-equivalent functions: 

^'Un ■ f^l:h^) = i^'^^i)' ■ f^2\^) ■ (VIII.109) 

But from here it is a simple exercise to find the direct link of both mass-equivalent 
functions: inserting the general form f lA.ip of the quadrupole potential ^''^'"(z/) in the 
definition flVI.22ap of the mass equivalent ^^^m\\\ yields 



III = am ■ ■f''^mi . 



(VIH.llO) 



Next, observe equation (1A.4I) and conclude from this 

4 



which then finally yields 



and thus 



mi = - am 
y 



9 



Wmm = ^ (i^WO' 



{c} / \ ^ {c} / \ 
f^2 i'^) = ^Kll l'^) 



(VIH.lll) 



(VIH.112) 



(VIH.113) 
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So we can ultimately express the spectral function (u) in terms of the three known- 



If desired, the three functions epot{i^), fi}ll (i^), fJ^^ (i^) may be alternatively used here in 
their closed analytic form, cf. (lD.3ap - (lD.3bp and flE.49p . The subsequent table (see p. ll45p 
represents a collection of the numerical results: 

i) The groundstate {jo = 0, first line) is not doubled (App.s A and B) and re- 
ceives a relatively bad numerical prediction. The conventional value, cf. fll.isp . 
is E^j^ ~ —6, 8029 [eV] whereas the present RST octupole approximation yields 
—7, 59507 [eV] which is a deviation of —11,6%. Why the groundstate has such 
a large negative deviation in contrast to the positive deviations of the excited 
states {jo = 1,2,3,.. .) must be clarified elsewhere. For instance, the first excited 
state jo = 1, ^^^jz = ±1 (second line) has a deviation of +7,4%. 

ii) The exceptional case {no = 3, jo = 2, ^^^jz = ±1}, which misses the level dou- 
bling in the quadrupole approximation (see the discussion of the table below equation 
(IVII.22P and App.D), is found to undergo now also the level doubling in the octupole 
approximation (third line). The energy difference of the emerging dimorphic pair is 
found as 0, 68566 - 0, 68124 ~ 0, 00442 [eV] which amounts to (roughly) 0,65% of 
the binding energy of the dimorphic partners. The corresponding energy difference 
for the other two dimorphic pairs, i.e. jo = 3 and jo = 4, show up in the same 
order of magnitude (1,65%). 

iii) It is also interesting to see how the predictions change when one passes beyond the 
quadrupole approximation (fourth column) to the octupole approximation (sixth 
column). As a convenient measure of this change one may introduce the relative 



functions epot('^) (jYLSSD, Ai^iV(z/) flVll.151) . and uriiy) i^M)- 




epot(^) + (I WmO' ■ /if,?(^) + (^^W,)^ . /.f (z.) 



(VIII.114) 
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quadrupole-octupole deviation through 

A{^> = ku^^^ loc ^ (VIII.115) 

see the seventh column of the subsequent table. Here it is interesting to note that the 
quadrupole predictions (fourth column) become corrected at most by (roughly) 0,5% 
or less. This may be understood as a signal that our separative perturbation method 
can be expected to converge very rapidly. In any case, the level splitting due to the 
ortho- dimorphism seems to be at least twice the present octupolc corrections; and 
therefore the doubling of the ortho-levels can hardly be interpreted as an artefact 
occurring only in the low perturbation orders (i.e. quadrupole and octupole). Of 
course, this question should be studied in greater detail by an extra treatment. 
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Jo = 


Jz 




e}~J) [eV] 


( ^b>m'2)2 


E\f [eV] 




4°) 


{no - 1) 


flVI.22a|). flVIII.112D 


Quadrupole, (1VII.22I) 


(1VIII.84D 


Octupole. flVIII.107) 


(1VIII.115!) 


(11.181) 








^ = 0,0625 


-7, 53786 


(^)2^ 30, 6595 


-7.59507 


-0, 75 % 


-6,8029... 


1 


±1 


^ = 0,015625 


-1,57527 


(11)2^ 0,3785 


-1.57551 


-0,015% 


-1,7007... 


2 


±1 


= 0,015625 


-0.68112 


(fif)2 ^ 0,3785 


-0.68124 


-0,017% 


-0, 7558 . . . 






(11)2- 13,6264 


-0.68566 


-0, 66 % 


3 


±1 


= 0,015625 


-0.38089 


(11)2- 13,6264 

^ ZOO ' ' 


-0.38374 


-0, 74 % 


-0,4251... 




0,03515 


-0.39001 


(li)'- 0,8516 


-0.39019 


-0, 046 % 




±1 


0,03515 


-0.25010 


0,8516 


-0.25023 


-0,051% 


-0,2721... 


4 


(1^)2- 0,04515 


-0.25330 


(ii)'=^7,664 


-0.25443 


-0, 44 % 



Appendix A 



General Properties of the Quadrupole 
Approximation 



For getting some survey of all properties of the quadrupole approximation it is most 
instructive to consider also a further example, i. e. the case jo = 4. The subsequent table 
(at the end of this appendix) presents a collection of the corresponding relevant data for 
jo = 4 and thus complements the precedent tables as our basis of inductive reasoning. 



General Validity of the Ratios ( [ VI. 1 7a\ ) - (jVI-l 71^ 



The first and most striking observation refers to the fact that the ratios (]VI.17ap - 
(IVI.lTbp are valid for all considered cases ranging from jo = 1 up to jo = 4. Thus 
the self-suggesting conclusion from this limited number of cases is that these ratios could 
perhaps be valid for all jo (i- e. jo = 1, 2, 3, 4, . . . ). If this conclusion could be shown 
to be true, the quadrupole equation (IVII.4|) for the potential correction ■t*iA"'(r) would be 
generally valid, as well as the general form of the quadrupole energy e};;^ f lVII.12|l . 

Indeed, it is a rather simple matter to validate the first one of these ratios, i. e. 
f lVL17a|) . For this purpose, observe first that the general form of the angular potential 
factor flV.25ajl is 



am ■ 



cos {} 



(A.l) 



m^iii(^) 

with the constant am being given by 

am = -^ (l-3ir{y , (A.2) 
see equation (lV.28p . Next, calculate both constants ^''•'em and ^'^^ f\\\ (lVI.16aP by use of 
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the present form of ^^^A^^\i!}) dA.lh and find 



We,,, = ^ ■ af„ (A.3a) 
("Vm = ^ ■ «f„ • (A.3b) 



Consequently, this verifies the first ratio flVI.17ap . Observe here that this ratio is a purely 



geometric relation, independent of the specific angular density ^^^kQ{'d). 

A quite similar (albeit somewhat more complicated) reasoning yields for the constant 
('')m,ii (lVI.22aD 



^'W„, = ( ^ai„ ) (A.4) 



which then validates the second ratio flVI.17b[) . Thus both ratios f lVI.17ap - flVI.17b|) are 



found to represent quite general features of the quadrupole approximation and thus turn 
out to be independent of the chosen trial amplitude Cl{r) fIVI.lap . It is also a rather easy 
exercise to verify that all states with Jq = 4 (see end of this appendix) actually do satisfy 
the claimed relations (1 VI . 1 7aP - ( I VI . 1 7b"|) . Especially, these relations must hold also for the 



exact (unknown) solution of the eigenvalue problem in the quadrupole approximation. 

No Dimorphism for ^^^jz = 

The subsequent table for jo = 4, p. 1152^ suggests in combination with the precedent 
table on p. [TO] that for ^^^j^ = and all values of jo (i-e. jo = 0, 1, 2, 3, 4, ... ) there 
does occur no dimorphism; and furthermore the angular density (1IV.20P is always 

given by 

^%(^) = ^ • , (A.5) 
27r^ sm v 

so that the angular potential factor ■^^^'"(■(9) (IV.25aP must always appear for = as 

1" 



m^iii(^) 



cos^ 'd — 
3 



(A.6) 



Clearly, once the claim (lA.Sp for the angular density ^^^Q^-d) is verified, then the result 
(Ol) for the angular factor ^''U"'(?9) is an inevitable consequence thereof, see equations 
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(]V.28p - (]V.29|) . Therefore, we merely have to show that the claim (]A.5|) for ^^^jz = is 
correct. 

However, the desired proof is very simple: First, specify the original angular system 
(lIV.6ap - (]IV.6bp down to ^''^jz = in order to find 



dgni^) „„,„Q „ .2 



- cot ^9 ■ gR{i») = TJo sin' ^ ■ 9si^) (A. 7a) 



— -— + cot ?9 ■ gs{^) = ±jo ■ ■ A.7b 

dv sm V 

This system is to be understood as actually representing two different coupled systems 
according to whether one takes the upper (u) or lower (1) signs. Let the solution of the 
upper case be denoted by ^^^gni^), ^^^gs{i^) and the solutions of the lower case by ('^(/^('i?), 
^^^gs{i^)- Then it is found that the only difference between both sets of solutions refers to 
a change of the signs, i. e. 

^'W(^) = -^'^k(^) (A.8a) 
^%{^) = ^"W(^) ■ (A.8b) 

However, as the very definition of angular density ^^^kQ{'d) (lIV.20p says, such a difference in 
sign does not imply a difference in the angular density itself; and this then entails also that 
the angular potential ^''^'"(t?) is invariant under that change of sign (IA.8ap -f lA.8bp . Such 



a change may be brought about also in the first-order system (IA.7ap - flA.7bp by simply 



putting jo =^ ~io- Finally, the Schrodinger-like equation f llV.17p remains unchanged 
by the latter replacement {jo =^ —jo)] and thus the solution Q{r) cannot react to it 
either, nor the solution ^^^Ao{r,'d) of the Poisson equation (lIV.18p . The result is that the 
physically relevant objects are insensitive to the change of sign (lA.8ap -( lA.8bp and thus 
all the quantum states due to ^''^jz = are physically unique (-w no dimorphism) . 

Finally, it remains to explain why all these states (with ^''^j^ = 0) share the same 
angular density ^''^ko('d) (lA.Sp . independently of the quantum number jo- To this end, 
reconsider the definition of angular density ^^^kQ{i!}) flIV.20p and find by use of the present 
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first-order system (lA.7a|) - (]A.7b|) the differential equation 



+ cot^9-^%(i9) = . 



(A.9) 



Of course, the solution hereof (being normalized to unity according to the prescription 
(lIV.24a[) ) is just identical to our claim (lA.Sp . 



No Dimorphism for ^^^jz = ije> 

As suggested by the tabulated cases for jo = 2,3,4 p. [TO] there is also sufficient 
motivation to suppose that no dimorphism can occur quite generally for "^^^j^ = ±jc). This 
supposition can easily be validated by writing down the eigenvalue equations for angular 
momentum flIV.6ap - flIV.6bp in terms of the eigenvalues '^^''jz and jq: 



_ (ibl 



( Wj, + 1) cot ^ ■ gni^) = T Jo) sin^ ■ gs{^) 



dd 



+ + cot§ ■gs{^) = {^'^Jz± Jo)- 



g_Bm 



(A.lOa) 
(A.lOb) 



Now for ^^^jz 



jo this reduces to either 
dgni'd) 



dd 
dd 



-ijo + l)cot^-gR{^) = 
+ (jo + l)cot^-^5(^) = 2j, 



o 



gnW 

sin^ ^ 



(A.lla) 
(A.llb) 



or to 



dgnm 

dd 
dgsm 

dd 



- {jo + 1) cot t9 ■ gnii^) = 2jo sin^ i3 ■ gsi'd) 
+ {jo + l)coti3-gsii^) = 0, 



(A.12a) 
(A.12b) 



according to whether we take the upper or the lower signs in flA.lOap - flA.lObp . The 
dimorphism under consideration can occur if both systems flA.llap -f lA.llbp and f lA.12ap - 
f lA.12bp do admit normalizable solutions. This, however, can not be true. 
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Obviously, the critical equations are here (lA.lla|) and (]A.llb|) because these equations 



cannot admit simultaneously the required non-singular and normalizable solution. Indeed, 
the general solution of (lA.llaP is found as 

gR{^)=g. (sin?9r+^ . (A. 13) 

{g^ = const.) 

It is true, this is an absolutely reasonable solution, but if this is inserted on the right-hand 
side of the second equation flA.llbl) we get 

^ [(sin^9)^-«+i ■ gsid)] = 2ja (sm^Y^-' ■ gni^) (A.14) 

= 2g,jo {sin^y^^ , 

with the obvious solution being given by 

gsi^) = 2jo g. J-Jy^ I d^' (sin ^T'^ • (A.15) 



But this solution for the second angular function gs{^) is obviously singular for i) ^ tt 
which thus spoils the required normalization condition for the angular density ^''^ko('d) 
(IIV.24a|) ! Consequently, we have to reject the first system (]A.lla|) - (]A.llbp : and a regular 
solution can be expected exclusively from the second set of equations (lA.12ap -( ]A.12b|) . 



However, the rejected system (IA.llap -f lA.llb|l admits also an exceptional case, i. e 



jo = ^''^jz = 0. For this case, the second equation flA.llbp becomes homogeneous 

^^^ + cot{}-gs{^)=0, (A.16) 

and thus admits the trivial solution gs{^) = 0. Furthermore, the solution gR{'d) flA.lSp 
becomes specified to 

gR{d)=g,-^md , (A.17) 
and we end up with the angular density 

- i ■ S?^ - 1 ■ • (*-^«) 
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or in normalized form, resp. 

^%(^) = A • ^ • (A.19) 
Itt^ sm xr 

This is just the one single solution listed at the end of the subsequent table for jo = 4, 
whose general emergence for all values of jo has been proven by the above arguments 
below equation (lA.Sp . 

Thus, there remains to be considered the second system ( ]A.12ap - (lA.12b|) . Here, the 
second equation (lA.12b|) admits only the trivial solution gs{'&) = 0, and this leaves (lA.12aP 
for also as a homogeneous equation, namely just the precedent case (lA.llaP whose 

solution has already been specified by equation flA.13|) . Combining now both solutions 
gsi^) = and flA.lSp yields the angular density ^^^ko^-d) as 

An sm V An 

i. e. in normalized form 

{%(^) = 1 ^.^^^2ja~i (A.21) 

^ ' 27r2 (2jo - 1)!! ^ ' ^ ' 

Special cases {jo = 2, 3) of this result are found in the table on p. [TU for jo = 4 see the 
table below. 
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Jo 


Jz 




W^iii(^) 


fJ-A 

(^)eiii 




Wmiii 


4 


±4 




21 
40 


cos^ ^ - 1 




49 
2000 


147 
1000 


49 
400 


4 


±3 


A sin^ 


15 
32 


cos^ — I 




5 
256 


15 
128 


25 
256 


l^sin^^? (1 - f sin^i?) 


21 
160 


cos^ ^ — I 




49 


147 


49 
6400 


32000 


16000 


4 


±2 


^sin3^(l-|sin2^) 














^ (sin^ ^9 - 1 sin^ ^9 + § sin^ 


20 ^ 


cos^ -i? — 1] 


1 
500 


3 
250 


1 
100 


4 


±1 




9 " 
32 . 


cos^ ^ ~ |] 


9 

1280 


27 
640 


9 
256 


sin ^9 (1 - 4 sin^ + f sin^ - § sin*^ ^) 


51 
160 


cos^ — I 




289 
32000 


867 
16000 


289 
6400 


4 





1 . 1 

27r2 sini? 


§[■ 


20S^ t9 - |] 


1 

80 


3 

40 


1 

16 



Or 

to 

Angular Density ^^^ko{i9) i\IV.2a) 
and Angular Potential {''l4"'(i/) i\V.25a\) for jo = 4 



Appendix B 
Groundstate (no = 1 => jo = 0; ^^^3z = ^) 



The absence of the ortho-dimorphism for both '■''^j^ = and ^^^j^ = ijo imphes now 
that this effect of dimorphism can not occur especially for jo = and jo = ^- For the first 
case {jo = 0), the quantum number ^^^j^ can assume exclusively the value ^'^^jz = (because 
of ^jo ^ ^^^jz ^ jo)i and for ^^^j^ = we already have proven that the dimorphism cannot 
occur (see App.A). Nevertheless, this case jo = owns some peculiarities so that it may 
appear worthwhile to reconsider it in some detail. 

First, observe here that the angular eigenvalue system (IIV.6ap -( IIV.6bp becomes de- 
coupled because both constants io and io (lIV.15ap -f lIV.15bp must vanish together with 
^''^jz and jo- Thus, we are left with the following decoupled system 



-cot^- gni^) = (B.la) 
+ cot^- gs{^) = . (B.lb) 



The solution of the first equation (IB. lap is 

gR{^)=gR,.-sm^, (B.2) 
and similarly, the solution of the second equation (IB.lbp is obtained as 

gsW = gs,* ■ -\ (B.3) 

where gR^^ and gs,* are the integration constants. But actually, the coupling of these 
eigenvalue equations (lIV.6ap - (]IV.6bp does enforce the presence of only one integration 
constant. However, for the present exceptional situation (i.e. jo = ^''^jz = 0) we have to 
admit the existence of two independent integration constants and gs,* which become 
restricted by the normalization condition flIV.24ap to the constraint 

9l. + gh = - ■ (B-4) 
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This restriction admits the parametrization by only one constant ( 7*, say) 



9r,* = 
9s,* = 

so that the solution flB.2p -f lR3|) becomes 

9r{^) = 
9sm = 



— ■ cos 7* 

7r 



— ■ sm 7* , 

TT 



— cos 7* ■ sin 



2 sin 7^ 



(B.5a) 
(B.5b) 

(B.6a) 
(B.6b) 



71 sin^} 

with the angular density ^^^/cqI"*^) being found in agreement with the former claim flA.lQp . 
As a special case, one recovers here the former solution (1A.17|) for 7* = 0. 

Concerning now the energy of the groundstate of ortho-positronium {no = I, jo = 
0) ^^^jz = 0), one first observes that the angular density ^''^ko^'d) (IA.19P is the same as for 
the groundstate of para-positronium |2]; and therefore both the first angular potential 
correction flV.25ap and the constants (^^em, (''^/m f lVI.16ap . (''^mm flVI.22ap must 

then also be the same, i. e. 

cos^^-l (B.7a) 

(B.7b) 
(B.7c) 
(B.7d) 

Furthermore, these numerical results ensure again the validity of the quadrupole equation 
(1VIL4P which then leads us to the gauge field energy E^^ flVII.17p . This is to be added 
to the kinetic energy Ekin flVI.9ap - flVI.9bp for the groundstate {jo = 0) 





3 
8 




1 

80 


^''"Viii = 


3 

40 




1 

16 



kin 



2aB 

£km{l^ 



{2aB(3y 



4(2z/+ 1) 



(B.8) 



4(2z/ + l) 
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in order to obtain for the reduced form of the groundstate energy function, cf. (]VII.22|) 



2 



Eg''°^(z/) = -— (2Z/+1) 
as 



1 ^ 



(B.9) 



According to the principle of minimal energy, the (approximate) groundstate energy '"-^ 
of ortho-positronium is now given by the minimal value of the energy function 'E^''^\h') 
f lB.9|) and is found as 

Eg*'"^ ~ -7, 5378 [eV] . (B.IO) 

This groundstate energy is unique because the corresponding angular density ^''^ko{'d) 
turned out to be unique despite of the continuous set of solutions flB.6aP - flB.6bp . This 
confirms the former assertion that the ortho-dimorphism cannot occur for angular states 
due to ^''^jz = 0, i.e. especially for jo = 0. 

Furthermore, it happens that the ortho-groundstate energy is (in the 

quadrupole approximation) the same as the groundstate energy Ep'^ of para-positronium 
[2]. The reason for this is the simple fact that all ingredients of the groundstate energy 
function E^'v>(/3, u) are the same for ortho- and para-positronium, i. e. especially the 
angular-dependent functions ^%(^9) = ^P^ko{^), ^pW\^) = ^^'"(t?), ["Uo(r) = [PUo(r). 
The equality of these angular functions ensures then also the equahty of the coefficients 

^^^eiii = %u = ^ (B.lla) 
oU 

(^)/„, = %u = ^ (B.llb) 

(^'Wiii = = 4 • (B.llc) 

16 

And finally, the equality of the angular functions plus that of the associated coefficients 
admits then also the equality of the corresponding groundstate energies, i. e. 

Eg'^^ = e5,°'°^ = -7, 5378 [eV] , (B.12) 

cf. the process of partial extremalization flVII.18j) - flVII.22p . Indeed by means of the 
above mentioned equalities, both spectral functions S^^''^^{h') become identified for the 
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groundstate (no = np = ip + 1 = 1, ^^^j, = ^^^j, = 4 = 0) 



.{0,0} 



16 



(B.13) 



see the precedent paper [2] for the corresponding treatment of para-positronium. 
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Appendix C 



No Dimorphism for the first Excited State (no = 2^ 

Since the ortho-dimorphism can occur only if the quantum number ^''^j^ is in the range 
< \^^''jz\ < jo, that phenomenon of dimorphism can surely not emerge for the first 
excited state no = 2 jo = 0, 1; ^''^jz = 0, ±1). However, the mathematical mechanism 
for the exclusion is quite different. Whereas for ^^^jz = both configurations {1,-1} 
and {—1,1} of {io,^o} are allowed but yield one and the same physical state, in the 
case ^^^jz = 1 one has to reject the configuration {2,0} and to resort to {0,2}; and 
analogously for ^^^jz = — 1 one can conversely admit only the configuration {2,0} and 
has to reject {0, 2}. 

jo = 1, ("b'z = 

For {jo^^^^jz} = {1)0} one first deduces from equation (IIV.36P that Gsix) must be 
constant 

Gs{x) = (Tq = const. (C.l) 
Indeed, the second-order equation (]IV.29b|) is satisfied for this conclusion (IC.ip . The 



corresponding angular function gs{^) (IIV.28cP must then also be a constant 

9si^) = cTo , (C.2) 

and for the associated angular function gR{^) one deduces from the second first-order 
equation ( 1IV.6bp 



gji{^) = ^ sin^cos^ . (C.3) 

But both possibilities io = ±1 fllV.lSaP yield here merely a change of sign which may be 
compensated for by the integration constant ctq- Thus, one obtains for the (normalized) 
angular density ^^^kQ{'d) flIV.20p just the unambigous result flA.lQp . as claimed above. 
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jo = 1, ^'^j. = 1 



The question of unambiguity is somewhat different for ^^^j^ = ±1 (and still jo = !)• 
First, one concludes again from equation (1IV.35P that the angular eigenfunction gR{^) 
must look as follows 

gR{^)=P2-sm''^ . (C.4) 
If this solution is substituted in the first eigenvalue equation (]IV.6a|) . one finds 



ia.gsi^) = (C.5) 

so that either io or gs{'&) must vanish. Furthermore, the solution (lC.4p for gnii)) recasts 
the second eigenvalue equation (lIV.Gbp to the following form 



^^ + 2cot{} ■gsm=io-P2 (C.6) 
dxf 

whose solution is 

gsm = P2 —T-. - ^ — f ■ (C-^) 

[ sm V 2, sm v ) 

However, this solution must be rejected because it is singular at {} = n. As a consequence, 
we have to demand = 0, so that the equation (IC.6P for gs{i^) admits the trivial solution 
in place of the former gs{^) (10.7^ : 



gsm ^ . (C.8) 

This result says that the configuration = 2, = must be rejected which leaves us 
with io = 0, io = 2. Thus, both angular eigenfunctions do generate the unambiguous 
(normalized) density ^^^/co(''^) iii the form 

^%(^) = ^ (C.9) 
which validates again the above mentioned claim of unambiguity. 
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The remaining case ^'^^jz = —1 (and jo = 1) can be settled now by a very brief 
argument; namely by simply evoking the symmetry replacements flIV.47aP - flIV.47bp . This 
transcribes the present solution for ^^^jz = 1 to the corresponding unambiguous solution 
for = -1: 

io = 2 (C.lOa) 
io = (C.lOb) 
gRW=0 (C.lOc) 

gg(^) = J^ = const. (C.lOd) 



Indeed it is easy to see that this is a solution of both eigenvalue equations f lIV.6ap - fjIV.6bp 



with the angular density ^''^ko{'d) flIV.20p coinciding just with the former result flC.Qp for 
^^^jz = 1. Thus, both cases ^''^jz = and ^''^jz = ±1 for jo = 1 are actually not subjected 
to the ortho-dimorphism! 

Concerning now the energy E^'^"'^-'' of the excited states due to rio = 2, jo = 1, 
^''^jz = ±1, one can apply again the method of partial extremalization which is based 
upon the spectral function Sq\i') flVILlOp . For the present case {jo = 1, ''^^jz = ±1, 



no = 2) the spectral function adopts the following form 



(C.ll) 



2u+l 2u) 

This is the same spectral function as was found for the corresponding first excited state 
of para-positronium (f^ = ip = 1), see ref. [2], equation (IV. 161). Therefore the present 
unambiguous (!) ortho-state with quantum numbers {no = 2, Jo = I, = ±1} has 
the same RST binding energy as the para-state {n-p = 2, jp = i-p = 1, ^^^jz = iz = ±1}, 
namely [2] 

E^^'^^^ = -1,57527. . . [eV] . (C.12) 
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Appendix D 



Dimorphism vs. Elimination of jz- Degeneracy 

It is important to observe that the occurence of the ortho-dimorphism and the effect 
of degeneracy cUmination are two separate things which are independent of each other 
and have different origins. Nevertheless, they may occur in a combined way. For such 
a combined situation, there naturally arises the question of magnitude of those energy 
differences being induced by both kinds of level splitting (i. e. degeneracy elimination vs. 
dimorphism) . 

In order to elucidate this question a little bit more, it may be instructive to consider 
a somewhat larger value of the quantum number jo {jo — 4, say) where both effects 
can occur simultaneously. Here, the quantum number ^^^j^ can adopt nine values (i. e. 
— 4 < ^^''jz < 4) which, at first glance, would have to be hnked to five different angular 
densities ^^^fco(t?), namely any one density being associated with any one value of \'^^^jz\ (= 
0, 1, 2, 3, 4); and therefore one expects the occurence of five different energy levels due to 
jo — 4 (yielding the electrostatic fine structure). But the effect of dimorphism consists 
now in the additional splitting of three of these five energy levels (i.e. \^^^jz\ — 1, 2, 3) so 
that two of them (i.e. those due to j^^'^j^l = 0,4) remain unsplitted! Thus, the splitting 
effect says that there should ultimately be left over eight different energy levels for the 
quantum number jo = 4 (i.e. one level for any l^-^^j^l =0,4 and two levels for any 
\^^^jz\ — 1, 2, 3). The subsequent table presents a collection of the actual results. 
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jo = 4 



Jz 




9si^) 


^"%nf?^) flIV.20D 


4 







■sin'^^9 


3 






A ■ sin^ 7? 




\/ OOTT 


l^.sin^^ (l_|sin2^) 


2 


./f^-sin^^'^ (1-1 sin^^) 


O y "J'^ 


A-sin3?9 (1-1 sin2?9) 


^■sin^^ (1 - 1 sin2^9) 


■ sin^ cos "i? 

VOTT 


f sin^ ?9 (1 - f sin^ ?9 + § sin^ 


1 


6 . /l^ ■ sin^ i^cos^ (1-2 sin^ ^9) 


10 ./l^ sm^9 (1-1 sin^^?) 


A - sin?? (1 - f sin2'^ + 1 sin^??) 

TT^ V 9 3 / 


-2 ./i^- sin^T^cos^? (1-2 sin^^?) 

y TV \ / 


2./i^sm^9 (l-l sin^^) 


i§ - sin?? (1 -4 sin2 79 + f sin^^9- M sin^^?) 





■ sin ?9 (1-8 sin^ ?9 + 8 sin^ ^) 


■ COS?? (1-2 sin^^) 


1 . 1 


-1 


-10 sin^?? (1 - 1 sin^i?) 


6 cos?9 (1-2 sinH) 


^ - sin?? (1 - f sin2i9 + f sin^?9) 


-2y^-sm3t9 (1- 1 sin^t?) 


-2 y^-cos?? (1-2 sin^^) 


- sin?? (1 -4 sin2?9 + f sin^t?- § sin^?9) 


-2 


-|./|^-sm^^9cost9 


^■sin^ (1 - 1 sin^??) 






- sin?? (1 - 1 sin'^^) 


^ ■ sin=^?9 (1 - i sin^?? + § sin^?9) 


-3 


./if ■ sin^ ^ 


— -i/l^ ■ sin^ -(9 cos 
Y 0"" 


A ■ sin^ d 




T \J^- sin^ t9 cos ?9 


1^ ■ sin^t? (1 - i sin^??) 


-4 





■ sin^ ?9 





This table demonstrates clearly the symmetries of the ortho-spectrum: 

(i) the transition from ^^^j.^ to — '•^■'jz leaves invariant the angular density ^^^kQ{'d) OIV.20p . 
but not the angular functions gni^'d) and gs{i^) 

(ii) concerning the angular functions QrI^) and gs{^) themselves, the transition from 
^''^jz to —^^'^jz rearranges them according to the recipe (lIV.47aP - (lIV.47bp . apart even- 
tually from an overall change of sign (i.e. gR{i!^) =^ —gR{{}); gs{^) =^ ^gsi^))- 
Observe that the latter is admitted by the eigenvalue system (lIV.6ap - (]IV.6bl) 



(iii) observe also that the symmetry rearrangement (lIV.47aP - flIV.47bp generates a second 
solution for ^^^jz = 0, namely 

g^{^) = J—- sin^^^cos^? (1-2 sin^T?) (D.la) 

gg(^) = -.1. J— (l-8sm^^ + 8smH) , (D.lb) 
V TT sinw 

but this leads us to the same angular density ^^^kQ^i!)) and therefore yields no further 
level splitting. 

For a preliminary synopsis of those cooperative effects of degeneracy elimination plus 
dimorphism it is instructive to collect the results for the lowest energy levels (i. e. for 
principal quantum number < 5) in a further table, see below. Here it is especially 
interesting to compare the RST results of the spherically sjTumetric approximation, the 
results of the para/ortho dichotomy, and those of the dimorphism to the conventional 
results E^-' (ll.isp . Concerning the effect of degeneracy, the RST spherically symmetric 
approximation e!^^ {fourth column) respects the conventional degree of degeneracy (lI.lSp 

(n) 

but yields numerically only a moderate approximation to those conventional results 
{third column). 
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no,v = jo,v + 1 


Jz Jz 


4^^) [eV] 


E^"l [eV] 


(pWiii 


eJ,^'^ [eV] 


(^Vniii 


Eg^ [eV] 


fll.181) 


121 


[2] 


121 


flVI.22a) 


flVII.22D 


1 





-6,80290 . . . 


-7,2311... 


_i_ 

16 


-7,5378... 


j_ 

16 


-7, 5378 . . . 


2 


±1 


-1,70072... 


-1,5510... 


_1_ 

64 


-1,5752... 


J_ 

64 


-1,5752... 


3 


±1 


-0,75588 . . . 


-0,6692... 


_1_ 

64 


-0,68113... 


J_ 

64 


-0,6811... 


4 


±1 


-0,42518... 


-0,37369... 


9 


-0,3900... 


1 

64 


-0,3808... 


256 


9 
256 


-0,3900... 


5 


±1 


-0,27211 . . . 


-0,23909... 


289 


-0,2533... 


9 
256 


-0,2501... 


6400 


289 
6400 


-0,2533... 



All RST predictions in the quadrupole approximation are based upon the method of partial extremalizing, i. e. 



maximalization of the spectral function Sq^j, 



'pot 



(z/) + (^'^')m,„-^f,^(z/)l' 



1 1^ •^<7,p 

2u+l V 4 IaT 



1 



(D.2) 



Whenever both quantum numbers of para- and ortho-positronium do agree (i.e. jo = jp] ^^^jz = ^^^jz ''-'^-''miii = 
^Pj'miii), the corresponding energies are identical (E^-'^ = '^l^^)- The ortho-dimorphism arises because for the same 
quantum numbers jo = jp and ^^'^jz = ^^'^jz there can sometimes exist two values for ^^'>m\\\ (but not for *^P^miii), see 



the last two lines of the table, n = 4, 5. 



The functions epot{i^) (1VI.33P and |J,l^^^ (IVII.lSp . which enter the numerator of the 
spectral function (ID.2p . may also be written down in a closed analytic form: 



^pot(^) = ( 1 - ^.JM'^J^^ ^..2 ) ■ (2 ^ + 1)^' (D.3a) 



24^+2 [r {2u + 2)y 

Ie}, , 2z/2 + 5z/ + 3 10z/2 + 15i/ + 6 r(2i/ + 3/2) ,^ 

= uiAu^ - 1) - .(4.- - 1) ■ V^ -r(2. + 2) • ^^-^^^ 



Lowering of Groundstate Energy 

There is some curious effect with the groundstate (n = 1): whereas all the RST 
predictions e!^', E^-*^ for the groundstate energy {n = n-p = uq = 1) are lower 

than their conventional counterpart -E^onv (third column), the excited states {n > 1) are 
equipped by RST with higher energy. For instance, for jo = jv = 2, ^''^jz = ^^^jz = ±1 
the conventional prediction is Econv = -0,75588. . . [eV] whereas the RST prediction is 
Eg'^^ = Eg'^^ = -0,68113 [eV]. Thus one may expect that better trial configurations 
will further lower the RST energy of the excited states and thus shift it towards the 
conventional predictions. Of course, it is highly desirable to know the exact RST solutions 
in order to see more clearly how close the RST predictions do approach the conventional 
results. But those "better" trial configurations will then rather lower the groundstate 
prediction then raise it, quite in agreement with the true spirit of the principle of minimal 
energy! 

No Dimorphism for jo = 2 

The next effect to be discussed concerns the unexpected missing of the dimorphism 
for the level being specified by the quantum numbers no = 3, jo = 2, ^^^jz = ^^^jz = ±1- 
Indeed for this value (n = 3) of the principal quantum number no, the ortho-number 
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jo can adopt the values jo = 0,1,2 where it is clear that for the possibilities jo = 
^^'^jz = 0), jo = 1, and jo = = 2 the ortho-dimorphism cannot occur (see 

App.A+B). But for jo = 2 and \^'^^jz\ = 1 the dimorphism could principally occur; 
nevertheless it "incidentally" does not occur in the present quadrupole approximation! 
What is the reason for this? It seems near at hand that this is the very logical structure 
of the quadrupole approximation itself. 

In order to see this more clearly, reconsider the spectral function Sq\,{i') f lD.2p whose 
general form does apply to both ortho-((9) and para-(P) positronium. Since the RST 
energy predictions are obtained by extremalizing that spectral function Sq\,{i') with 
respect to the continuous variable the quantum numbers jo and being 
held fixed, the corresponding energy predictions can be different not only if at least one 
of the two quantum numbers jo and ^''^jz is different, but also if for the same pair jo, ^''^jz 
the parameter ^''\n\\\ (IVI.22aP can adopt more than one definite value. Recalling here the 
route back from ^^\n\\\ to the angular density ^''^kQ^i!)) via the relations f lA.4p . f lA.2p and 



l-3jdn cos2 ^ (t?) 1 ^ 

2 ' ^ ' ^ 

the two possible values for are seen to become induced by the two possible angular 

densities ^''^ko{'d) due to the same pair jo, ^''^jz- But the latter two possibihties for ^''^ko('d) 

are of course induced by the two possibilities for gR{^) and gs{^) (table on p. [70]) according 

to the prescription (IIV.20p . Especially for the present case jo = 2, ^''^jz = ±1 we find the 

following two angular densities 

1 



■ sin '(9 (D.5a) 

4t -sin?? ( 1 - - sin^?? 1 , (D.5b) 
n'' V 9 / 



see the tables on p. [7T]and[82l And the crucial point with these two densities flD.5aP -f lD.5bp 
is now that their substitution in the prescription flD.4p for ^^''P^m\\\ yields "incidentally" 
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one and the same mass-equivalent parameter ^^\n\\\: 

("Win = i- . (D.6) 
o4 

But clearly, when there is only one parameter ^^\n\\\ at hand for given values of jo and 
'^''^jz the associated spectral function Sq\i') is unique and provides us with one unique 
maximum, i. e. we end up with one single energy level. This is the reason why there is 
only one energy level associated with the quantum numbers Jo = 2, = ±1. 



Equality of Ortho- and Para-Levels 

Even when the ortho-dimorphism can occur (e.g. for the combinations jo = 3,4 and 
^''^jz = ±1, see the table on p. I163p . then the ortho- and para- levels are not necessarily 
different in the quadrupole approximation. For instance, take jo = = 3 together with 
^''^jz = ±1 and read from the table on p. I163l that one of the two ortho-parameters ('•''Vriin) 
of a dimorphic pair does agree with the corresponding para-parameter ^'''\n\\\, e.g. for 
jo = jv = 3: 

(D.7a) 

Jv = 3, ^^^Jz = ±1 ^ (^W,M = ^ , (D.7b) 

zoo 

or similarly for jo = j-p = 4:-. 

(D.8a) 

'"III 6400 

ocq 

jr = 4, = ±1 ^ (^Wm, = ^ . (D.8b) 

Consequently, one of the two spectral ortho-functions Sq\i') must agree with the cor- 
responding para-function S^^u) and therefore the corresponding ortho- and para- levels 
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must also agree, i. e. for jo = jv = 3, {no = n-p = 4): 



, Eg-'^ = -0,3808 . . . [eV] 
Eg''^ ^ <! (D.9) 

Eg'^^ = -0,3900 . . . [eV] = E^'^^ , 



and analogously for jo = = 4, {no = n-p = 5): 

I Eff'^^ = -0,2501 . . . [eV] 
Eg'^>^<^ ^ ^ (D.IO) 

[ t^^'^^ = -0,2533 . . . [eV] = eJ,^''^ , 

see Fig. VII. A. Thus, the result of this ortho-splitting is that the lower one of the two 
ortho-levels does agree with the corresponding para-level whereas the solitary ortho-level 
being left over is higher and therefore farer away from to the conventional prediction 
Ec^\ However, the width of this ortho-splitting is clearly smaller than the deviation of 
the ortho- and para-levels from their conventional counterpart: in (ID.Qp the dimorphic 



splitting width is (0,3900 — 0,3808) = 0,0092 [eV] whereas the average deviation of the 
ortho/para- levels from their corresponding conventional level is (0,425 ... - 0,385 . . .) [eV] 
~ 0, 04 [eV], which is (roughly) four times larger than the dimorphic splitting width! See 
also Fig. VII. A. This result says that the dimorphic splitting effect is of the same order 
of magnitude as the degeneracy-elimination effect which for para-positronium has been 
estimated by 0,42003 [eV] - 0,37369 [eV] ~ 0, 046 . . . [eV], see the table on p. 96 of ref. [2]. 

Thus, the difference of ortho- and para-energies conies about (in the quadrupole ap- 
proximation) solely through the dimorphic effect! 
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Appendix E 



Octupole Solution 

For the calculation of the octupole energy "Ef ^ flVIII.91aD - flVIII.91bD and its mass 
equivalent "M^^c"^ flVIII.91cP we obviously have first to determine the octupole potential 



(r) as the solution of the octupole equation f lVIII.83b|) . whereas for the amplitude 



field Cl{r) we will resort to our variational ansatz f lVI.lap - flVI.lb|) . For these purposes, it 



is very helpful to pass over to dimensionless variables; i. e. we put again 

y = 2(3r (E.la) 

.^,(,) . . ^ . (E,lb) 

i^'^>m2 2 pas 

so that the octupole equation flVIII.83bp adopts its dimensionless form as follows: 

A, ^^yp.iy) - 20 = V^-i e-y . (E.2) 

Furthermore, the octupole energy "E^'^^ reappears in terms of the dimensionless variables 



as 



^ aB ^ ' 11025 Vr(2z/ + 2)7 J \\ dy ' ^ ' 





2 



+20 



y 



Comparing now this result to the more concise one flVIII.95aj) one deduces thereof the 
"anisotropic" potential function 6\^\i') as 



oo 





By a similar calculation one finds the mass-equivalent function of equation flVIII.95b|) 



as 



oo 

= ^ [ dy -^'^^{y) . (E.5) 

^ ^ 11025 r(2z/ + 2)2 J ^ ^ 
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The claimed octupole identity in coefficient form ( 1VIII.96P is now easily verified by partial 
integration in flE.4p - (]E.5|) and use of the octupole equation ( ]E.2|) . 



Boundary Conditions 

Before turning to the elaboration of the desired exact solution ^''^^{y) of the octupole 
equation (]VIII.83b|) . or (lE.2p . resp., it is necessary to regard the right boundary conditions. 



The right-hand side of this octupole equation is the "octupole source" and does vanish 
very fast at infinity (r =^ oo, or y ^ oo, resp.), because we adopt the exponential 
trial function Cl{r) flVI.lap - flVI.lb|) . In contrast to this exponential decay at infinity, our 



octupole solution may be assumed to decay much slower (via some power law y'"^, w > 0). 
But also for sufficiently large value of the trial parameter u, the trial amplitude Q{r) may 
be assumed to be much smaller near the origin (y = 0) than the octupole solution ^^^^{y). 
Thus we have two regions of three-space where the solution ^^^^{y) is to be expected to 
obey the homogeneous octupole equation: 

A,^''hP4(2/)-20^^^^ = 0. (E.6) 

From here it is easy to deduce the following limiting behaviour in these asymptotic regions 

hm ^'^VM = (E.7a) 
\\m^'^,{y) = bI'^ . y\ (E.7b) 

where both and B^^^ are constants. 

The constant B^ characterizes the behaviour at infinity {y — oo) and may be de- 
termined by a first integration of the octupole equation (IE.2p . This can be realized as 
follows: First, rewrite the octupole equation flE.2p in the following form 

1 r.„.<l/l!^Ai (E,, 



y^ dy [ dy \ 
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Next, integrate here from y to oo with observation of the boundary conditions (lE.7aP - 
f lKTbj) and find 



.10 



d 

dy 



dy' y 



dy y 



(E.9) 



Now take this equation at infinity {y — oo) and observe also the boundary condition 
flE.TaP which then yields 



lim < y 



J/— >oo 



.10 



d 

dy 



| = -9i?W = _r(2z/ + 6). 



(E.IO) 



Thus, the asymptotic constant flE.7al) becomes hereby fixed to 



(E.ll) 



In order to try out whether perhaps the first integration step flE.9p contains also 
some information about the second boundary condition flE.7bl) one further carries out the 
integral on the right-hand side of that equation, i. e. 



dy' y 



dy y 



n=0 



r(2z/ + 7 + n) 



(E.12) 



By use of this result, the first step of integration flE.Dj) adopts the following shape 



d 
dy 



-r(2i^ + 6) e'^Y^ 



°° y2v-i+n 



n=0 



r(2z/ + 7 + n) 



(E.13) 



The conclusion from this result is that we must demand > 2 in order to have the bound- 
ary condition (]E.7b|) satisfied! But more detailed conclusions, concerning that boundary 
condition, cannot be drawn from the present result (lE.lSp . 
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Exact Octupole Solution 



In order to finally get the solution ^^^^{y) of the first-order equation (]E.13|) it suggests 
itself to try a power series of the following form: 

= r{2u + 6) e-^ Vp„ ■ y"^-'+- . (E.14) 

This ansatz satisfies the equation (IE.13P provided the coefficients Pn do obey the following 
recurrence formula 

^"^^ ~ 2u-2 + n ~ (2z/-2 + n) ■T{2u + 8 + n) ' 
with the lowest-order coefficient being given by 

= -(2.-3) .r(2. + r) ■ ^"^-''^ 



The solution of flKlSjl is {n = 1, 2, 3,4,.. .) 



_ 1 / 1 r(2z/ -3) 1 1 

^" ~ 9 \ r(2z/ + 7 + n)~ r{2u + 6) ' r{2u - 2 + n) j " ^ 

Thus, the (preliminary) solution appears ultimately as 

oo 

^'^4y) = r{2u + 6) e-yJ2pn ■ 2/'^^+^+" (E.18) 

n=0 

with the coefficients Pn (lE.16p - (lE.17p . However, observe now that this solution behaves 
in the vicinity of the origin {y = 0) like 

lim ^^^4(2/) = r(2z/ + 6) po y^"^^ ' (E.19) 

in contrast to the second boundary condition (]E.7bp (recall here the former demand 
1/ > 2). But it is possible to add some fourth-order term (~ y^) to the solution flE.lSp 



because the latter is a solution of the homogeneous version (IE.6P of the octupole equation 
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For the determination of the correct fourth-order term, to be added to the prehminary 
solution (]E.18|) . we simplify for a moment the problem by resorting to integer values of 
the variational parameter 2z/, i. e. 2v = 5, 6, 7, 8, . . .. For this simplified situation, the 
preliminary solution (lE.lSp can be transcribed to the following form: 



r(2z/ + 6) 



2u+5 



n=0 



n=0 



r 



(E.20) 



But in this form the deficiency of the preliminary solution (lE.lSp becomes now evident: 
for approaching infinity {y =^ 00) we have for integer 2z/ in flE.20p 

1 



hm i'yp,(y) = --Ti2u-3)-y' 
j/-5>oo y 



(E.21) 



which says that the preliminary solution flE.20p does not obey the required boundary 
condition flE.7aP at infinity! We can now easily remedy this deficiency of the preliminary 
solution by simply adding to it the negative of the false limit term (IE.2ip which itself is a 
solution of the homogeneous form flE.6p of the octupole equation flE.2p . Thus, the correct 
octupole solution for integer 2i/ can be read off from equation (]E.20p as 



{2u + 5)\ 



2u+5 



yO 



2i/-4 



r 



n=0 J n=0 

Obviously, this "integer- valued" solution ^^^i{y) consists of two parts: 



(E.22) 



^'^,{y) = ^'^t\y) + ^'^f\y) 



(E.23) 



with 

i^)vt\y) 



l(2. + 5)! Ji_e-'gV 



9 y^ 



n=0 



i2u + 5)\ 
9 7/5 



n=2u+6 



r 



(E.24a) 



2u-4: 



^f\y) = 1 (2. - 4). y' ^ ^ = U E 



n=0 



n=2u—'i 



r 



(E.24b) 
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such that the first part i^yp^°°^ (y) (lE.24aP is responsible for satisfying the boundary condi- 
tion (lE.7aP at infinity; and the second part {y) f|E.24bp is responsible for satisfying 
the other boundary condition flE.7bp near the origin. Both parts 



do cooperate in order that their sum ^'^^i{y) flE.231) can actually satisfy the octupole 
equation flE.2p . namely via the following differential relations: 

{^y - p) ^'^^^\y) = ^ {(2^ - 3) .y'-' - e-y .y'^] (E.25a) 

(a. - = -\ {{2u + 6) e-y -y"^-' - e"^ -y"^} . (E.25b) 

Finally, it remains to recast that "integer-valued" solution flE.22p to the general case 
of arbitrary (but real) values of the variational parameter v. This may simply be done by 
adding the negative of the false limit term flE.2ip to the more general (but preliminary) 
solution (lE.lSp which then ultimately yields the correct solution as 

^ oo 

^'^4(1/) = g r(2i. - 3) ■ + Y{2v + 6) e-yY^Vn y^''^^^'' . (E.26) 

n=0 

Mass-Equivalent Function fi\^\h') 

Once the exact octupole solution ^'^^P^^y) for arbitrary value of the variational param- 
eter u is now at hand, cf. flE.261) . one can go the last step and insert this solution in the 
former result flE.Sp for the mass-equivalent function This then fixes the latter 



object definitely as follows: 

Id,, 64 r(2z/ + 6) fr(2z/-3) 1 ^fPuM^, x\i „x 

ufUu) = ^ < ^ ^ + — — ■ > ^-^^r 4z/ + 3 + n)] } E.27 

^ ' 11025 r 2i/ + 22] 9 ^ 24^+3 ^ V 2" ^ ^ ^ ' I \ ^ ' 



where the coefficients Vniy) have already been specified by equations flE.16p - flE.17p . For 



integer values of 2z/, one substitutes the solution ^''^4(1/) flE.22p into the right-hand side 
of equation flE.Sp and thus finds the "integer-valued" version of f]E.27p as 

^f{y) = (-)4'^>(z/) + ^(^) (E.28) 
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with the first part jjl^^u) , being due to ^^^[°°\y) (lE.24ap . being given by 



^ ^ 99225 



64 (2^ + 5)! - 4)' - V (^^^^1±Z^^ 

64 (2z/ + 5)! ^ (2z/-4 + n)! 



(E.29) 



99225 f(2i^ + 1)!12 ^ n! 22'^-3+« 



and, analogously, the second part jj,^^ {v) , being due to ^'^^f\y) (lE.24bl) . by 



(o)„W^,^ ^ 1 (2^-4)! ^ (2// + 5 + n)! 



99225 [(2z/ + 1)!]2 ^ n! 22^+' 



^j2z^-4)^/ f (2^^+5 + n)! 

99225[(2z/+l)!]2 I ^ ^ ^' ^ ^i! 22^+6+n 



(E.30) 



The identities emerging in both equations flE.29p and flE.SOp may easily be validated by 
tracing them back to the following identity, needed here for z = 2v — A and ^ = 2z/ + 5: 

n=0 

Relative Magnitude of Quadrupole and Octupole Corrections 

In order to better understand the numerical differences between the quadrupole and 
octupole corrections (see the table on p. I163p it is helpful to oppose both perturbation 
orders / = 1 and / = 2 to each other on the same notational footing. In this sense, we 
recast the quadrupole potential A!2{r) to its dimensionless form ^^^2{y) by putting: 

^'^^2(0 = f^^^ (2/3«s ) ■ ^'^2{y) , (E.32) 



in close analogy to the octupole case (lE.lap - llE.lbp . By this rescaling, the quadrupole 



equation (1VII.4I) or (]VIII.83ap . resp., becomes simplified to 

(^A, - ^'yp.iy) = . (E.33) 

which says that the dimensionless quadrupole potential ^^^2{y) feels the same source (~ 
right-hande side) as does its octupole counter-part ^^^^{y), cf. flE.2p . 
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Indeed, there are many similarities between the former octupole and present quadrupole 
modes (presumably these similarities can be generalized to arbitrary / = 1, 2, 3, 4, . . .). For 
instance, the former octupole boundary conditions (lE.7ap - (]E.7bl) transcribe to the present 
quadrupole case as 



lim ^'''n^2{y) 



^{2} 
-tJoo 



(E.34a) 
(E.34b) 



The constant may again be determined by a first integration step for the quadrupole 
equation ( 1E.32p . quite analogously to the octupole situation ( 1E.9P : 



d 

dy 



dy y 



dy y 



2I/+3 . g-y 



(E.35) 



From this result, one concludes for the quadrupole limit at infinity flE.34al) 

= i ■ r(2z. + 4) , 



(E.36) 



quite analogously to the octupole limit flE.lip . 

The exact solution of the quadrupole equation f lE.32p can also be determined in form 
of a power-series expansion which is found to look as follows [2J 



1 °^ 

^''^2{y) = -n2u-l).y' + -e-yJ2K 



n=0 



5 ' ' " 5 

with the coefficients &„(z^) being given by 

r(2z/ + 4) T{2u-1 



i") ■ y 



2u+l+n 



(E.37) 



V(2u + h + n) V{2u + n) 



(E.38) 



Thus, the constant b}^^ (]E.34b|) is deduced hereof as 



5, 



{2} 



(E.39) 
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which thus turns out as the quadrupole analogue of the octupole hmit (]E.7b[) . And finally, 
the "integer-valued" {2u = 2, 3, 4, . . .) version of the general result (IE. 361) reads 



(2Z. + 3)! ^""^"'S,^ , (2^-2) 



2v-2 , 

2 -y y 
■y ^ l^—; 



(E.40) 



n=0 



This is obviously the quadrupole analogue of the former octupole result (IE.22P and consists 
again of two parts 

^'^2{y) = ^'W\y) + ^'^f\y) (e.41) 

with the first part ^^^'^\y) being again responsible for the behaviour (lE.34bp at infinity 



^'W\y) 



2v+Z 

1 - e"^ V ^ 
(2Z. + 3)! t,"-' 



5 y^ 
and the second part ^^^^^[y) near the origin (y = 0): 



(E.42) 



(2z/-2)! 



ra=0 



(E.43) 



Clearly, the sum ^^^2{y) f!E.40p must again obey the quadrupole equation flE.32p which, 
however, is ensured just by the separate differential relations for the long-range part 

and the short-range part ^^^f^?/) f lK42|l 



^ e-^ .y'^-' - i e-y ■y'^ 
5 5 



(E.44) 



y y2 



2^ + 4 _ 2u-l , ^ -w 2v 

^ — e ^ ■?/ +-ey-y 



(E.45) 



cf. the octupole case ( 1E.25ap -( 1E.25bp . 

The quadrupole energy "E^"^ = "eJ'^ + "Ef ^ (IV111.90al) - (IVIII.90bp and also its mass- 
equivalent "M^'^c^ f lVIII.90cp may now be re-expressed in terms of the dimensionless 
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quadrupole potential ^^^2{y) and thus read as follows 



"eJ'^^ = -— (2/3aB) (^'Wi)' ■ e^^\v) (E.46a) 
as 

= _^ (2/3aB) (Wmi)' ■ , (E.46b) 

as 



with the electrostatic function £2 i^^) being given by 



2' 



4 r(2z/ + 2)2 7 I V 

and the corresponding mass-equivalent function fi\^\i') by 

00 

A'^H^) = f YW^2 1 dy y'"^' e-y .^'"tm . (e.48) 



Just as in the octupole case (IVIII.96p . there is here also a corresponding quadrupole 
identity 

5^ ^ ^J'^^ (E.49) 

which may easily be verified by integrating by parts in flE.47p and using the quadrupole 
equation flE.33p . Observe here that these separate identities flVIII.96p and flE.47p can occur 
only because the separative method admitted us to decouple completely the quadrupole 
and octupole modes! 

Now that both modes are elaborated in great detail, one can proceed to the intended 
study of their numerical relationships. 

The qualitative features of these multipole solutions ^^^2{y) and ^^^^^{y) may perhaps 
be best seen from their sketch, cf. Fig. E.I below. For y ^ 20 one actually observes the 
asymptotic behaviour flE.34ap and f lE.7ap for ^^^2/A{y)- In a similar way, the behavior 
flKTbl) and flE.34bp near the origin {y = 0) is also neatly displayed by Fig. E.I. Further- 
more, one expects that any multipole solution ^''^21' (y) of higher order /' is smaller (for 
fixed value of u) than its counterpart ^^^2i{y) due to a lower multipole order Z (< /'). 
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For I' — 2,1 — 1 this expectation is also clearly expressed by Fig. E.I below. Observe, 
however, that the magnitudes of the multipole solutions ^''■^2,4 (y) are not directly respon- 
sible for the magnitudes of the corresponding multipole energy corrections. 
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Fig.E.I Quadrupole Solution ^^^2{y) H^E. S6\) and Octupole 
Solution ^^^4{y) (KM) foru = 3 
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Both multipole energy corrections "M.2'^^c^ (]E.46bp and "m\'^^c^ (]VIII.95b|) do essen- 
tially appear as the product of some reference energy jil^^u) and fi^^^u) times the mass- 
equivalent parameter {^^^miY or {^^^7712)'^, resp. Here, the reference energies fiyl^lu) flE.Sp 
and flE.47p do not depend on the quantum number ^^^jz of angular momentum but are de- 
termined solely by the spherically symmetric trial density $^(r) fIVI.lap - flVl.lbl) . Thus the 
"anisotropic" energy corrections are composed of two different effects, one of "radial" and 
the other of "angular" type. The latter one is measured by the mass-equivalent param- 
eters and ■"^^Vn'g which then play the part of quadrupole and octupole "strengths" 
and are displayed by the table on p. 11451 From the entries of that table it becomes 
clear that the octupole strength ^^^m2 is larger than its quadrupole counterpart 
(i.e. ^^^1712 > ^''^rrii) which is of course due to the specific anisotropy of the charge dis- 
tributions, see Fig. IV. A. But this does not yet admit the conclusion that the octupole 
corrections are larger than their quadrupole counterparts. 

Indeed, it is necessary to take into account also the second influence on the "anisotropic" 
energy corrections, i.e. the reference energies jJ-^fii^)- The subsequent Fig. E. II presents a 
plot of both functions fi^^^u) = | /if|'|-'^(z/) (IVII.lSp and jj^l^^u) (lE.27p . Observe here that, 
just as for yuf|'f"'^(z^) flD.3bp . there does exist also a closed analytic form for flE.271) : 

{e},._ 64 4z/^ + 28z/3 + 71z/2 + 77z/ + 30 

64 (12z/^ + 44z/3 + 81Z/2 + 67z/ + 20)r(2z/ + 3/2) ^ ' ^ 

~ 33075 z/ (2 z/ - 3) (z/ - 1) (4 z/2 - 1) r(2 u + 2) ' 
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Fig.E.II Quadrupole Function /i,2^^(z/) = |Atfn^(i^) ( 1 VJJ. 1 5l ) and 
Octupole Function ^E.27\) 
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